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COST-EFFECTIVE OPTIMIZATION OF RUBBLE-MOUND 
BREAKWATER CROSS SECTIONS 


PART I: INTRODUCTION 


Objectives 


1. The primary objective of this report is to introduce a systematic 
method by which planners and designers of rubble-mound breakwaters, specifi- 
cally those in the US Army Corps of Engineers (Corps) District offices, can 
formulate an optimum cross-section configuration and verify its effectiveness, 
both in terms of structural integrity and functional performance. Rubble- 
mound breakwaters, the most common coastal structures worldwide, are built to 
provide protection from direct wave attack to boat harbors (Figure 1) and to 
port facilities. Recent advances in coastal oceanography have greatly im- 
proved the understanding of wave generation, propagation, and transformation 
into shallow water. These advances, along with greater availability of mea- 
sured and hindeast wave data, have allowed procedures for design of rubble- 


mound structures to become much more complex than in previous years. The 


cee 


Figure 1. A rubble-mound breakwater protecting a boat harbor 


guidance available in the Shore Protection Manual (SPM) (1984) provides the 
basic tools for planning and designing breakwaters. This paper is intended to 
supplement that guidance by providing a practical perspective to the wide 
variety of environmental data now available to coastal engineers for rubble- 


mound breakwater design. 
Scope 


2. A brief review is presented of past and present criteria development 
procedures, design techniques, and related practical considerations, followed 
by a more detailed discussion of breakwater damage prediction and estimation 
of wave transmission characteristics. A systematic procedure is proposed to 
formulate alternative cross-section designs, evaluate their structural and 
functional effectiveness, and determine detailed dimensions which realize max- 


imum net incremental benefits. 
Definition and Purposes of Rubble-Mound Breakwaters 


3.  Breakwaters and, to some degree, jetties and groins are designed as 
barriers to sea waves, providing calmer water in their lees. Wave barriers 
can be constructed in many different ways, including vertical-sided concrete 
caissons, sheet-pile walls, wooden crib structures, and floating bodies. The 
oldest and most common type of wave barrier is the rubble-mound breakwater be- 
cause of its typical economy and constructibility in harsh coastal conditions. 
The long history of rubble-mound breakwaters has proven them quite reliable in 
a wide range of environments (Bruun 1985). A rubble-mound breakwater consists 
of sloped layers of stone or concrete shapes that are sized to withstand wave 
attack, excess settlement or loss of fill material, and to prevent scour, as 


shown in the typical cross section in Figure 2. Their inherent flexibility 
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tends to prevent catastrophic failure, even in the event of underdesign. The 
design parameters for rubble-mound breakwaters are rather inexact compared to 
those of most rigid civil engineering structures; thus, conservative overde- 
sign is quite common. 

4. Rubble-mound breakwaters can have a number of secondary purposes 
that are related to their primary purpose as a wave barrier. A breakwater 
protecting a harbor entrance and mooring area from wave attack might serve 
to divert currents and longshore transport of sediments. Also, it could be 
designed to provide access by people and equipment to the outer or deeper 
portions of the harbor. A breakwater protecting port facilities where cargo 
is being discharged and loaded might have these additional purposes and could 
even serve as a foundation for the port facilities themselves. This paper 
concentrates on considerations surrounding the wave barrier function. Fur- 
thermore, the perspective of the Corps as a public works agency is maintained 
since, in this case, the owners of the structure and the beneficiaries of its 
protection are the same (i.e. the taxpayers). The discussion to follow could 
also easily apply to a rubble-mound breakwater financed by private enterprise 
for commercial purposes, since tangible public benefits can, in many in- 
stances, be translated as profits. Many features of the planning and design 
procedures discussed later in this report can be extrapolated to planning 
and design of facilities other than rubble-mound breakwaters. The emphasis 
and most computational aspects will apply specifically to rubble-mound 


breakwaters intended as wave barriers. 


The Need for Optimization 


5. The construction cost for rubble-mound breakwaters is usually on the 
order of millions of dollars for smaller harbor or shore protection projects 
and on the order of tens of millions of dollars for larger harbor or port 
projects. The consequences of a dramatic structural failure include costs for 
repair of the breakwater which may approach the order of magnitude of the 
original construction costs due in part to expensive mobilization. Also, such 
consequences may include costs from property damage and inconvenience to port 
and harbor operations which occurred during the storm that damaged the break- 
water. These latter costs would typically be of a lower order of magnitude 


than the breakwater construction costs. All of these costs of rubble-mound 


breakwater failure are minimized by the tendency for this type of structure 
not to fail catastrophically. Catastrophic failure of flood control struc- 
tures (dams and levees) causes tremendous adverse consequences for the prop- 
erty and people in their flood plains, often including loss of life. The 
costs of these consequences can easily exceed the order of magnitude of the 
construction costs for the flood protection. This comparison illustrates 
that, in comparison to some other civil engineering works, a certain small 
risk of failure for rubble-mound breakwaters can be tolerated. 

6. Federal public works agencies in the United States have the statu- 
tory constraint for project authorization that the tangible benefits realized 
by the proposed plan must exceed all the life-cycle costs. This constraint 
has been further defined to apply to the incremental benefits and costs of 
each major feature of a proposed project. A rubble-mound breakwater built as 
a part of a federally funded project must "carry its own weight" in terms of 
its incremental net benefits. Recent administrative policies have provided 
additional restrictive criteria for federal financing of public works projects 
by requiring cost sharing with regional or local governments. These policies 
force planners to carefully consider the financeability of a project as well 
as its overall economic feasibility. Local sponsors of federally funded navi- 
gation projects commonly have severe limits on what costs they can share. A 
proposed breakwater project may be theoretically justified by a wide margin, 
but if it is not affordable it will not be built. Conversely, a sponsor may 
have the luxury of ample funding sources for cost sharing, but if a breakwater 
plan does not achieve enough incremental benefits, federal participation will 
not be possible. It is therefore critical that rubble-mound breakwaters be 
designed to provide the optimum trade-off between life-cycle costs and incre- 
mental benefits. This paper will deal with methods of formulating such an op- 
timum plan without extending planning schedules and budgets beyond reason. A 
commitment, both in time and money, is necessary, however, to address enough 


key questions for systematic optimization to be possible. 


Organization of the Report 


7. This introduction will be followed by a review of design principles 
for structural stability, including some of the many practical considerations 


involved in rubble-mound breakwater design. Current references offering more 


detailed discussions of various specific design considerations are given 
wherever possible, and readers are urged to consult these works. Review of 
design procedures is necessary in this paper to place an appropriate perspec- 
tive on simplifying assumptions made in this and other discussions of optimi- 
zation procedures. An introduction to a number of methods now in use to pre- 
dict damages to rubble-mound structures will be presented as tools to estimate 
future maintenance and repair costs for a breakwater design. Similarly, a 
discussion of methods to predict the wave transmission characteristics of 
breakwaters will follow to show how the structure's functional performance may 
be evaluated. The main paper will be concluded with a procedure to guide 
planners and designers of rubble-mound breakwaters from the choice of design 
criteria to determination of final dimensions. Appendixes will document the 


software available to accomplish some steps of this procedure. 


PART II: BASIC DESIGN PRINCIPLES 


Design Criteria 


8. There is a well-known tendency for subjective judgments to creep 
into supposedly systematic project planning endeavors in the earliest phases. 
A proven method to order your thinking in the conceptual phase of a project is 
to first thoroughly define the problems and opportunities at the site in terms 
of desirable goals to be achieved. This has long been the first step in the 
civil works planning process as practiced by the Corps. Two types of design 
criteria or "planning objectives," as stated in Corps planning guidance (Board 
of Engineers for Rivers and Harbors 1985 and Water Resources Council 1983), 
can be identified at this point relative to the function of a breakwater as a 
wave barrier. The first, and most familiar, is a criterion which defines the 
structure's ability to withstand the effects of extreme storms without itself 
suffering significant damages. This type of criterion can be referred to as 
the "structural integrity" or "survival" criterion. The second type, referred 
to as the "functional performance" criterion, deals with the effectiveness of 
the structure at its intended function which is to provide protection from 
waves. 

9. The structural integrity criterion determines the breakwater's life- 
cycle costs to the extent that a certain level of investment is necessary to 
prevent damages from an extreme event. There will always be a finite proba- 
bility that any storm, no matter how extreme, will be exceeded in intensity, 
so this criterion also determines the expected repair costs during the proj- 
ect's life. The most extreme sea state in which a particular breakwater de- 
sign will suffer no damages cannot, in practice, be precisely defined, as will 
be discussed later. The statement of a structural integrity criterion should 
be phrased with this in mind. It should be stated in terms of the desired 
effect, that is, prevention of breakwater damages (and associated repair 
costs). An example would be "damages to more than 5 percent of the breakwater 
armor will occur with less than 2 percent probability per year." There are, 
of course, numerous complications in achieving such a goal, including defini- 
tion of the types of possible damages and determination of the combined proba- 
bility per year of the physical parameters (wave height, wave period, wave 


direction, water level, storm duration, and others) which could cause them. 


Nevertheless, this is a workable statement in terms of an objective which is 
adaptable to more than one means of determining structural dimensions. 

10. The functional performance criterion determines the incremental 
economic benefits of a breakwater design since it defines the structure's 
level of effectiveness as a wave barrier. It also affects the cost since a 
certain additional increment of investment may be necessary to achieve a given 
level of effectiveness. This level of effectiveness can usually be stated in 
terms of a maximum transmitted wave condition during a given extreme event. 
The probability of exceedance for this event can in turn be related to 
property damage and other economic losses. Probability of exceedance is 
usually stated in terms of any single year, but it can also be stated in terms 
of all or some portion of the life of the project. A workable statement of a 
functional performance criterion might be that "10 percent of transmitted 
waves in any storm will exceed 1 m with less than 5 percent probability per 
year." This statement assumes that "10 percent of transmitted waves" can be 
related to some level of unacceptable property damage or operational 
disruption inside the breakwater. An even more general statement might be 
that "navigational delays and property damages from transmitted waves shall 
occur with less than 5 percent probability per year." 

11. Criteria of both types need to be defined for each section of the 
breakwater where either the environment (water depth, wave exposure, or other 
factors) or the required level of protection significantly differs. These 
sections can essentially be treated independently until a point when economy 
of breakwater materials, related constructibility constraints, and the transi- 
tion requirements become apparent. Usually the breakwater head, any elbows, 
and one particular section of trunk will take precedence over other sections. 
Head and trunk designs do not as yet lend themselves to reliable analytical 
methods and typically require much subjective judgment and extensive physical 
modeling. Most remarks in the rest of this paper will refer to the critical 
trunk section, with the understanding that other less critical trunk sections 


may have different design criteria. 
The Hudson Formula 


12. Investigations into the stability of rubble-mound coastal struc- 


tures were performed in the decade before the second World War by a Spanish 


10 


engineer named Cavanilles Iribarren. Iribarren (1938) presented the first 
widely used empirical formula for estimating a stable armor unit weight, given 
incident wave height, seaward slope of the structure, density of the sea wa- 
ter, and certain characteristics of the armor material. He assumed that 
stones on the outer slope were subject to gravity and wave forces, the latter 
of which included buoyant, impact, and friction components. The Iribarren 
formula was intended to predict the minimum weight stone which would remain in 
place when subject to waves of a given height. This height, defined in scale 
model tests as the level of "incipient damage," indicated that over the entire 
slope no more than 1-5 percent of the stones was displaced (d'Angremond 1975). 
The Iribarren formula is coming back into use in its original and in modified 
forms and will be discussed again later in this report. 

13. During and after World War II, the approach of Iribarren was con- 
tinued by Robert Hudson, a Corps investigator at the Waterways Experiment Sta- 
tion (WES) in Vicksburg, Mississippi. Hudson performed a great number of 
scale model tests on a variety of rubble-mound breakwater configurations. He 
also published a paper (Hudson 1958) which presented an armor unit weight pre- 
diction formula with many of the same features and assumptions as those of 
Iribarren. This formula is still in almost universal use by coastal engineers 
because of its relative simplicity and the many experimental and prototype 
tests of its reliability. The Hudson formula is 

p gH? 


3 (1) 
A Ky cot 6 


where 


W = weight of armor unit at the level of incipient damage* 


p,, = mass density of the armor material 
g = acceleration of gravity 

H = incident wave height 

AS SAO BD. 

p,, = mass density of the water 


Ky = an empirical stability coefficient 


@ = the angle from horizontal of the seaward slope of the structure 


* For convenience, symbols and abbreviations are listed in the Notation 
(Appendix E). 


14. Table 7-8 in the SPM (1984) presents the values for Ky recom- 
mended by the Corps for use in the Hudson formula. Values are presented for a 
variety of quarried materials and artificial concrete shapes. Each value is 
associated with a number of factors, including: 


a. Shape characteristics of the armor units (i.e., smooth, rough, 
round, or elongated rock). 


b. Position of units on the trunk or head of the breakwater. 


c. Wave form (i.e., whether or not the wave is breaking directly 
on the structure). 


d. Slope or range of slopes (in some cases). 


e. Method of placement (random versus special individual 
placement). 


f. Number of layers of armor units to be placed on the slope. 
g. Relative gradation and smoothness (for quarried rock). 

15. An important point to note about the Ky, values in the SPM (1984) 
is that 58 percent of them were derived from monochromatic wave model test re- 
sults, while the rest are interpolated values. Another factor of importance 
is that some of the armor unit types for which Ky values are presented have 
actually been used in only a small number of prototype breakwaters. All of 
the units lack systematically documented prototype verification of their rela- 
tive stability, though efforts are currently under way to consolidate histori- 
cal performance of Corps constructed breakwaters. Uniform rough angular quar- 
rystone, riprap (graded rough angular quarrystone), and dolosse have been most 
extensively tested in scale models and currently have the best documentation 
of prototype experience (Jackson 1968a and Carver 1983). 

16. The coefficient Ky» as applied in the Hudson formula with its 
basis in the assumptions of Iribarren, does not directly account for as many 
as 20 or more design conditions (Ligteringen and Heijdra 1984) that are now 
known (in at least a qualitative sense) to affect breakwater stability. Some 
investigators (Brorsen, Burcharth, and Larsen 1974 and Burcharth 1979) have 
questioned whether the Hudson formula is reliable for predicting stability of 
dolosse and other slender concrete armor units. In the future, these units 
may require variable Ky factors related to slope and other conditions not 
now inherent in the values presented in Table 7-8 of the SPM (1984). Some of 
the other conditions of concern include: 


a. Influence of wave period or the steepness of individual waves 
(Ahrens and McCartney 1975 and Losada and Gimenez-Curto 1979). 


Influence of wave groupiness in natural irregular seas (Bur- 
charth 1979) 


ec. Effect of the foreshore or the breakwater toe on wave transfor- 
mation (Bruun 1979 and Kjelstrup 1979) 


d. Effect of oblique waves (Losada and Gimenez-Curto 1982 and 
Christensen et al. 1984). 


e. Interaction of waves with monolithic crest elements or densely 
packed underlayers, such as resonance of reflected waves with 
incident waves (Jensen 1983). 


In 


f. Friction of outer armor material with underlayers (Hedges 1984) 


g. Mechanical strength (resistance to tension, compression, im- 
pact, fatigue, etc.) of individual armor units (Poole et al. 
1984 and Groeneveld, Mol, and Zwelsloot 1983). 


h. Potential settlement, foundation failure, and related geotech- 
nical problems (Thorpe 1984). 


i. Seismic stability. 


17. The Hudson formula can be applied to interpret scale tests of pro- 
posed designs to measure the "actual" Ky of an armor unit in a particular 
breakwater configuration, in which case many of the above factors would be ad- 
dressed. A series of successive tests on the same configuration with varying 
monochromatic wave period can determine the critical period when waves of that 
height would break directly on the face of the armor slope. Likewise, this 
"sensitivity analysis" approach (vary one parameter while holding others con- 
stant) can provide estimates for the reliability of the point of incipient 
damage and damage rates for more severe wave height and period combinations. 
Tests with irregular waves are also possible and should be considered, even 
though the procedures involved and interpretation of results in terms of Hud- 
son formula parameters are less standardized. Physical modeling is an essen- 
tial step in the cost-effective design of rubble-mound breakwaters and should 
not be neglected for any except the smallest, most inconsequential structures 
(Paape and Ligteringen 1980). Some specific techniques for verifying armor 
stability and damage rates by scale model testing will be discussed later in 


this report. 


Alternative Stability Relations 


18. The Iribarren formula, as mentioned earlier, has recently been re- 


ceiving renewed attention worldwide because of some spectacular failures of 
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large rubble-mound breakwaters in the last 10 years (Stickland 1983). The 
Iribarren formula in its original form appears as follows (d'Angremond 1975): 


No ust? 
W = (2) 


A(u cos @ + sin @)3 


where 


N = empirical coefficient related to the armor material character- 
istics (comparable to K, in Equation 1) 


u = coefficient of static friction between individual armor units 
(equivalent to the tangent of the angle 6 at which armor would 
slide from gravity alone; values found by Graveson, Jensen, and 
Sorensen (1980) are presented in Table 1) 


Table 1 


Values for the Coefficient of Static Friction u 


_Type of Armor ~~ ~+Coefficient u =~=~=~=~~—~—~— Angle of Repose 90 
Round seastones 120 45 
Quarrystones eet 48 
Concrete cubes tee 50 
Concrete tetrapods ~1.4 ~55 
Concrete dolosse 2a if ~70 


19. The original Iribarren formula has only one additional parameter u 
with N being essentially equivalent to Ky . This additional explicit pa- 
rameter appears to have little advantage to offer, except that static friction 
has recently been investigated as a potentially critical factor in the overall 
stability of complex artificial shapes such as dolosse (Price 1979). It al- 
lows the Iribarren formula to account for the marginal stability of materials 
placed at their natural angle of repose. This factor might be used also in 
the future as a measure of seismic stability of rubble-mound structures. 

20. Engineers at the Danish Hydraulics Institute (DHI) have proposed a 
modification to the Iribarren formula for application with scale model tests 
using irregular waves (Graveson, Jensen, and Sorensen 1980). This DHI- 
Iribarren formula is 

3he 


po gy H_L 
Z Sapa a ate 
W fae (3) 


K oC cos @ - sin a)? 


where 


H, = significant wave height of the incident irregular waves 
(average height of highest one-third waves at site) 
Ly = wave length at the site corresponding to the period of peak 


energy density for the incident irregular waves 


K, = alternate stability coefficient = Lp /NH, 


21. The principal modification of the original Iribarren formula is the 
substitution of an alternate stability coefficient on the basis that the orig- 
inal stability coefficient (N in the numerator of Equation 2) is a function 


of the wave steepness H,/L A number of other investigators have proposed 


p° 
similar stability relations (Rybtchevsky 1964, Jensen 1984, and Ahrens 1984). 
A similar modification to the Hudson formula could be made by substituting Ky 


= KqH,/Lp . Ahrens (1984) found that stability by "reef type breakwaters," or 


low-crested breakwaters without traditional multi-layered cross sections (ba- 
sically homogeneous rubble-mounds), was reflected with greater confidence us- 
2 
stp 

original Hudson formula (Figures 3 and 4). 


ing a modified Hudson formula with H in the numerator than with the 

22. Engineers at Delft Hydraulics Laboratory in The Netherlands re- 
cently performed an extensive series of scale model tests of the stability of 
rock slopes under random wave attack (Van der Meer and Pilarezyk 1984). These 
tests resulted in the formulation of a set of stability formulae for quarry- 
stone armor of breakwaters and revetments. Their tests also gave information 
on how to predict damage rates as a function of the number of incident waves. 
Armor layer gradation was found to have a lesser effect than that found by 
other investigators (Ahrens and McCartney 1975). Slope angle was found to 
have an effect on stability similar to that predicted by the Hudson formula. 
Wave period effect was investigated as a function of the "Iribarren number" or 


surf parameter as follows: 


ae tan 6 
a 2 
Hy (4) 
Ie 
fe) 


where Lo e gTe/2n, based on the average wave period T, 
23. The influence of wave period was found to correspond roughly with 
the traditional distinction between breaking and nonbreaking waves. The ef- 


fect of variations in the incident wave spectral shape, as measured in various 
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Figure 4. Model data plotted by spectral stability number 
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ways to reflect both irregularity and groupiness, was found to be minimal. 
This result differs from the conclusions of other tests relative to the 
influence of wave groupiness on stability (Burcharth 1979). A major influ- 
ence by core permeability was found. The stability formulae proposed for 
rubble-mound (quarrystone) structures with permeable cores (D5 armor /De¢ core 


= 3.2 , as tested) for breaking waves (&— < 2.5 - 3.5) was 


Os22 
H, 33 -0.54 
AD = 5\58 73 (5) 
n50 N 
or, equivalently, 
H, p Whe S, 1/3 
n50 n50 N 


The formula proposed for nonbreaking waves (& > 2.5 - 3.5) with cot @ < 3 
was 


H S 


1/2 2 Oje 
= 1.65 (cot 6) — E (7) 
y/e 


whereas for nonbreaking waves (& > 2.5 - 3.5) with cot @ > 3 the formula 


was 
ils »2 ae 0.1 
AD = 2.86 72 E (8) 
n50 N 
where 
H. = the significant wave height of the incident spectrum 


Dn50 = the nominal diameter, based on the mass of the 50th percentile 
W50 from the armor material mass distribution curve 


| 
= (Weo/0,)'/3 


So =a dimensionless damage level, defined as the number of 
equivalent Dn50 cubes eroded over a width of Dn50 


= 2-3 for incipient damage (as with the Hudson formula) 


= 8 to 17 for armor layer "failure" (significant exposure or 
underlayers) 


N = number of incident waves 
The range of & values from 2.5 to 3.5 for the transition from breaking 
to nonbreaking wave conditions apparently represents the difficulty in de- 


scribing an irregular sea state as either breaking or nonbreaking, since both 
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breaking and nonbreaking waves can occur in the same sea state. Others who 
have investigated breakwater stability as a function of the surf parameter 
(Equation 4) include Gunbak (1976) and Losada and Gimenez-Curto (1980). 

24. A relatively complete list of rubble-mound breakwater stability 
formulae, proposed by various investigators over the years, was published by 
the Permanent International Association of Navigation Congresses (PIANC) 
(1976). The variety of model tests and prototype experience inherent in these 
formulae and those developed since 1976 is but a small fraction of the many 
thousands of monochromatic wave tests conducted to determine Hudson formula 
parameters used to design hundreds of breakwaters all over the world. Use of 
these other stability relations should, therefore, be applied only in conjunc- 
tion with traditional procedures using the Hudson formula for comparison. A 
conservative choice can then be made between the stable armor weights and 


damage rates predicted by the Hudson formula and these alternate methods. 


Practical Considerations for Stability 


25. The analytical methods available for predicting rubble-mound break- 
water stability have been shown not to include many important considerations 
that could cause a structure to fail. Breakwater design has always involved a 
great deal of subjective judgment and probably always will. Some of the most 
pertinent practical considerations that must be made in determining rubble- 
mound breakwater material characteristics and dimensions are reviewed below. 
Comprehensive review of both practical and analytical considerations is avail- 
able in the SPM (1984), Angerschou et al. (1983), Institute of Civil Engineers 
(1984), Jensen (1984), and Bruun (1985). 

Incident wave conditions 

26. The incident wave conditions are traditionally defined as the wave 
height at the seaward face of the structure with a further distinction as to 
whether or not the waves are breaking. This breaking versus nonbreaking cri- 
terion has been argued extensively over the years. The convention remains in 
practice, however, due to obvious differences in design conditions for rubble- 
mound structures built in shallow water, where wave heights are depth limited, 
and in deeper water (depth > ~15 m), where waves have not transformed to the 
point of breaking in front of the structure. The natural irregularity of sea 


states can be fairly well represented by a2 single height and period related to 
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some specified exceedance value, but it must be acknowledged that both height 
and period will vary in any storm. Consequently, some incident waves will be 
breaking on the structure, and others will not. The succession of high and 
low waves (wave groupiness) and of breaking and nonbreaking waves can be a 
critical factor. The potential effects of wave groupiness or multiple con- 
verging wave trains (multi-peaked spectra) are difficult to assess without a 
substantial amount of field data and scale model testing with irregular waves. 

27. The alternative stability coefficients for the Hudson and Iribarren 
formulae discussed above which include wave steepness H/L provide one means 
of making a more explicit description of incident wave conditions. Other de- 
scriptive parameters that have been investigated include the surf similarity 
parameter in Equation 4 (Bruun and Gunbak 1978, Burcharth 1979, Losada and 
Gimenez-Curto 1980, Van der Meer and Pilareczyk 1984, and Bruun 1985) and the 
Stokes or Ursell parameter HL®/h3 (Carver 1983). Estimated values of these 
wave form parameters can be used as a more systematic means of classifying 
individual waves as breaking in the critical plunging mode, as spilling, or as 
nonbreaking. Irregular sea states require further definition in terms of ei- 
ther time domain characteristics or spectral (frequency domain) parameters. A 
number of useful parameters for characterizing irregular waves are discussed 
by Rye (1977). 

28. Wave transformation effects caused by the proposed construction 
works themselves cannot be neglected. Breakwaters with shallow slopes or with 
extensive toe development can also change the wave conditions at the waterline 
on the seaward face by "tripping" the waves. Scale-model tests are necessary 
to quantify these effects on the armor layer (Jackson 1968b). Relatively 
steep and impermeable structures may partially reflect incident waves such 
that resonance of incident and reflected waves causes scour near the toe. De- 
termination of the sensitivity of a structure to these effects from oblique 
waves requires scale-model testing in a three-dimensional wave basin. These 
potential problems make physical modeling critical for reliable estimation of 
the stability of a proposed rubble-mound breakwater. 

29. The duration of a storm at sea is a real world parameter that 
should be considered in any design effort or laboratory stability analysis. 
Figure 5 illustrates the time-history of significant wave height, peak spec- 
tral wave period, and predominant direction of wave propagation for a storm 


in the Gulf of Alaska simulated from synoptic weather data at 6-hr intervals. 
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This rise to and fall from peak conditions over many hours, sometimes days, is 
typical for severe storms in most areas of the world. The peak condition is 
typically applied in extremal analyses, but the duration of conditions above a 
threshold related to the stability of a proposed structure is also important. 
Simulation of many hours (or many thousands of waves) is performed as standard 
practice for breakwater stability tests by a number of prominent laboratories 
(Owen and Allsop 1984 and Van der Meer and Pilarcyzk 1984). The effect of 
duration on breakwater stability is discussed by Graveson, Jensen, and 
Sorensen (1980), Jensen (1984) and Bruun (1985). 
Foundation considerations 

30. The weight of a rubble-mound breakwater and the hydraulic effects 
it causes near its foundation are potential factors which can lead to a struc- 
tural failure. Investigation of gravity related stability problems, such as 
slip failure of the foundation or excessive (possibly differential) settle- 
ment, requires the attention of a geotechnical specialist. Hydraulic problems 
such as scour at the toe must be addressed in the earliest stages of design. 
The suitability of a natural foundation and the possibilities for preventive 
measures can ultimately determine the feasibility of constructing an entire 
breakwater. Excavation of poor foundation materials and replacement with fill 
or artificial improvement of the strength of natural materials can amount to a 
substantial fraction of the project cost. The need to place filter materials 
or other scour protection along a breakwater can also substantially constrain 
the geometry of the armor and underlayers. Seismic stability analyses in 
areas subject to earthquakes should be performed. All of these geotechnical 
considerations require extensive field data consisting of numerous borings 
supplemented by acoustic surveys and penetrometer tests. 
Primary armor 

31. During the past 40 years many lengthy journal articles, textbook 
chapters, and conference papers have been written on the subject of armor de- 
sign for rubble-mound breakwaters. A discussion of the entire multitude of 
practical considerations applicable to armor design would be beyond the scope 
of this report. A comprehensive review is available by Baird and Hall (1984) 
in which many of the most important factors in armor design are discussed. 
Rubble-mound breakwaters have a tendency to be designed from the top down be- 
cause the exigencies of design and construction of those portions exposed to 


direct wave attack tend to constrain all other features. The stability 
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formulae, presented as Equations 1 through 8, apply only to the resistance 

to displacement of individual armor units. The use of concrete armor units 
also requires the investigation of mechanical strength related to the 
interaction betw rn the units in the armor layer and the associated impacts, 
fatigue, and creep (static) effects that occur. Quarrystone can be subject 
also to fracturing without displacement, but experience shows rock and the 
bulkier conerete units (such as plain or modified cubes) develop less of this 
sort of damage than do more slender concrete units (such as dolosse). A num- 
ber of proof tests and other quality control procedures have been proposed to 
account for mechanical strength limitations in concrete armor units (Burcharth 
1981 and Price 1979) which should be considered for application in any project 
involving these units. Large concrete armor units should be designed with the 
advice of a specialist in concrete engineering, particularly where fiber re- 
inforcement is contemplated. The availability of existing forms should be 
investigated before fabrication of expensive specialized concrete forms is 
undertaken (Owen 1985). Design formulae indicate a minimum size armor unit, 
but the availability of existing forms and other practical factors may make 
slightly larger units more economical. 

32. Design considerations related to the geometry of the armor layer 
are of particular interest in discussions of optimization since armor units 
are typically the most expensive materials used in a rubble-mound breakwater. 
The extent to which primary armor extends below the still-water level on the 
seaward face is typically set subjectively at 1.5 to 2.0 design wave heights 
(Figure 1). A berm of secondary armor or underlayer material at the toe of 
the primary armor is considered good practice, enhancing both the accuracy of 
underwater placement of the primary armor units and their resistance to slid- 
ing failure near the toe. The primary armor is usually extended below the 
waterline on the leeward side by 0.5 to 1.5 wave heights, depending on the de- 
gree of overtopping anticipated. If a monolithic wave screen is planned for 
construction on the crest (as illustrated in Figure 6) and virtually no over- 
topping is to be allowed, armor on the lee side need only be sized to remain 
stable in the ambient wave climate on that side of the breakwater. Wave 
sereens and monolithic crest structures are sensitive and highly specialized 
features (Jensen 1983) and will not be dealt with in this paper. 

33. The allowances above, along with the crest width, crest height, 


and number of armor units comprising the thickness of the primary armor layer, 
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CREST ELEMENT 


Pots 


Figure 6. Typical breakwater cross section with a monolithic crest element 


determine the total volume of primary armor per unit length of structure. 
The conditions determining these dimensions will change along the length of a 
breakwater. Transitions should be gradual with conservative allowances for 
the limited confidence in the predicted variations in design conditions. All 
of these considerations must account for both extreme high water conditions 
and the possibility of a low tide condition which could greatly complicate the 
stability of features near the toe. 

34. The dimensions of the armor layer generally are formulated as func- 
tions of the primary armor weight. The armor thickness and crest width are 


related to the weight of the armor unit by the following relation: 


1/3 | 
ES We 
r= nk, (4) (9) 


where 
r = total average layer thickness or crest width 
n = number of armor units comprising the thickness or width 
(usually 2 for the thickness and 3 for the crest width) 
K, = "layer coefficient" (see Table 7-13, SPM 1984), an empirical 


measure of the thickness compared to tnat of the same number of 
equivalent cubes 


35. The weight of the individual armor units, as determined by the 
Hudson formula (Equation 1), is a function of the slope, the armor material's 
density, the K, factor, and the wave height cubed. A small increase in de- 
sign wave height makes a substantial difference in the armor weight, i.e., a 
10 percent increase in H corresponds to a 33 percent increase in W . The 
armor thickness will increase only 10 percent. The in-place unit price of 
armor material (both quarrystone and concrete) will vary directly with the 


total weight of the units relating also to the practicalities of quarry 
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development, concrete unit forming, and difficulties in handling. A reduced 
slope (increased cot @ in Equation 1) will reduce the armor weight require- 
ment but will change the runup characteristics of the seaward face in a non- 
linear manner. Overtopping and the associated transmitted wave characteris- 
tics are then affected, which in turn affects the required crest elevation for 
acceptable wave attenuation. The overall volume of a roughly trapezoidal- 
shaped breakwater (in cross section) increases as the square of the increase 
in the crest elevation. A significant effort is therefore necessary to deter- 
mine the most economical combination of slope, armor type, armor weight, and 
crest elevation for every pair of functional and structural design criteria, 
even when first cost is the only consideration. 

Other breakwater features 

36. The constraints involved in primary armor layer design can 
sometimes overshadow other considerations for design of the secondary armor 
layers, underlayers, core, foundation filters, and scour protection (Fig- 
ure 2). The terminology of the SPM (1984) refers to a secondary armor layer 
as material placed on the face of the breakwater below the primary armor 
layer. An underlayer is placed between the armor on the exposed face and the 
core in the interior of the structure. Underlayers serve basically three 
functions: to keep the core in place through filtering action, to further 
dissipate wave energy that has penetrated through the primary armor, and to 
act as a foundation for the primary armor. These functions also apply to 
underlayers between the primary armor and the natural foundation (sea floor). 
Multiple underlayers may be required to satisfactorily accomplish all these 
functions. Material with small enough voids to hold finer core material in 
place may be too fine to stay in place itself under the larger voids in the 
primary armor layer. The primary armor also needs a relatively rough surface 
under it to discourage sliding. A coarser underlayer also provides some pro- 
tection from waves during placement of the primary armor (Hedges 1984). 

37. Filtering criteria developed for water quality or seepage control 
purposes, such as Dis (filter) < 5Dgc (foundation) (Sowers and Sowers 1970), 
can be restrictive, thus the size gradation should be a consideration in evalua- 
tion of borrow areas for core material. Efficient use of quarry materials is 
encouraged in the SPM (1984). Given the practical problems of accurate place- 
ment of complex underlayers in the field (especially underwater), this goal 


may not always prove as economical as relaxing gradations of the various 
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layers such that their number, complexity, and associated construction quality 
control requirements are minimized. Unfortunately, breakwater specialists do 
not agree on a precise filter criterion for rubble-mound breakwater underlay- 
ers (Jensen, Graveson, and Kirkegaard 1983), and physical modeling of scour of 
core material is complicated by scale effects (Hedges 1984). 

38. A densely packed core can reflect a significant amount of wave en- 
ergy back through the underlayers and reduce the stability of the armor or in- 
crease scour near the toe of the breakwater. A core and underlayer system 
that reduces wave energy through turbulence and frictional loss is preferred 
to a more reflective system. A core that is too permeable can transmit waves 
as much as 80 percent of the incident wave height (Kogami 1978), and it may 
pass littoral materials. Some useful experiments with wave transmission 
through porous rubble-mound breakwaters were performed by Madsen and White 
(1976) and continued by Seelig (1980a). Their methods are helpful in predict- 
ing wave transmission characteristics and will be discussed again later in 
this report. The effects of variations in permeability are discussed further 
in Bruun (1985). 

Head and elbow construction 

39. The inevitable lateral flow across round heads and elbows and the 
reduced interlocking and compaction in these areas complicate just about every 
facet of breakwater design. Practical methods to deal with these compli- 
cations consist primarily of conservative adjustments to analyses as applied 
to sections of the breakwater trunk. This type of adjustment has limited 
confidence as evidenced by the frequent need to repair heads and elbows of 
conventionally designed rubble-mound breakwaters. Model testing in a three- 
dimensional wave basin is at present the only reliable means of improving this 
confidence. This is particularly important with slender concrete units (such 
as dolosse), which may have little or no increased stability over rock or 
bulky units in lateral flows (Burcharth and Thompson 1982). It is this fact 
that has caused some investigators to question the reliability of the Hudson 
formula and the associated K, factors published in the SPM (1984) for use in 
head or elbow design (Angerschou et al. 1983). The detail design of heads and 
elbows will very likely remain a highly subjective and empirical process for 
some time. 

Toe construction 


40. A number of practical problems related to the toe of rubble-mound 
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breakwaters have already been mentioned. This area of transition from a 
hopefully stable static environment (the breakwater) to the natural, often 
dynamic, sea floor is critical to the overall stability of the structure. Toe 
features not only protect the bottom from scouring (which can lead to under- 
mining) but also support the weight of the armor material above. The need to 
provide primary armor 1.5 to 2.0 wave heights below the still-water level can 
conflict with the need to filter foundation sediments at the tow. This is 
particularly true in high tidal ranges where low water conditions can expose 
the toe to more extreme wave effects. The support of armor materials is most 
reliably accomplished with a substantial berm of secondary armor or underlayer 
material at the toe of the armor slope. This berm should have at least sev- 
eral units or a minimum 3-m top width. Wide differences in the size of the 
bottom sediments and the breakwater material near the bottom may require 
excavation of a trench along the toe to accommodate a toe berm with an ade- 
quate filtering underlayer, as illustrated in Figure 7. Geotextiles can be 
used also in some instances to reduce the height of toe features and the as- 
sociated exposure to more severe wave energy. The concurrent physical mod- 
eling of armor stability and toe scour is complicated by scale effects, but 
model tests can reveal trends which could suggest a compromise of either the 
filtering criteria or the extent of primary armor. One radical concept in toe 
design is the "wave reducing berm" (Delft Hydraulics Laboratory 1983) which 
provides artifically shallow depths for dissipation of wave energy. Sugges- 
tions for design of more conventional toe features are discussed by Eckert 
(1983) and Jensen (1984). 
Construction equipment and techniques 

41. The constructibility of a rubble-mound breakwater design is an 
extremely important and practical consideration that can control its over- 
all feasibility. Smaller breakwaters can often be constructed with con- 
ventional land-based construction equipment and techniques by building from 
the shore outward. Detached breakwaters can be constructed in this fashion 
only if a temporary causeway to the permanent portion is constructed and later 
removed. Larger or more exposed breakwaters often include features which make 
construction exclusively with land-based equipment difficult. For example, 
placement of large armor units in relatively deep water near the toe of a 
shallow slope (perhaps at the head) may be too far to reach for a mobile crane 


on the breakwater crest. Another example is the occasional need to build up 
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Figure 7. Typical toe trench 


underwater features using floating equipment, particularly toe berms, prior to 
placement of core material, underlayers, and primary armor. The sequence of 
operations, specific placement techniques, and the associated equipment avail- 
able to perform this work usually constrain the range of alternate breakwater 
configurations to some degree. No breakwater configuration should be con- 
ceived without thorough attention to its method of construction. More de- 
tailed discussions of these considerations are available in Bruun (1979), 
Kjelstrup (1979), Maquet (1984), and Bruun (1985). 


Physical Modeling for Stability 


Guiding principles 

42, The specific techniques applied in physical modeling of rubble- 
mound breakwaters by various hydraulic laboratories differ in detail, but the 
guiding principles of similitude offer the same basic constraints in all 
cases. Scale models of breakwaters for hydraulic stability are designed ac- 
cording to the Froude scaling relation which requires that the Froude number 


of the model be equal to that of the full-scale prototype in its intended 


Cat 


natural setting. This relation is expressed as 
ve ua 
gL} ~ \eL eg®) 
a g D 


g = acceleration of gravity 


where 


V = flow velocity 


L = a linear dimension associated with the flow 
These scaling criteria provide that the linear dimensions of the model are all 
geometrically similar to those of the prototype. Typical rubble-mound break- 
water model scales range from 1:5 to 1:70. The Froude number theoretically 
represents the ratio of inertial to gravitational forces, an appropriate mea- 
sure in situations where gravity is the predominant force. It is widely 
accepted that this is usually the case for rubble-mound breakwaters (Hudson 
et al. 1979). 

43. Another scaling law sometimes applies, however, which requires that 


the Reynolds numbers of the model and prototype be equal, or 


LV LV 
m p 


where v is the kinematic viscosity of the fluid. The Reynolds number the- 
oretically represents the ratio of inertial forces to viscous forces. Viscous 
forces in the primary armor layer, underlayers, and core are now thought to 
have greater importance than they did in the pioneering days of rubble-mound 
breakwater design. The Reynolds criterion conflicts in many instances with 
the Froude criterion in sizing structural materials for models (particularly 
in smaller, more economical models), and compromising measures are usually ne- 
cessary. Other scale effects can come into play when model waves are so short 
that surface tension has a significant effect (seldom a real problem in prac- 
tice) or when the mechanical strength of armor units is critical. Dealing 
with these conflicting criteria makes physical modeling of rubble-mound break- 
waters a highly specialized practice. Proper execution of a rubble-mound 
breakwater scale model study requires both specialized equipment and extensive 
experience available only at a handful of hydraulic laboratories around the 


world. 


Operational procedures 


44, The representation of the sea state in scale models continues to 
improve in modeling facilities because of enhancements of wave generating 
equipment and improved understanding of the physics of water waves. The ear- 
liest wave generators were capable only of a sinusoidal motion generating 
monochromatic waves. The last decade has seen these facilities replaced in 
many laboratories with wave generators capable of producing irregular waves 
which simulate specified prototype energy spectra or irregular time series. 
The techniques for application of monochromatic waves are somewhat standard- 
ized, but presently there are widely differing opinions on the most appro- 
priate application of irregular waves in scale models of rubble-mound stabil- 
ity. The transformation of the waves from deep water to shallow water must be 
arranged to be equivalent to that in the prototype for both monochromatic and 
irregular waves. The shallow-water waves of interest for stability are usu- 
ally taken to be those naturally transformed waves that would exist at the 
site without the structure in place. This convention usually involves a cal- 
ibration of the model facility before the model structure is placed in a flume 
or basin. 

45. Complications with reflected waves arise after the structure is in 
place. Techniques are available for analysis of model wave data which resolve 
incident and reflected waves (Goda and Suzuki 1976). Some facilities are ca- 
pable of compensating for reflected waves by modified motion of the wave gen- 
erator. It is necessary in facilities without this capability to reduce wave 
reflection as much as possible by various other means. Use of irregular model 
waves can also result in spurious long-period waves (Jensen and Kirkegaard 
1985) which must be compensated for by the generator or in the interpretation 
of measured results. 

46. Model breakwater materials must reflect a number of prototype con- 
ditions, including geometry, density, surface roughness, and orientation in 
the structure. A number of recent tests have also involved attempts (the re- 
sults of which remain in question) to estimate mechanical stresses within 
armor units (Timco 1981 and Delft Hydraulics Laboratory 1985). Geometry, den- 
sity, and surface roughness are controlled by careful choice of model mate- 
rials and preparation of the units. Minor density scale effects due to use of 
fresh water in a model of a saltwater site can usually be compensated for by 


small adjustments to the weight of the model breakwater units. Orientation 
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in the structure is accomplished with a variety of manual and automatic tech- 
niques designed to simulate the realities of full-scale field placement. The 
placement tolerances of model rubble-mound breakwaters often are smaller than 
their prototype counterparts, however. The hydraulic characteristics of un- 
derlayers and the core can be especially difficult to model by the Reynolds 
criteria since the shape of the units, their surface friction (particularly 
between layers), and the shape of the interstices are critical. Erosion of 
fine foundation material at the toe of breakwaters is also a problem, gener- 
ally yielding only qualitative conclusions. An account of these and other 
scale effects is necessary for reliable interpretation of model results 
(Jensen and Klinting 1983). 

47. Scale modeling operational procedures associated with the design 
of rubble-mound breakwaters can be classified in three general groups: 


(a) cross-section design tests run in two-dimensional flumes (Figure 8), 
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Figure 8. Scale model testing in 
two-dimensional wave flume 


(b) tests of heads, elbows, transitions, offshore hydrographic effects, and 
oblique waves in three-dimensional wave basins, and (c) tests of breakwaters 
at various stages of construction (in either flumes or wave basins). The 
first of these is of primary interest to discussions of analytical optimiza- 
tion, since it is this type of scale model testing which has generated most of 
the analytical relations used by designers. These tests of proposed cross- 
section designs are intended to verify the predictions of analytical proce- 


dures and to refine detailed features of the cross section. They are often 
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more than a fail/no fail "proof test" of a design and should be arranged to 
provide the maximum information of use for similar future designs. 

48. A procedure used for many years to verify the stability of proposed 
cross-section designs involves subjecting a model breakwater to a short series 
of monochromatic waves at the design (stability) wave height at various wave 
periods above and below the design period. The water level is also varied 
within the range of possible levels predicted for the prototype site. This 
sensitivity analysis approach is intended to reveal the breakwater's response 
to the severe condition when plunging breakers are directly impacting the sea- 
ward face, as seen in Figure 8. Displacement of some fraction of the armor 
layer is measured by before and after soundings of the model structure. This 
procedure is relatively economical and provides an indication of the design's 
resistance to armor unit displacement by a group of waves with a "worst case" 
combination of period and water level. Some statistical confidence is lost 
since the design criteria for wave period and water level are not held con- 
stant in the modeling procedure. Subsequent changes to the cross section in 
response to unacceptable damages in the model contribute to further departure 
from initial design criteria and any associated risk analysis. Design cri- 
teria must then be reformulated and associated analyses repeated with the new 
criteria. 

49. Tests of cross-section designs with irregular waves typically 
involve a longer series of waves, since a significant number of waves (100 or 
more) are necessary to adequately resolve a specified energy spectrum. The 
added test condition parameters related to reconstructing a specific spectral 
shape in a wave flume discourage the sensitivity analysis method described 
above. Hydraulic laboratories differ in their approach to tests for the 
effect of wave groups with irregular waves, however. Some favor manipulation 
of spectral shape parameters to enhance wave groupiness, while others prefer 
spectra that are as natural as possible. Recorded spectra are reproduced in 
some instances to assure a completely natural incident wave condition in sta- 
bility tests. Durations of individual tests also vary from relatively short 
tests of around 100 waves (30 to 45 min) to tests of thousands of waves and 
many hours simulating the growth and decline of a storm, as illustrated in 
Figure 5. Further discussion of model tests with irregular waves is available 
in Jensen (1984) and Bruun (1985). 


50. Evaluation of damages after a test is a critical step which 
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requires special care and can involve sophisticated techniques and equipment. 
Color coding armor units in their initial placement is a simple way of illus- 
trating the degree of overall displacement of the armor layer. Soundings on a 
small grid before and after a test will measure the overall volume of material 
which was moved, though net profile changes can hide more drastic gross move- 
ments which may have occurred during the test. Photographic or video proce- 
dures have been used to follow actual movements, including rocking in place, 
of individual units with good success (Delft Hydraulics Laboratory 1985). 
Detection of rocking is especially important in testing dolosse or other 
slender concrete armor units since it is known that they experience signifi- 
cant breakage in place from impacts between individual units. Testing for 
damage rates of these units is therefore a highly subjective process because 
the excessive mechanical strength of model units prevents evaluation of the 
stability of a design after some of the armor units have broken. Reduced 
strength model units (Timco 1981) may eventually provide a better means to 
measure stability of slender units, but model units fully similar to their 
prototype units in mechanical strength are not currently available. Stresses 
in prototype concrete armor units are far from fully understood, but research 
in this area is under way at most leading hydraulic laboratories. 

51. A number of other characteristics are sometimes measured in con- 
junction with stability tests of breakwater designs, including reserve sta- 
bility and wave transmission. Reserve stability refers to the extent of dam- 
age that occurs when the breakwater is subjected to waves in excess of the 
design condition, an important consideration in risk analyses. Wave trans- 
mission characteristics require additional tests to be fully defined, par- 
ticularly when the functional performance design criterion (in terms of wave 
transmission) is substantially different from the structural integrity design 
criteria. Runup is a useful parameter to measure in conjunction with wave 
transmission tests, since the ratio of runup to freeboard seems to be the most 
sensitive parameter in analytically predicting wave transmission by overtop- 
ping. Runup is difficult to gage precisely on rough permeable slopes, and 
traditional visual methods are still common. Measurement of volumetric over- 
topping rates is also occasionally of interest, but a special setup with pro- 
visions for containing overtopped water is necessary. Techniques for measur- 
ing and evaluating the detailed relationship of. runup, volumetric overtopping, 


and transmitted waves to incident waves in terms of wave-by-wave effects and 
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time series analysis parameters need a great deal of further development. 

52. Tests in wave basins to determine the overall susceptibility of a 
rubble-mound breakwater to direct and oblique wave attack (with attention to 
the head and elbows), to transitions between cross sections, and to the hydro- 
graphic features offshore of the structure, are necessary for most projects. 
Wave basin facilities are larger and more complex than wave flumes (Figure 9) 
and therefore are more expensive to use. Tests of this nature not only reveal 
unique information about breakwater stability and other characteristics but 
also provide important confirmation of conclusions from flume tests. Basin 
models are typically at smaller scales than most flume models; thus, Reynolds 
scale effects are exaggerated. Basin model testing confirms the location of 
the most critical cross section of which more precise stability tests should 
be performed in a flume at larger scale. Wave transmission by diffraction 
through the entrance channel or other breakwater gaps is one of the most im- 
portant measurements in a basin test. Long-period oscillations resulting 
from the enclosure of a harbor area by a breakwater are also important to 
detect. Model tests including tidal fluctuations can reveal circulation pat- 
terns inside a proposed breakwater (Headquarters, Department of the Army 
1984). 

53. The last category of breakwater model test is most important for 
large breakwaters requiring complex construction procedures and many months of 
construction time. Provisions for interim protection of partially completed 
breakwaters must be tested to justify what can be a significant additional 
cost to the project. Wave basin testing is more often appropriate for this 
work, but flume testing can be quite helpful also. 

54, The modeling procedures discussed above are the true basis of 
virtually all the analytical tools available to rubble-mound breakwater de-, 
signers. Quantitative measurements of prototype breakwater performance are 
just now becoming available and have yet to be applied toward reliable ana- 
lytical design procedures. Each application of analytical procedures is an 
interpolation or extrapolation of limited prior experience. More often than 
not, refinements which reduce cost and improve performance result from model 
tests of a proposed design. Vital confirmation of analytical assumptions, 
both explicit and implicit, is provided by even the simplest model test. The 


expense and time are worth it in every case. 
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PART III: ESTIMATING DAMAGE RATES 


55. A key step in identification of an optimum among alternative 
rubble-mound breakwater plans is to estimate the expected damages and life 
cycle costs of related maintenance and repairs. The concept of designing a 
rubble-mound breakwater for zero damage is unrealistic because a finite risk 
always exists for the stability criteria to be exceeded in the life of the 
structure. The stochastic nature of the physical phenomena affecting coastal 
engineering structures requires that a probabilistic approach be applied, if 
these maintenance cost estimates are to be more than guesses. The incident 
wave climate can be characterized by estimating probability distribution func- 
tions by a number of relatively well accepted methods (Battjes 1984). The 
crucial problem for rubble-mound breakwater designs is in relating a given 
level of damage and associated repair costs to specific incident wave condi- 
tions. The rate at which this damage accumulates must also be predicted in 
order to tentatively schedule maintenance and related cash flows. The fol- 
lowing section will review some techniques proposed for making these predic- 
tions. Their relative merits will be discussed and areas of ongoing or needed 


future research identified. 


Damage Assessment 


56. Damages to rubble-mound breakwaters have been quantified in many 
ways by researchers and field engineers. The current issue surrounding break- 
age of concrete armor units has led to a number of recent publications pointed 
at systematic assessment of damages of all kinds. One useful characterization 
of prototype damages in terms of displaced primary armor units was proposed by 


Groeneveld, Mol, and Den Boer (1984) and is presented in Table 2. 


Table 2 
Classification of Breakwater Damage 


Type of Failure Displacement, 7% Description 


Minor 0-3 A few individual units of top layer dis- 
placed, but no gaps in top layer larger 
than 4 units; bottom layer intact 


(Continued) 
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Table 2 (Concluded) 


Type of Failure Displacement, b Description 


Moderate 3-5 No gaps in top layer larger than 6 units; 
only slight displacements of bottom units 


Major 5-30 Top layer removed over a large area; bottom 
layer over not more than 2 units 

Total Over 30 Primary armor and underlayers removed over 
a large area with exposure of core 
material 


57. This classification of prototype damages is realistic as far as 
field reconnaissance of a damaged breakwater is concerned, but it departs 
somewhat from the convention of detecting incipient damage in model tests. It 
does not take into account any concrete armor units which have broken in 
place. This inadequacy is compensated for, in part, by the displacement of 
intact primary armor units being accompanied, in most instances, by concurrent 
displacement of broken pieces. It is the exposure and, ultimately, the ero- 
sion of underlayers and core that spell the actual failure of a rubble-mound 
breakwater in the functional sense, with the exception of the case when a 
monolithic crest element has been rendered ineffective. Field investigators 
should also search for evidence of other modes of failure besides hydraulic 
' displacement, including sliding due to toe failure, excessive foundation set- 
tlement, and seismic displacements. Classification of damages as a function 
of both cause and effect is discussed in detail in Bruun (1985). 

58. Laboratory investigations, as pioneered by Iribarren (1938) and 
Hudson (1958), typically attempt to identify the point of incipient damage. 
Kogami (1978) defined this criterion as "...the condition in which the number 
of armor units clearly recognized to have been moved or rocked on the cover 
layer surface by wave actions was less than 1% of the total of the units on 
the forward cover layer...." The account of rocking implies that a precise 
method of measuring the extent of rocking is available. Another interpreta- 
tion relates to the point at which displacement has reached a depth in the 
armor layer equal to the equivalent cube dimension of the armor units (Losada 
and Gimenez-Curto 1979). Techniques developed by WES in the 1950's for mea- 
suring model breakwater displacements with before and after soundings have 
been estimated to have a resolution (repeatability) of +2 percent (Carver 


1983). Identification of incipient damage with this commonly used method 
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is therefore only meaningful in the range of 0-3 percent primary armor dis- 
placement. Nielsen and Burcharth (1983) have indicated that measurements 
of very low levels of displacement or rocking (0-3 percent) are less reliable 
than those of higher levels. This trend relates to the resolution of mea- 
surement techniques as well as the repeatability of the experimental results 
themselves. 

59. Given a relatively consistent and precise method of measuring 
displacement, Ahrens (1984) has proposed a useful dimensionless parameter for 


systematic quantification of breakwater damage: 


4D 
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Where Ap is the average eroded cross-sectional area for a specific length of 


D' (12) 


model breakwater (Figure 10). Van der Meer and Pilarezyk (1984) applied the 
following dimensionless damage parameter Sp in their model tests of quarry- 
stone, which was mentioned previously (Equations 5 through 8) in the discus- 


sion of their conclusions regarding stability: 
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It is also important to identify erosion of the underlayers or core that may 


S (13) 


coincidentally occur with erosion of the armor layer. 


Analytical Damage Prediction 


60. Scale model studies reported by Jackson (1968a) and Carver and 
Dubose (in preparation) have addressed, to a limited degree, the level of dam- 
age to breakwater armor layers experienced when the design wave height is ex- 
ceeded. This information was applied to formulate Table 7-9 in the SPM (1984) 
which predicts the percent damage %D for various armor types as a function 
of the design wave exceedance ratio H/Hy where H is a monochromatic inci- 
dent wave height which is greater than the design wave height Hy - The re- 
serve stability trends, or tendency for damage levels to increase with design 
wave exceedance ratio, can also be characterized by a function of the follow- 


ing form: 
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7 [s (H/H4-1)] 
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where 


%D(H/H,) = damage experienced by a particular armor type from an 
d yi 
incident wave height H , given a design wave height Hy 


%4D(H,) = level of incipient damage detectable in the model tests 

d g 
which identifies the damage trend for a particular armor type 
(i.e. when H/Hq = 1) 


S, = an empirical coefficient fit to the scale model test results 


16 
for a particular armor unit type 

61. A higher S, coefficient means that an armor unit type experiences 

higher damage levels for the same increase in H. Table 3 gives the values 


Table 3 
Coefficients for Analytical Prediction of Breakwater Damage 


Armor Wave 


Unit Type Condition BD (Hq) Sr Data Source 
Quarrystone Nonbreaking 3:0 6.95 Jackson (1968a) 
(rough) 

Quarrystone Breaking 2.0 3.165 Carver and 
Dubose (1985) 

Quadripods Nonbreaking 3.0 6.00 Jackson (1968a) 

Tribars Nonbreaking 3.0 4.87 Jackson (1968a) 

Dolosse Nonbreaking 2.0 1.68 Carver and 
Dubose (1985) 

Dolosse Breaking 2.0 3°55 Carver and 


Dubose (1985) 


for #%D(H4) and S,, found for the armor unit types which have been tested 
at WES. These coefficients may be used with caution in Equation 14 to pre- 
dict breakwater damage. The variation in mD(Hq) between armor unit types 
reflects improvements in the accuracy of damage measurements as well as dam- 
age trends which may be related to armor unit characteristics. One or more 
statistical outliers representing more severe damage than predicted by Equa- 
tion 14 (presumed to have been caused by weaknesses other than hydraulic sta- 
bility) were excluded from the analysis of data for each armor type. The dam- 
age predicted by Equation 14 at selected levels of design wave exceedance and 


the associated upper 95 percent statistical confidence limit are presented in 
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Table 4. Equation 14 predicts the statistical mean trend of the experimental 
data. Since this is the most probable damage level for a given H/Hqy ratio 
based on the empirical evidence available, it is appropriate for application 
in estimates of expected damage. Designers should be sure also to consider 
the damage predicted by the upper 95 percent confidence limit of the pertinent 
model test (as shown in Table 4) and report these predictions in their docu- 


mentation of the design analysis. 


Table 4 


Damage Level Predictions at Selected Design Wave Exceedances 
in Percent Displacement of the Armor Layer* 


Mean Trend/95% Confidence Limit 


4 Quarrystone Quarrystone Quadr ipods Tribars Dolosse Dolosse 

eis (Nonbreaking) (Breaking) (Nonbreaking) (Nonbreaking) (Nonbreaking) (Breaking) 
1.00 3.0/24.7 2.0/10.2 3.0/18.6 3.0/10.4 220/423 2.0/7.4 
1.05 4 .2/25.7 2.4/10.6 4.0/19.5 S2O/Alileal 2ee/4e5 2.4/7.8 
U0 6.0/27.2 2.9/11.0 5.5720. 7 4.9/12.1 2.4/4.6 2.9/8.2 
sass 8.5/29.6 325/11-36 7.4/22.6 6.2/13.4 2.6/4.9 3.4/8.7 
1.20 12 733a1 4.1/12.2 10.0/25.1 7.9/15.1 2.0/5ai 4.1/9.4 
acs) We 1/3852 5.0/13.0 13.4/28.6 10.1/17.3 8 .10/51.3 4.9/10.1 

130 24 ..2/45.5 6.0/14.0 18.1/33.3 12.9/20.1 3.3/5.6 587 1 

1.35 34 .2/55.8 We? 1Dice 24 .5/39.8 16.5/23.7 3.6/5.9 6.9/12.2 
1.40 48 .4/70.4 8.6/16.7 33 .0/48.5 21.0/28.4 3.9/6.3 8.3/13.6 


* Displacement of more than 30-40 percent of the armor layer will often involve erosion 
of underlayers, which in practice requires a repair effort of greater scope than re- 
placement in kind. 


62. Tables 3 and 4 include predictions for only four types of armor 
units, two of which do not include breaking wave conditions. This is unfortu- 
nate, but it leaves the designer with no option but to apply subjective judg- 
ment to choose damage coefficients which are close to those of the most simi- 
larly shaped armor unit in the same wave conditions. Slender concrete armor 
units, including nearly all concrete types more complex than plain cubes, are 
subject to breakage in place from impacts between individual units in the 
armor slope. This breakage would presumably be accompanied by displacement of 
the broken pieces during an extreme storm. An increase in S, of 50-100 per- 
cent would provide some allowance for this likelihood, but there are no data 


currently available with which to more precisely predict breakage or its 
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above empirical results are from monochromatic model tests of limited duration 
which do not account for the natural irregularity of ocean waves nor the ef- 
fect of variable duration of exposure. The many untested or otherwise unre- 
solved questions about breakwater damage modes should not, however, prevent 
designers from applying the information that is available. The need to con- 
firm analytical predictions of breakwater stability and performance by scale 
model testing prior to construction cannot be overemphasized. 

68ha59..8 characterization of damage as a function of incident wave height, 
with the features of Equation 14, allows the "expected" or long-term average 
damage to be estimated. The statistical definition of expectation for con- 


tinuously distributed variables is 


E{x} = f xt dx (15) 


where f(x) is the probability density function (pdf) of x (DeGroot 1975). 
A function of x , g(x) may be substituted for x in Equation 15 without 


changing the definition, thus 
Bg (x) = f eet) dx (16) 


The long-term distribution of wave heights formulated for most current design 
exercises to represent the incident wave climate is derived as a cumulative 


probability distribution (cpd) F(H) where 


dF(H) 


£ CH) a= qH 


(17) 


The expected annual damage can then be estimated from a damage function 


*D(Hq) , Such as Equation 14, and a epd for wave heights F(H) by using the 


px =f w(F -)| ae (18) 


where i = the Poisson parameter or average number per year of extreme 


following equation: 


events represented by H values. This formulation assumes that the number 
of storms per year is a random variable and can be represented by a mean 


value. It assumes further that this number is independent of the H_ values 
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which represent the intensity of the individual storms. 

64. The availability of synoptic hindcast data base of wave data for 
most of the US coastline (Corson et al. 1981) accommodates the technique for 
formulation of F(H) where only the significant wave height H, values 
(representing the intensity of a severe storm) above a threshold value are 
addressed (Battjes 1984). Recent applications of hindcast wave data at WES 

(Andrew, Smith, and McKee 1985) have yielded good results with a cpd function 
for significant wave heights above a threshold using the following extremal 


(Fisher-Tippet) Type I distribution: 


=H.) 
ll - 4] (19) 


FH = 
s 


dF(H,) FHA) [(e-H,)/0] 
dH. = ae e (20) 


i] 
where E(H,) is the cumulative probability that a significant wave height H, 
in a sample is equal to or less than some specified Hy , OF PH, < Hed WE 
and are parameters fit to the data by regression. The traditional return 


period RT can be estimated as (Borgman and Resio 1982) 


1 


es MPT - FHS) (21) 


65. Another commonly applied cpd, traditionally used for annual ex- 
tremes, is the following Weibull distribution: 


c 
F(H,) = 1 - [#24] 


(22) 
where C is an additional empirical parameter which must be fit to the data. 
This distribution is equivalent to a Rayleigh distribution when C = 2 and 
reduces to an exponential distribution when C=1 and e¢« =0 (Petraukas 

and Aagaard 1970). 

66. Either of these cpd functions could be applied to estimate the 
expected damages, given a damage function such as Equation 14. These cpd 
functions are typically applied to present the probability of exceedance 
for a specified H, (i.e. P(H, bs )). Assignment of a representative 


unit price for repair of displaced armor units allows the expected cost of 


42 


damages E{$D/yr} to be estimated for a breakwater design, which is the same 
as the "equivalent annual amount" that might be derived by discounted cash 
flow analysis. An interactive FORTRAN computer program called "BWDAMAGE" has 
been developed at WES. This program estimates E{$D/yr} given values of 

%D( Hq) and S, for Equation 14, ce and o for an Extremal Type I ecpd of 
significant wave heights (Equation 19), representative armor repair unit 
prices and the volume of the armor layer. This program is documented further 
in Appendix D of this report. Its intended use is for comparison of alter- 
native plans, and for this purpose the limited statistical confidence of the 
applied formulae is acceptable. Substitution of a measured damage function 
from model tests of a particular design would greatly improve the reliability 
of the program's estimates. 

67. A number of refinements to the above scheme of analytical predic- 
tion of rubble-mound breakwater damage are conceivable. The effects of wave 
period and storm duration on stability have been recently investigated by a 
number of specialists. The effects of wave period and storm duration were 
incorporated directly into the stability formulae proposed for quarrystone 
rubble-mound breakwaters by Van der Meer and Pilarezyk (1984) (Equations 5 
through 8). The joint effect of wave height and period on armor damage, in 
the form of the surf parameter (Equation 4), and risk analysis in terms of a 
probability distribution of wave steepness ‘is discussed in Bruun (1985). The 
DHI-Iribarren stability formula (Equation 3) directly incorporates the effect 
of wave period as the corresponding wave length (Graveson, Jensen, and Soren- 
sen 1980). This latter work also addresses the effect of storm duration by 
focusing on the rate at which damage occurs for variations of the other sta- 
bility. related ‘factors -(W , Hig oT ';.cos“®. ,.cte.).: The, following relation, of 


damage was derived from the data of Graveson, Jensen, and Sorensen (1980): 


K Sian 
B MUO M23) 
De. = "0! 0622 7,000 
where 
D, = damage rate D/t , in percent armor displacement per hour 
o = the DHI-Iribarren stability coefficient 


= 0 Bw HEL, /WA9 (1 cos 6 - sin 9)3 from Equation 3 


K is a function of wave height and period, so all three parameters (H,T , 


t) are also included in the DHI approach, since 
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ke ae 
fe) 
DDCHAeL it) m= OOG22E 1,000 (24) 
where t can be taken as the average duration of exceedance of Reus in 

Equation 3. This duration is difficult to assess in practice. Investigations 
of the long-term joint probability distribution F (H , T , t) are needed for 
a more precise definition of this parameter. If the pdf f(H , T , t) could 
in turn be estimated, a particular rubble-mound breakwater design could be 


evaluated for its expected annual damages by 


E = : fff ou, T, t) £(H, T, t)dHdTdt (25) 


68. The practical problem in applying Equation 25 is estimating the 
joint pdf f(H , T , t) . An interim approach to account for duration would 
be to assume an average t for all storms exceeding the design condition, 
based on evaluation of hindcast statistics or other long data records. Like- 
wise, characteristic peak periods and water depths d can be associated with 
extreme storms in most cases without rigorous definition of the joint pdf or 
epd. This is already common practice, since a design wave period and water 
surface elevation have always been necessary for accomplishment of wave trans- 
formation analyses and estimates of runup, overtopping, and wave transmission. 
Methods for estimating the joint long-term probability distribution of H and 
T are discussed by Sigbjornsson, Haver, and Morch (1976) and Ochi (1980). A 
practical approach to estimating expected damage by use of the DHI-Iribarren 
formula, given appropriate wave data, is proposed in Jensen (1984). Assump- 
tions concerning the mean direction of wave propagation 9 and the associated 
directional spreading o are also inherent in current practice for defining 
the wave climate at a site. The effect of wave direction on rubble-mound sta- 
bility and damage rate is discussed by Losada and Gimenez-Curto (1982) and 
Christensen et al. (1984). Estimation of expected damages in terms of D(H, 
T, t, d; ¢, o) and f(H , T, t,d, 4, o) will be possible only after much 


additional theoretical, laboratory, and field investigation. 
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PART IV: ESTIMATING WAVE TRANSMISSION 


Wave Transmission by Diffraction 


69. Waves are transmitted by rubble-mound breakwaters in three ways: 
around, over, and through. The first way plainly refers to diffraction of 
incident waves around the heads of breakwaters at the entrance channel or 
through other gaps in the structure. Wave transmission by diffraction, the 
most substantial of the three modes, can be limited by careful orientation of 
the breakwaters. Diffracted waves combine with waves transmitted over and 
through a breakwater within the area influenced by diffraction. All three 
modes must be addressed in this area. Methods to define the limits of pene- 
tration of diffracted waves, including estimates involving directional irregu- 
lar incident waves, are presented in the SPM (1984) (see also Goda, Takayama, 
and Suzuki 1978). 

70. Many projects, such as boat harbors and ports where the breakwater 
is relatively extensive and the principal physical feature providing wave pro- 
tection, can deal with optimization of the breakwater in plan as a separate 
measure. This optimization can precede the optimization of the breakwater 
cross section and include layout of all the other major features associated 
with the proposed coastal development. Procedures for systematic optimization 
of breakwater lengths and orientation with respect to wave penetration by dif- 
fraction are discussed by Groeneveld et al. (1983) and in EM 1110-2-1615 
(Headquarters, Department of the Army 1984). 


Wave Transmission by Overtopping 


71.  Rubble-mound breakwaters are designed usually with the intention 
that waves do not overtop the structure except in the most extreme incident 
wave conditions. Traditionally this has been a matter of estimating runup on 
the seaward face for an extreme wave height and period combination and setting 
the crest just above the maximum runup. This method can provide a crude ap- 
proximation of crest elevation for concept formulation, but it should not be 
carried further into the design process. More precise techniques for esti- 
mating wave transmission by overtopping were devised by Cross and Sollitt 


(1971) and later refined by Seelig (1980b). 
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72. The height of a wave transmitted by overtopping has been found to 
be a function of incident wave height, period, freeboard (vertical distance 
from the crest to the mean water level), slope, crest width, and surface char- 
acteristics affecting runup. Water depth and bottom slope at the toe of the 
structure also affect wave transmission by overtopping to the extent that they 
affect the characteristics of the incident wave. The reflection characteris- 
tics and permeability of the structure also have an effect. Figure 11 illus- 
trates incident wave energy being partially reflected, partially dissipated in 
turbulence at the seaward face, and partially dissipated by viscous effects. 
Energy not reflected or dissipated in these ways either passes through or over 
the breakwater, or both. The explicit method developed by Seelig (1980b) for 


predicting wave heights transmitted by overtopping is as follows: 
He = Keo (Hj) (26) 


where 


Hy = transmitted wave height 


Kio = transmission coefficient (by overtopping) 
F 
= c(1 - *) (27) 
C = an empirical coefficient 
B BL 
= 0.051 - 0.11 ny (fer h, < 3.2) (28) 


F = freeboard 
R = potential runup, as if the seaward slope were infinitely high 
B = crest width 

h, = total height of the crest above the sea bottom 


x 
" 


incident wave height 
73. Runup can be estimated by a number of methods, but the method de- 
veloped by Ahrens and McCartney (1975) is particularly useful for analysis by 


the wave transmission formula above. It is expressed as 


td 


t ag 
a Gi be) (29) 


where a and b are empirical coefficients associated with the particular 
type of armor unit in place. In this case, the surf similarity parameters €& 


(Equation 4) is related to the incident wave height, the equivalent deepwater 
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wave length of the incident wave period, and the slope of the seaward face. 
Values of a and b_ have been derived by Seelig (1980a) from monochromatic 
laboratory data for riprap revetments (graded quarrystone) on an impermeable 
surface, uniform quarrystone on both highly permeable and conventional multi- 
layered breakwaters, and for dolosse on conventional multilayered breakwaters. 
From experiments with conventional multilayered breakwaters, additional values 
have been fit for this report to monochromatic runup data taken by Jackson 
-(1968b). The values of these runup coefficients are presented below in Ta- 


ble 5 along with the linear correlation coefficient r to the data from which 


Table 5 
Runup Coefficients 


Armor Unit a b ig Data Source 

Riprap 0.956 0.398 -- Ahrens and McCartney (1975), 
(revetments) impermeable base 

Quarrystone 0.692 0.504 -- Hudson (1958), 
(breakwaters) highly permeable core 

Quarrystone 0.115 0.361 -- Gunbak (1976), 
(breakwaters) multilayered 

Modified Cubes 0.95 0.69 0.91 Jackson (1968a), 
(breakwaters) multilayered 

Tetrapods 1.01 0.91 0.76 Jackson (1968a), 
(breakwaters) multilayered 

Quadripods 0.59 0.35 0.83 Jackson (1968a), 
(breakwaters) multilayered 

Hexapods 0.82 0.63 0.78 Jackson (1968a), 
(breakwaters) multilayered 

Tribars 1284 1.5% 0.78 Jackson (1968a), 
(breakwaters) multilayered 

Dolosse 0.988 0.703 -- Bottin, Chatham, and Carver (1976), 
(breakwaters ) multilayered 


they were derived. The relation of the runup predicted using these units as a 
function of & is illustrated in Figure 12. An important feature to note is 
that some armor unit types may have runup advantages over other types in that 
they can be more efficient energy dissipaters with respect to runup. Some of 
this effect may be due to variations in underlayer material size and porosity 
that are functions of the primary armor unit weight as well as the total depth 


of the primary armor. This means that breakwaters built with certain heavier 
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Figure 12. Relative runup versus the surf parameter, 
(after Ahrens and McCartney 1975) 


armor units may have lower crest elevations and, in some cases, less overall 
volume and cost. This aspect of armor unit characteristics has not been very 
well explored to date. 

74. A variety of interpolation schemes based on other armor unit param- 
eters (including stability coefficient, layer porosity, layer coefficient, and 
combinations of these parameters) failed to yield results similar to the a 
and b values directly fit to runup data. The marginal correlation of a 
and b for some armor units to the Ahrens and McCartney (1975) runup equation 
is an indication that further carefully controlled runup experiments are badly 
needed. Runup is difficult to measure precisely with instruments on a rough 
permeable slope; therefore, the above data were measured primarily by manual 
means. Losada and Gimenez-Curto (1980) also investigated runup on breakwaters 
as a function of & , which is expressed as 


BE 


R ; 
Tie A(1 - e 7) (30) 


49 


where A and B are empirical coefficients yielding runup trends very simi- 
lar to those proposed by Ahrens and McCartney (1975). Their regression showed 
similar trends and correlation, as indicated in Table 6 and Figure 13. 
Table 6 
Runup Coefficients of Losada and Gimenez-Curto (1980) 


el 


Correlation 
Armor Unit A B Coefficient 
Riprap 1.789 -0.455 0.96 
Quarrystone We oy -0.523 0.81 
Quarrystone 1.370 -0.596 0.61 
Tetrapods 0.934 -0.750 0.74 
Dolosse 1.216 -0.568 0.74 
Quadripods 1,538 -0.248 0.86 
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Figure 13. Relative runup versus the surf parameter 
(after Losada and Gimenez-Curto 1980) 


75. Irregular runup can be predicted, based on either of the above re- 
lations, by applying a joint cpd for the sea state (Longuet-Higgins 1975 and 
Ochi 1980) to predict the runup of each wave as a function of its steepness 
H/Lo . This process applies the principle of equivalence, first proposed by 


Saville (1962), which assumes that the effects of each wave in an irregular 
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sea state may be represented by the effects of an equivalent monochromatic 
wave of the same height and length. Losada and Gimenez-Curto (1980) have ap- 
plied this principle and Equation 30 to several joint distribution functions 
for H and T to derive distributions of runup that compare well with experi- 
mental data. The SPM (1984) proposes a more expedient method which assumes the 
runup heights will have a Rayleigh distribution. An alternative expedient 
method has been proposed by Andrew and Smith (in preparation) which assumes a 
Rayleigh distribution of wave heights and a constant wave period equal to the 
period of peak energy density. The resulting distribution of runup heights is 
not Rayleigh distributed, in keeping with the joint effect of height and pe- 
riod as predicted by the runup formulae above. Interactive programs written 
in BASIC for microcomputers are available from WES to estimate both runup and 
wave transmission by overtopping by this technique.* An example of trans- 
mitted wave height exceedance probabilities estimated by this method is pre- 
sented in Figure 14. 


76. The principle of equivalence may not remain the key to prediction 
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Figure 14. Predicted transmitted wave 
height exceedance probabilities 


* US Army Engineer Waterways Experiment Station, Coastal Engineering Re- 
search Center, in preparation, "Wave Runup on Rough Slopes: Computer Pro- 


gram WAVRUNUP (MACE-14)," Coastal Engineering Technical Note, Vicksburg, Miss. 


, in preparation, "Wave Transmission by Overtopping: Computer 
Program WAVTRANS (MACE-13)," Coastal Engineering Technical Note, Vicksburg, 
Miss. 
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of wave transmission by overtopping since investigators have noted that over- 
topping tends to generate waves of much shorter period than the incident wave 
(Jensen and Sorensen 1979 and Jensen 1984). Investigations of wave transmis- 
sion over a natural reef and associated laboratory experiments by Gerritsen 
(1981) resulted in development of a theoretical approach to the redistribution 
of energy that occurs with wave breaking on, and spilling over, a low-crested 
or submerged reef. The transmitted waves were found to be fairly well repre- 
sented as a collection of "solitons" or wave energy packets generated by inci- 
dent breaking waves represented as long waves or bores (analogous to hydraulic 
jumps). The phenomenon of "surf beat" or wave grouping was found to be criti- 
cal to higher levels of energy transfer. The methods of Gerritsen (1981) 
might yield useful results if applied to wave transmission by rubble-mound 


breakwaters. 


Wave Transmission Through Permeable Breakwaters 


77. The tendency of wave energy to permeate through the interior of 
rubble-mound breakwaters can be important for structures with relatively 
coarse core material. Keulegan (1973) performed laboratory experiments of 
this phenomenon which led several others to further theoretical and laboratory 
investigations. Sollitt and Cross (1976) and Madsen and White (1976) devel- 
oped semiempirical techniques to predict wave transmission through permeable 
rubble-mound breakwaters. Wave transmission by this mode was assumed by these 
authors to be a function of wave steepness H/L , structure permeability, 
structure width, and the capacity of the structure to reflect wave energy or 
to dissipate it in turbulence. The theory of long waves was applied to formu- 
late expressions for wave transmission since it was assumed that the waves of 
significant consequence would be much longer than the width of the structure. 
Laboratory experiments indicate this as a practical assumption for most break- 
water sites. 

78. Madsen and White (1976) also developed a computer program for 
predicting wave transmission through multilayered rubble-mound breakwaters. 
This program was refined by Seelig (1980b) who successfully tested its pre- 
dictions against an extensive set of laboratory results to account for com- 
bined wave transmission from overtopping and permeation. It was further 


modified for interactive use (US Army Engineer Waterways Experiment Station, 


By 


Coastal Engineering Research Center (WES, CERC) 1984a) and to incorporate the 
estimation of wave transmission by overtopping for irregular waves as proposed 
by Andrew and Smith (in preparation). This program, titled "MADSEN," is an 
extremely useful tool to analytically predict wave transmission for planning 


purposes where diffraction is not a significant factor. 
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PART V: COST-EFFECTIVE OPTIMIZATION 


The Principle of Optimization 


79. Optimization is referred to as "trade-off analysis" in some Corps 
of Engineers planning guidance (Board of Engineers for Rivers and Harbors 
1985) in the sense that identification of an optimum plan usually requires one 
desirable goal to be compromised or "traded off" against one or more other de- 
sirable goals. The basic trade-off in public works economics can be stated as 
a contest between minimum costs versus maximum benefits. The desired effect, 
such as elimination of damages by wave attack, must be balanced against the 
desired goal of no cost. To eliminate the remotest likelihood of damages, a 
structure might be astronomically expensive to build and maintain. A struc- 
ture in which all but some very remote likelihood of damages is eliminated 
might be much more affordable. The damages or other economic losses and inef- 
ficiencies which are undesirable in their unmitigated state can be associated 
with a level of cost to those who are suffering the losses. The tangible ben- 
efits realized by a public works project are the sum of the incremental reduc- 
tions in that level of costs directly attributable to the functional perform- 
ance of the project. The construction and maintenance costs of the project 
are added costs to the beneficiaries, however. These project costs can be 
considered as negative benefits, thus the optimum plan is the combination of 
features which achieves the maximum net benefits. These maximum net benefits 
must be positive; that is, the benefits must exceed the project costs for fed- 


eral participation to be possible. 


An Idealized Approach 


80. Figure 15 illustrates the principles discussed above in an ideal- 
ized arrangement. The horizontal line labeled "user (total) cost without 
structure" refers to the expected annual economic losses that exist without 
any mitigation. This cost may increase over time due to population increase, 
inflation, or other factors, but it can be represented by an equivalent annual 
amount for the sake of evaluating project alternatives. The representative 
amount of economic losses without the project will be the same for each alter- 


native. The ensemble of alternatives will individually reduce these 
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economic losses by varying degrees. The line labeled "user cost with struc- 
ture" represents the economic losses at the reduced level, and its shape in 
this case indicates that the alternative plans along the x-axis are ordered by 
increasing benefits. The benefits themselves are the difference between the 
"without project" and "with-project" conditions, as indicated by the line la- 
beled "user benefits." 

81. Project alternatives which are built so soundly as to preclude de- 
terioration of any kind would obviously have a tremendous first cost. Most 
projects therefore accept some minimal level of predictable deterioration and 
associated maintenance costs in order to reduce the first costs. to an afford- 
able level. A number of authors have treated this problem as an independent 
matter, taking for granted that a specific level of benefits is to be achieved 
by all alternatives. This approach overlooks the situation in public works 
development in which the "user" or beneficiary is the same agency which must 
pay the life cycle project cost. A true optimum plan must minimize all costs, 
i.e., the economic losses and the structure life cycle costs. 

82. Additionally, Figure 15 shows the hypothetical ensemble of alterna- 
tives to be ordered in terms of increasing first cost, as indicated by the 
shape of the line labeled "structural first cost." The increasing first costs 
are taken in this idealized representation to correspond to reduced mainte- 
nance liability as indicated by the line labeled "structural maintenance cost" 
which slopes in the opposite direction. The sum of these costs for each al- 
ternative is shown as "total structural cost," with a minimum in the vicinity 
of alternative 5. The sum of the total structural cost and the user cost with 
structure is shown as the line labeled "total cost with structure." This line 
dips below the "user cost without structure" line at a point where the bene- 
fits first exceed the costs. The region where benefits exceed the costs has 
been shaded and labeled as "project feasibility." The alternative with the 
maximum vertical spread in this shaded area has the maximum net benefits, in- 
dicated by the optimum point on the line labeled "net benefits." This point 
corresponds to the point of minimum total cost with structure, somewhere 
around alternative 9. 

83. It is useful to note that the optimum can be identified without 
knowledge of the without-project condition. Port and harbor projects are of- 
ten justified in terms of transportation savings over some alternate route or 


through some other existing port. The user cost with project would in these 
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cases be compared to the user cost through the alternate route, for definition 
of project benefits. Whatever economic philosophy or administrative policy is 
applied, an estimate of tangible economic benefits must be made. This esti- 
mate, either with the "user cost with project" or with the benefits them- 
selves, can be applied in the manner of Figure 15 to optimize the major fea- 
tures of the project. 

84. The idealized nature of Figure 15 is useful to illustrate the 
concept of cost-effective optimization, but it is misleading in its impli- 
cation that a set of alternatives will follow such a smooth comparison of 
costs and benefits. A typical set of plans could not, in most cases, be 
ordered by both increasing benefits and first cost. Neither is it the case 
that increasing first cost always means reduced maintenance. The two types 
of design criteria--functional performance and structural integrity--are es- 
sentially independent of each other, and both have an effect on first cost. 
Practical applications require that an ensemble of alternatives be compared 
without reference to the order of their benefits, first cost, and maintenance 
cost. A systematic approach to criteria development as a means of initially 
identifying alternatives is important in this respect. An alternative is thus 
known in the optimization process by its governing design criteria rather than 


its resultant physical features. 


A Practical Approach 


85. The analytical and practical aspects of rubble-mound breakwater 
design have now been reviewed. The discussion above concerning the principles 
governing optimization indicates that the first cost, maintenance cost, and 
user cost with project must be estimated for each plan in a set of alterna- 
tives. The total costs with project need to range at close intervals from 
well below to well above the unknown minimum for reliable identification of 
the optimum alternative. A procedure is proposed below that accomplishes this 
optimization exercise using information commonly available and already incor- 
porated in most coastal engineering planning and design efforts. Potential 
future refinements are mentioned where appropriate. 

Step 1--define site conditions 
86. The physical conditions and other constraints affecting the design 


of a rubble-mound breakwater, such as water level, tidal currents, foundation 


oT 


characteristics, and wave climate, must first be quantitatively defined. It 
is assumed that many of these conditions have been already defined in a master 
planning effort which identified the tentative need for a breakwater and its 
most promising alignment. The water level and wave climate are the most crit- 
ical considerations for this optimization procedure, specifically the estima- 
tion of the annual cumulative probability distribution F GH Oy TV iG te ene 
o , d). Current practices typically require planners and designers to esti- 
mate F(H,) for a limited range of wave directions » affecting the site of 
the breakwater or, at best, to define a wave rose and then deal with the mar- 
ginal distribution of wave heights for one sector. Design values of wave pe- 
riod T , storm duration t , directional spreading o , and depth d asso- 
ciated with a given wave height are typically subjectively determined. The 


mathematical estimation of F(H Ty) is becoming more common, however, with 


’ 
the availability of hindcast oe bases of wave information (Corson et al. 
1981). 

87. The Extremal Type I distribution for F(H,) , based on H, values 
above an extreme threshold value, is recommended in this procedure for designs 
where hindcast information or other comparably long records of wave data are 
available. A Weibull distribution of extremes is a workable alternative. Ap- 
plication software program WAVDIST for estimating Extremal Type I and Weibull 
Significant wave height distributions has been documented in a Coastal 
Engineering Technical Note (CETN),* and an example of its use in a design 


problem is presented in Andrew, Smith, and McKee (1985). 


Step 2--estimate ex- 
pected economic losses 


88. Estimation of losses or "user costs" due to wave attack requires 
derivation of a site-specific relation in which losses are a function of inci- 
dent wave height. The typical harbor mooring area or cargo transfer area is 
unaffected by waves below a certain height Hp, Which might be on the order 
of 1m. The total disruption of the port or harbor area at the other extreme, 
by the worst conceivable wave attack, is also possible to estimate as a prac- 


tical upper limit to losses These two values are useful in that they 


$Linax : 


* US Army Engineer Waterways Experiment Station, Coastal Engineering Research 
Center, in preparation, "Estimation of Extremal Significant Wave Height Dis- 
tributions: Computer Program WAVDIST (MACE-17)," Coastal Engineering Tech- 
nical Note, Vicksburg, Miss. 
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do not require historical information for their estimation. They can be based 
on a current engineering and property valuation assessment of the facilities 
to be protected. Historical information relating specific levels of economic 
loss (in dollars) to the measured or hindcast wave height of the associated 


storm can then be used to derive a function of the form, as follows: 
A(H_-H i 
S, Lo 
$L(H,) = s,_,[1 =e (31) 


where A is a coefficient determined by regression. This function, illus- 
trated in Figure 16, can then be used to estimate the expected annual economic 


losses, or user costs without project, according to Equations 16 and 17 by 


dF(H_) 
$L\ _ sided Sa 
Ea E WE i, fe, (32) 


where 0: ee) is the cumulative probability distribution of significant wave 
heights derived in Step 1. A joint distribution BCHe : Tp) should be ap- 
plied where operations and facilities are particularly sensitive to a certain 
range of periods. Software has been developed to estimate both $L(H,) and 


E{$L/yr} given an estimate of $L Hig » Extremal Type I F (HE) coeffi- 


max ’ 
cients, and at least one historical data point [H, : $L(H.) | . The program 
SLyyax - MAXIMUM CONCEIVABLE LOSSES FROM WAVE ATTACK 


$L (HJ =$lLyax [1-e4"4s~ Heo! | 
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Figure 16. Economic loss function versus incident wave height 
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"BWLOSS1" has been documented in a CETN*, and a sample is provided in Appen- 
dix A of this report. 


Step 3--formulate en- 
semble of alternatives 


89. This step is highly subjective and will control the scope of the 
overall optimization effort since it determines the number of individual al- 
ternative breakwater configurations that must be investigated. The applica- 
tion of practical judgment can reduce this number, but too few alternatives or 
a conservative bias could also preclude identification of an optimum plan. A 
proposed method of organizing an ensemble of alternatives is illustrated in 


Table 7. 
Table 7 


Selection of Alternative Design Criteria for Return 
Periods of Storms Causing the Stated Conditions 


Functional Performance Structural Integrity 
xh H, > H* zD < SD(Hq) 
<0 £30), 40; 50, 60, *70, 80; “9021100 
20 <30;<40, 50, 60, 70, 80, 9077100 
30 <30; 40,5060; 70,80; °90)'"100 
HO <305140;1/50,, 60; 70, 80; 90, 100 
50 <30, 40, 50, 60, 70, 80, 90, 100 


90. Table 7 lists a comprehensive set of potential functional and 
structural design criteria combinations which may be abbreviated by carefully 
considered subjective judgments. The first column in Table 7, "functional 
performance," refers to an exceedance value x% of transmitted wave heights 
H, greater than some critical wave height H* . H* might conveniently be 
taken as the Ho value applied in the loss function of Step 2, but this is 
not necessary. This column includes a range of functional performance design 
criteria which could be addressed in terms of wave transmission. A wave 


height of 1m, for example, might be a threshold value for damage to vessels 


* US Army Engineer Waterways Experiment Station, Coastal Engineering Research 
Center, in preparation," Estimation of Expected Annual Economic Losses Due 
to Wave Attack--Computer Program BWLOSS1 (MACE-15)," Coastal Engineering 
Technical Note, Vicksburg, Miss. 
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moored behind the breakwater. The last functional performance criteria would 
thus be that x%Z of the waves transmitted by the breakwater during a 50-year 
storm would be in excess of 1m. The return period convention is in keeping 
with traditional practice, though the phrase "with 2 percent probability per 
year" would be a more accurate description of the storm of interest. Esti- 
mated probability per year might be a more appropriate increment in terms of 
providing even steps of cost between alternatives, but either convention will 
serve. The value of x% should relate to some consideration of the actual 
number of waves of H* or greater necessary to cause a measurable effect. A 
storm whose peak conditions lasted 3 hr with Tp = 10 sec would include 
roughly 1,080 waves. A small value of x% is appropriate, on the order of 

1 percent, which for the example condition would include 10 or 11 waves. 
These waves would not likely occur in sequence, but a few of them might. 

91. The shorter return periods of 20 or <10 years might be too risky 
for a small boat harbor where relatively fragile vessels and mooring facili- 
ties are planned immediately on the lee of the breakwater. These criteria are 
reasonable, however, when losses due to cargo handling inefficiencies or ves- 
sel transit time are all that is at stake. The 50-year storm is, on the other 
hand, a very conservative criterion for wave transmission. At least four 
functional performance criteria should be addressed to assure identification 
of an optimum design. 

92. The second column of Table 7 includes choices for structural in- 
tegrity criteria in terms of the damage to the armor layer, as might be esti- 
mated by Equation 14. The hD(Hq) value chosen should be consistent with 
the incipient damage level, as measured in model experiments pertinent to 
the breakwater design at hand. Hy, is the wave height applied in analyti- 
cal stability relations. Return periods of 30 years or less for the storm 
represented by Hq will plainly involve substantial expected damage and 
therefore should be investigated only for minor breakwaters where repairs can 
be easily accomplished or postponed without significant adverse consequences. 
Long return periods greater than 50 years are important to address, however, 
since rubble-mound breakwaters require such a tremendous commitment of equip- 
ment and materials to repair. The risk of affordable quarrystone being un- 
available 30 or 40 years in the future might be great, even though it may be 
readily available at present. Repair of breakwaters in remote areas involves 


high mobilization and demobilization costs, even for small repair efforts. 
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Another important consideration in favor of addressing these longer return 
periods is the uncertainty of the future funding capacity of local sponsors 
for repair efforts. 

93. The final choice of alternatives should contain a minimum of 
15-20 pairs of functional performance and structural integrity criteria pairs. 
A single pair of these criteria will define each alternative breakwater con- 
figuration throughout the optimization process. Consistency in application of 
these criteria in analytical design efforts is critical to maximizing the re- 
liability of the procedure. New alternatives should not be added without car- 


rying the new ones through the entire procedure. 


Step 4-identify apparent optimum com- 
bination of armor size and type, slope, 
and crest elevation for each alternative 


94. Each pair of design criteria will have several combinations of 
features that will provide the same performance and stability. An acceptable 
method of choosing an apparent cost-effective combination for each plan is to 
consider a standard parameterized cross section, as illustrated in Figure 17. 
The Hudson formula (Equation 1), the relation of armor thickness and crest 
width to armor weight (Equation 9), and the wave transmission relations (Equa- 
tions 26 through 29) can then be used to approximate all the dimensions of 
this standard cross section for a range of armor type and slope combinations. 
The relative advantages of armor unit hydraulic stability and of runup dissi- 
pation are both measured by this approach. The relative cost per unit length 
of breakwater trunk for each slope and unit type combination can also be esti- 
mated by incorporating representative unit prices for each armor type and 
size. 

95. This method does not deal with the variation of reserve stability 
between armor types which would involve a substantial amount of extra input 
and computational effort. The question of reserve stability is addressed 
later in this proposed procedure, but at this stage it is neglected as a 
time-saving measure. "BWCOMP," an interactive computer program, has been 
developed to estimate the volume and first cost per unit trunk length of the 
parameterized cross section of Figure 17, given the two design criteria (as 
incident H, and Tp values and an H* maximum transmitted height) 
along with the other information discussed above. The program is documented 
in Appendix B of this report and in a CETN (WES, CERC 1984b). 
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Step 5--design detailed 
cross section for each alternative 

96. This is the second highly subjective step in the proposed optimi- 
zation procedure where coastal engineers should, for each pair of design cri- 
teria, prepare a cross-section design with all the detailed features appropri- 
ate for the site conditions and other constraints. Practical considerations 
discussed in Part II of this report should be incorporated. All the special- 
ized experience and intuition available should be applied in this step, but it 
must be applied consistently to each alternative. It is critical that bias be 
studiously avoided at this stage. An estimate of the construction cost for 
each alternative detailed cross section should be prepared at the conclusion 
of this step. 
Step 6--estimate wave transmission 
characteristics of each alternative 

97. An analytical procedure should be performed at this point to esti- 
mate the wave transmission characteristics as a function of incident waves 
H,(H; ) for each alternative. The program MADSEN (Seelig 1980a and WES, CERC 
1984a) is useful for this purpose. The program accounts for the relative size 
and permeability of each layer of the breakwater cross section and the rela- 
tive runup characteristics of the armor layer. Wave transmission by overtop- 
ping (Equations 26 through 29) and permeation (Madsen and White 1976) is es- 
timated. The program is not as well verified for concrete armor units as for 
quarrystone, but it serves well at this stage for comparative purposes. A 
range of incident wave conditions should be simulated to obtain a substantial 
set of H, (H; ) points, including several more severe than the design condi- 
tion. The incident wave conditions need to correspond to height and period 
combinations predicted for the site in Step 1. Wave period is a sensitive 
factor for wave transmission, as applied in the program MADSEN. An appropri- 
ate wave period (such as the peak spectral period Ty) must therefore be as- 
sociated with each (significant) incident wave height, as suggested in Step 1. 
Transmitted waves are not Rayleigh distributed, as discussed in Part IV and 


Andrew and Smith (in preparation), but can be represented by a single height 


such as the root mean square wave height ans or H43 54° MADSEN predicts 
the Hrms Of waves transmitted by the combined effects of both permeation and 
overtopping. 
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Step 7--estimate economic losses 
with the breakwater for each alternative 

98. The climate of transmitted waves behind the breakwater can now be 
approximated as a cumulative probability distribution F(H,) given a set of 
H, (A; ) points from Step 6 and the cumulative distribution of incident waves 
F(H; ) from Step 1. The loss function estimated in Step 2 can be used to es- 


timate the expected annual economic losses E{$L'/yr} for each alternative by 


' dF(H, ) 
e{ $"| Siig) SLC) Uda (33) 


yr dH, t 
"BWLOSS2," a computer program, has been developed to perform these computa- 
tions. It has been documented in a CETN*, and it is included in Appendix C of 


this report. 


Step 8--estimate expected annual 
breakwater damages for each alternative 
99. The methods discussed in Part III can be applied to relate a damage 


funetion %D(H/Hq) to each alternative. The incident wave climate defined 
by F(H) from Step 1 can in turn be applied to estimate the expected annual 
damages E{$D/yr} given representative unit repair prices $/vol and the 


volume of the armor layer Vol by adapting Equation 18 as follows: 


p22 | yore (5) Sa ee (34) 


lyr vol dH 


100. This quantity is useful for comparative purposes, but it does not relat 

directly to a programmed cash flow for repairs. It is better that Equation 18 

be applied to each alternative in its unmodified form to predict the expected 
annual %D in order to make some judgment if and when a repair project should 

be scheduled. The average time to reach a threshold level of unacceptable 

damage %D* can be estimated by simply dividing that value by E{%D/yr} . The 
return period of %D* could also be estimated by solving for H(ZD*) in the 

damage function (Equation 14) and applying Equation 21 to determine the asso- 


ciated return for that particular storm intensity. A computer program titled 


* US Army Engineer Waterways Experiment Station, Coastal Engineering Research 
Center, in preparation, "Estimation of Expected Annual Economic Losses from 
Waves Transmitted by a Breakwater--Computer Program BWLOSS2 (MACE-16) ," 
Coastal Engineering Technical Note, Vicksburg, Miss. 
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"BWDAMAGE" has been developed which applies Equation 14 and the information 
of Table 3 to estimate E{ZD/yr} , E{$D/yr} , and the repair interval by both 
methods discussed above. This program has been documented in a CETN,* and it 
is presented in Appendix D of this report. Once a repair interval and the as- 
sociated extent of repairs have been estimated for an alternative, discounted 
cash flow methods can be used to estimate the equivalent annual amount which 

can be substituted for E{$D/yr} . The damage functions, as stated in Part 
| III, are currently the least reliable of the analytical tools available for 
rubble-mound breakwater design and should be used with circumspection. 
Step 9--tabulate and sum 
costs for each alternative 

101. This is the final analytical step of the proposed procedure, fol- 

lowed only by laboratory verification of the analytical predictions. The min- 
imum sum of the three costs identifies the cost-effective optimum alternative, 


as indicated in the following equation: 


' st 
| Zz S| el g{ $2 (35) 
yr yr 
The first cost must be transformed from a present worth value to an equivalent 
annual amount E{$18* yr} by discounting: prior to the summation. Incremental 
benefits E{$B/yr} can be estimated by subtracting E{$L'/yr} from 
E{$L/yr}: 
( 
e{ 8} - oft - of Eh (36) 
yr yr yr 
Net benefits E{$B,.,/yr} can in turn be estimated by subtracting E{$18* yr} 
and E{$D/yr} from E{$B/yr} as follows: 


oft] = BB) ofA] - of) an 


This method of estimating benefits may not be appropriate for some projects, 


however, as discussed at the beginning of Part V. 


* US Army Engineer Waterways Experiment Station, Coastal Engineering Research 
Center, in preparation, "Estimation of Rubble-Mound Breakwater Expected 
Damages--Computer Program BWDAMAGE (MACE-18) ," Coastal Engineering Technical 
Note, Vicksburg, Miss. 
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Step 10--verify by physical 
modeling damages and wave transmis- 


sion of apparent optimum alternative 


102. This step is necessary to assure that all the compounded assump- 
tions and analytical inaccuracies are within acceptable limits. This is the 
case with any analytical design procedure for rubble-mound breakwaters since 
the empirical relations have been shown to all have limited confidence. Each 
laboratory test of analytical assumptions applied to a specific design will 
narrow the confidence limits and improve the reliability of future analytical 
efforts. A simple proof test with monochromatic waves of varying period con- 
stitutes a minimum effort in this direction, but it is inadequate to test the 
accuracy of an optimization procedure such as that proposed above. 

103. The damage function %D(H/Hq) must be verified by model testing, 
including simulation of conditions for incipient motion and a range of more 
severe conditions. The design conditions should be simulated as accurately as 
possible in order to include the effects of the numerous physical parameters 
not explicit in the analytical stability formula that was applied. Wave pe- 
riod, wave groupiness, storm duration, and static stability, among other fac- 
tors, should be considered. The static friction factor yw from the Iribarren 
formula (Equation 2) should be measured by sliding tests, as proposed by Price 
(1979) and Graveson, Jensen, and Sorensen (1980). 

104. The fully described incident’ wave conditions cannot be simulated 
with monochromatic waves. Either an average (for example JONSWAP) spectral 
shape or one adjusted to be similar to measured spectra for extreme storms 
near the site can be applied in flume tests of the apparent optimum cross sec- 
tion. Simulation of a gradual rise to peak conditions, then 1,000 waves or 
more at the peak (stability criterion) condition, followed by a gradual de- 
crease of wave energy, would be most useful for tests to verify damage func- 
tions. A test or tests at the design condition should be followed by tests at 
more extreme conditions related to the extremal distribution of wave heights 
(and periods) derived in Step 1. Enough %D(H/Hq) points must be measured 
to verify or refine the %D(H/Hy) analytical function that was applied in 
Step 8. A minimum of three tests would be useful, including the mD(Hq) 
point and at least two more severe conditions. More stability tests should 
be conducted if agreement with the predicted damage function is not good. 


Techniques to detect gross rocking motion should be applied in identifying 
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incipient motion. Actual damage should be measured by before and after sound- 
ings on a fine grid, but some judgment must be made as to the additional dam- 
age that might have occurred in prototype from armor unit breakage. 

105. Wave transmission characteristics of the apparent optimum cross 
section must also be verified. Tests of these design conditions simulating 
the fully described forecast conditions at the site as accurately as possible 
should be performed for the functional performance criteria and a number of 
more extreme conditions. Again, these extremes should relate to the F(H, A 
Ty) derived in Step 1. At least three Hee (Hoa) points should be measured 
in order to verify or refine the economic loss function derived in Step 7. 
Operational techniques should include efforts to accurately model reflection 
and wave transmission by both overtopping and permeation. Transmitted waves 
should be measured as time series comparable to time series measured of inci- 
dent waves. Coherence and cross-correlation analyses should be performed for 
the incident and transmitted time series along with computation of more common 
spectral parameters. Individual runs of 100 or more waves are recommended for 
the wave transmission tests in keeping with the widely accepted assumption of 
stationarity in natural sea states. 

106. The measured %D(H/H,) data and H,(H;) data should be applied 
in Steps 6 through 9 for the apparent optimum cross section. All its associ- 
ated costs should then be adjusted according to the revised expected damages 
and economic losses with the breakwater in place. Model tests often make sig- 
nificant refinements to a design cross section obvious, and any such refine- 
ments should be incorporated. Drastic changes to the original apparent opti- 
mum cross section may require similar changes to be made to all the alterna- 
tives and for Steps 5 through 9 to be repeated for these cross sections as 
well. If the original apparent optimum is still indicated as the optimum 
cross section, then no further model testing will be necessary. A new appar- 
ent optimum should have its #D(H/Hg) and H, (H; ) functions verified in as 


thorough a manner as the first. 
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PART VI: SUMMARY AND CONCLUSIONS 


Summary 


Optimization 

107. Optimization has been demonstrated as a systematic process of max- 
imizing net tangible economic benefits or of minimizing the total costs (in- 
cluding economic losses) to the beneficiaries of a public works project. Op- 
timization of rubble-mound breakwaters addresses the incremental net benefits 
of these structures which are often major features of a larger coastal devel- 
opment. Federal laws and policies currently require that incremental net ben- 
efits be positive for all major features of projects proposed for federal 
funding. Furthermore, cost sharing policies have placed a substantial burden 
for financing these projects on local and regional governments. Financeabil- 
ity of civil works projects is now an important question outside that of posi- 
tive net benefits. Rubble-mound breakwaters must achieve the maximum benefits 
for the least cost in order to be affordable as well as economically feasible. 
Arbitrary conservatism in design of rubble-mound breakwaters is no longer 
affordable, and coastal engineers must use all the tools and information 
available to assure the optimum alternative has been proposed. 
Design criteria 

108. Alternatives for rubble-mound breakwaters should be optimized ac- 
cording to two criteria: functional performance and structural integrity. 
The functional performance criterion refers to the structure's effectiveness 
as a wave barrier as measured by its wave transmission characteristics. The 
structural integrity criterion refers to the structure's ability to survive an 
extreme storm without significant damage and the rate it suffers damage from 


storms more extreme (less probable) than the structural design event. 


Analytical design and 
laboratory verification 
109. The analytical tools available to designers of rubble-mound break- 


waters have been reviewed in some detail. They have all been shown to be the 
products of a finite set of laboratory experiments, with very little quantita- 
tive prototype verification. Current research continues to refine the preci- 
sion of these empirical relations, but this precision is not yet sufficient to 


warrant construction of rubble-mound breakwaters without verification of 
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analytical predictions by scale model tests. Nevertheless, analytical proce- 
dures are available for prediction of armor unit hydraulic stability (resis- 
tance to displacement by waves), armor layer damage rates, and breakwater wave 
transmission characteristics. These tools, with laboratory verification, can 
be used to systematically select an optimum alternative. 
The proposed procedure 

110. A systematic optimization procedure has been proposed which makes 
use of the analytical tools currently available to coastal engineers for 
rubble-mound breakwater design. The procedure begins with definition of the 
site conditions and formulation of an ensemble of alternative design criteria 
pairs. These steps are followed by estimates of first costs, maintenance 
costs, and user costs with the breakwater in place for each alternative. The 
concept of statistical expectation is applied to measure the costs of all al- 
ternatives on the same basis. The process is concluded by physical model 
tests to verify the analytical predictions for structural stability and wave 
transmission characteristics of the apparent optimum alternative. The entire 
procedure is summarized in Table 8, with references to pertinent formulae, 


software, and documentation. 


Table 8 
Summary of Optimization Procedure 


Pertinent Available 
Step Procedure Equations and Tables Software 
| Define site Equations 19* or 22* WAVDIST (WES, 
conditions CERC (in 
preparation) ) 
2 Estimate economic Equations 31* and 32* BWLOSS1 (WES, 
losses without CERC (in 
breakwater preparation) 


and Appendix A) 


3 Formulate an en- Table 6 -- 
semble of alterna- 
tive functional and 
structural criteria 
pairs 


(Continued) 


Note: * indicates the equations which are applied in the referenced 
software. 
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Procedure 


4 Identify optimum 
armor, type W cot 
@ , and crest ele- 
vation for each 
alternative 


5 Design detailed 
cross section for 
each alternative 


6 Estimate wave 
transmission char- 
acteristics of each 
alternative 


if Estimate economic 
losses with break- 
water for each 
alternative 


8 Estimate breakwater 
damages for each 
alternative 


9 Tabulate expected 
costs for each 
alternative and 
identify apparent 
optimum 


10 Verify predicted 
damage and wave 
transmission by 
scale modeling 


Table 7 (Concluded) 


Pertinent 
Equations and Tables 
Equations 1* (or 2-8), 
29* (or 30), and 
Table 5* or 6 


Equations 1-9 


Equations 19 or 22, 
26*, 27%, 287, and 29* 
(or 30) and Table 5* 
or 6 


Equation 19% or 22,13 1*; 
and 33 


Equations 19* or 22, 14%, 
(or (5,5. °65. (7 sxand )8.c0r 
23 and 24), 18*, and 34* 
and Table 3* 


Equations 35, 36, and 37 


Equations 1-8, 10-14, 
and 26-29 


Available 
Software 


BWCOMP (WES, CERC 


1984b and Ap- 
pendix B) 


MADSEN (Seelig 
1980a and WES, 
CERC 1984a) 


BWLOSS2 (WES, 
CERC (in prep- 
paration) and 
Appendix C) 

BWDAMAGE (WES, 
CERC (in prep- 
aration) and 
Appendix D) 


Uthke 


Conelusions 


The investigation which was conducted in order to develop the 


above optimization procedure led to the following conclusions regarding 


rubble-mound breakwater design: 


a. A systematic optimization procedure should be applied in any 


rubble-mound breakwater design to assure that an alternative 


with maximum cost effectiveness is proposed. 


7) 


Io 


Ke) 


Rubble-mound breakwater designs should not be constructed 
without physical model testing of some kind due to the limited 
confidence of available analytical methods. 


The confidence of the key analytical tools for rubble-mound 
breakwater design would be improved if current research were 
continuously concentrated in the following specific areas with 
probabilistic applications in mind: 


(1) Site conditions--Estimation of the long-term joint proba- 
bility distribution .F(H:, T , t., d 9.eé¢so0)e0 for vay site 
should be developed for application in estimating ex- 
pected breakwater damages and the long-term distribution 
of transmitted wave characteristics. 


(2) Armor stability--Standardized methods should be developed 
for scale model testing of rubble-mound stability in nat- 
ural irregular sea states. Improved analytical stability 
prediction should be the goal of tests conducted by these 
methods, explicitly including the effect of wave period, 
storm duration, and other factors. Prototype verifica- 
tion of analytical predictions should be attempted also, 
particularly for new constructions where the design as- 
sumptions are most thoroughly documented. 


(3) Mechanical stength of armor units--Prediction of armor 
unit breakage by scale model tests should be developed in 


order that both incipient damage and reserve stability 
can be more accurately defined. 


(4) Breakwater damage prediction--The reserve stability of 
a wide range of rubble-mound breakwater configurations 


should be comprehensively tested by methods similar to 
those developed to detect incipient damage. Improved 
analytical prediction of reserve stability should be the 
goal of these tests. 


(5) Runup on rubble-mound breakwaters--Improved instrumen- 
tation and testing methods need development for measure- 


ment of irregular runup on rough permeable slopes. A con- 
certed effort should be made to define runup coefficients 
for Equations 31 and 32 while concurrently investigating 
means for improved analytical prediction of irregular 
runup. The possibility of armor units designed both for 
enhanced hydraulic stability and for efficient attenua- 
tion of runup should be explored. 


(6) Wave transmission--The characteristics of irregular 
waves transmitted by rubble-mound breakwaters should be 
investigated. Improved analytical prediction of trans- 
mitted wave characteristics as a function of incident ir- 
regular wave characteristics should be the goal of this 
research. 
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APPENDIX A: COMPUTER PROGRAM BWLOSS1 


Estimation of Economic Losses as a Function of Wave Height 


Program purpose 


1. The program BWLOSS1 is intended to aid planners of coastal struc- 
tures which provide protection from wave attack by deriving an empirical math- 
ematical expression relating a given level of economic losses to the responsi- 
ble incident significant wave height. This loss function can be used to de- 
fine the "without-project" condition with respect to the incremental economic 
benefits provided by artificial wave protection. The program optionally pro- 
vides an estimate of expected annual economic losses due to wave attack, given 
the coefficients of an Extremal Type I cumulative probability distribution 
function of significant wave heights for the site. 

Program capabilities 

2. BWLOSS1 is written in FORTRAN IV as implemented on the Honeywell 
DPS-8 mainframe system at the US Army Engineer Waterways Experiment Station 
(WES). A BASIC version written for the IBM PC is also available. The least 
squares method is applied to historical data on economic losses associated 
with the significant wave heights of the storms that caused the losses. A 


loss function is derived from the following form: 


A(H_-H._) 
$L(H,) = Li 23 3 Lo? | (A1) 


$L(H,) = economic losses as a function of significant wave height H, 


= maximum conceivable economic loss from wave attack (at any 
intensity) 


A = site-specific coefficient derived by regression 


Hp = maximum significant wave height for which economic losses are 
negligible 
3. The regression requires at least one point for H, , $L(H,) » but 


it can deal with up to 100. The coefficient A is presented along with the 
nonlinear correlation coefficient and the sum of the square residuals. A ta- 
ble of residuals is optionally presented. Losses can be optionally predicted, 


given a specific significant wave height, or the significant wave height 


Al 


corresponding to a given level of losses can be predicted. The form of this 
function is illustrated in Figure 16 in the main text. The program will also 
apply an Extremal Type I cumulative probability distribution of significant 


[e-H,/¢] 


wave heights as follows: 


F(H_) = e617 
(A2) 
“where 
F(H,) = cumulative probability distribution of events 
Where ~H_ < H 
s Ss 
e and » = site-specific coefficients derived by regression of 
historical wave data 
to estimate the expected annual economic losses by 
H dF(H_) 
$L\_ So Ss 
e\t = a $L(H,) ai. dH, (A3) 


Lo 


where 


A} = the average number of extreme events per year above the 
threshold H. value originally used to derive e« and 9 
(must be input by the user) 


H = a practical upper limit taken as the H value whose 
probability of exceedance is 0.0000001 


4, This formulation assumes that the number of extreme events per year 
is random and can be represented by a mean value and is independent of the 
Significant wave heights representing the intensity of the individual storms. 
The lower limit of integration is H;, , below which the expected losses are 
taken as zero. Extrapolation of F(H,) to H, values below the threshold 
value applied to data used to originally derive e and 9 is probably con- 
servative, but this question will be the subject of further study. A thresh- 
old H, value set equal to H;, would presumably resolve any problems if 
adequate statistical confidence can be maintained. The integration is accom- 
plished by a numerical application of Simpson's Rule with 100 intervals. 

5. The majority of the expected losses statistically occur during 
storms whose H,; is just above H,, where the probability density is sub- 
Stantial. The higher H, values occur on the tail of the probability density 
function and may even be precluded by depth limitations. The program does not 


deal with depth limitations and assumes the Extremal Type I function fully 


A2 


represents the wave climate at the site. A potential improvement of BWLOSS 1 
is the incorporation of period and depth effects for an estimate of $L(H, , 
Tp , d) given the joint probability distribution FHS P Tp wd)! ov heeurther 
improvement would also incorporate the storm duration t for an estimate of 
$L(H, i Ty ; tsa) , given GES ; Ty , t , d) . These enhancements will 
involve a much more rigorous computation than is now performed by BWLOSS1. 

The program is now completely interactive and easily adaptable to execution by 


microcomputer systems. 


Sample Execution and Output 


6. Below is a sample execution and output for computer program 
"BWLOSS1." 


INPUT THE MAXIMUM CONCEIVABLE LOSS IN MILLIONS OF DOLLARS 
=20. 

INPUT THE MAXIMUM SIGNIFICANT WAVE HEIGHT FOR WHICH 
LOSSES ARE NEGLIGIBLE - USE CONSISTENT UNITS 

=2, 

HOW MANY SIGNIFICANT WAVE HEIGHT VS LOSS 

DATA POINTS DO YOU HAVE? 

=4 


ENTER SIGNIFICANT WAVE HT.,COMMA,LOSS IN MILLIONS OF DOLLARS 
AND RETURN FOR EACH POINT 


23.,.5 
=4.5,1. 
#605205 
winavi 5 


DATA ON EXPONENTIAL CURVE... 
CURVE HAS FORM: $L(Hs)=$Lmax#{1-explA#(Hs-HLo) J} 


$Lmax= 20. 0000000 

HLo= 2.00008 

A= -8.8437137 

$L(Hs)= LOSSES 

Hs= SIGNIFICANT WAVE HEIGHT 
NON-LINEAR CORRELATION IS @.9735292 
SUM SQR RESIDUALS..... 1.9485841 


A3 


PRINT RESIDUAL TABLE(Y/N)? 
=Y 


XVALUE YVALUE YEST DIFF 

2.0008 Q. 0.2000 0.0000 
3.0008 0.5800 0.8554 8.3554 
4.5008 1.0008 2.0785 1.0785 
6.0080 2.5008 3.2084 0.7084 
12.0000 7.5000 7.0823 0.4177 


DO YOU WANT TO MAKE SOME LOSS PREDICTIONS 

FROM SIGNIFICANT WAVE HEIGHT DATA(Y/N)? 

=y 

INPUT SIGNIFICANT WAVE HEIGHT 

=10. 

PREDICTED LOSS IN MILLIONS OF DOLLARS IS 5.98 
DO YOU WISH TO MAKE ANOTHER PREDICTION(Y/N)? 

=N 


DO YOU WANT TO PREDICT SIGNIFICANT WAVE HEIGHTS 
FROM LOSS DATA(Y/N)? 

zy 

INPUT LOSS IN MILLIONS OF DOLLARS 

215. 

PREDICTED SIGNIFICANT WAVE HEIGHT IS 33.71 

DO YOU WISH TO MAKE ANOTHER PREDICTION(Y/N)? 

=N 

DO YOU WANT TO PREDICT EXPECTED ANNUAL LOSSES(Y/N)? 
=y 

SELECT A DISTRIBUTION... 

EXTREMAL TYPE I...1 

WETBUL( catdokacaaiee 2 

LOG-EXTREMAL......3 

SELECT 1, 2, OR 3 

=1 


INPUT EXTREMAL TYPE I EPSILON, AND PHI 
=-2,27,3.216 


INPUT AVERAGE NUMBER OF EXTREMAL EVENTS PER YEAR, 

THE POISSON ‘LAMBDA’ PARAMETER 

=4 

EXPECTED ANNUAL LOSS IN MILLIONS OF DOLLARS IS 2.4141522 


J 


AY 


Program Listing 


7. Below is a program listing for computer program BWLOSS1 (FORTRAN 
version). 


10@C PROGRAM "BWLOSSI". 11/85 VERSION 

20C DESIGN BRANCH-COASTAL ENGINEERING RESEARCH CENTER 
3OC U.S. ARMY ENGINEERS WATERWAY EXPERIMENT STATION 
4ec P. 0. BOX 631 

98C VICKSBURG, MS 39180-0631 

60C FOR FURTHER INFORMATION CONCERNING THE APPLICATION 
7OC OF “BWLOSS1", CALL.. 


BOC ORSON P. SMITH (601)-634-2013 FTS:542-2013 OR 
98C ROBERT B. LUND (601)-634-2068 FTS:542-2068 OR 
1@@C DOYLE L. JONES (601)-634-2069 FTS:542-2069 

11@C 

120C FORTRAN 4 HONEYWELL DPS-8 


130C REF: "COMPUTER PROGRAM WAVDIST" CETN-I- 
140C REF: "PROBABILITY AND STATISTICS" BY MORRIS DEGROOT 
158C REF: "COST EFFECTIVE OPTIMIZATION OF RUBBLE-MOUND BREAKWATER 


16@C CROSS-SECTIONS" BY ORSON P. SMITH 

170C REF: "EXTREMAL STATISTICS IN WAVE CLIMATOLOGY" BY BORGMAN AND RESIO 

188C 

19@C N THE NUMBER OF DATA POINTS 

208C X THE ARRAY OF SIGNIFICANT WAVE HEIGHTS 

218C YH THE ARRAY OF LOSSES CORRESPONDING TO EACH SIGNIFICANT WAVE HEIGHT 


THE TRANSFORMED Y ARRAY USED IN THE METHOD OF LEAST SQUARES 
Hlo, THE MAXIMUM WAVE HEIGHT FOR WHICH LOSSES ARE NEGLIGBLE 
$Lmax, THE MAXIMUM CONCEIVABLE LOSS IN MILLONS OF DOLLARS 
A, THE REGRESSION COEFFICIENT A<Q 

THE NON-LINEAR CORRELATION OF THE LOSS FUNCTION 

THE SUM OF THE SQUARE RESIDUALS 

THE ARRAY THAT CONTAINS THE RESIDUAL FOR EACH DATA POINT 
THE LOSSES AS ESTIMATED BY THE LOSS CURVE 

THE PROBABILITY DENSITY FUNCTION OF EXTREMAL WAVES 

THE CUMULATIVE DISTRIBUTION FUNCTION OF EXTREMAL WAVES 


nN 
oe 
Ss 
a 
(sp 
oO 
pe) 
ps) 
ow 0 8 ob th ob eb ow ow we oe 


328C G(X) THE LOSS FUNCTION 

338C C = SIMPSON'S RULE COEFFICIENTS 

348 

350C INITIALIZE VARIABLES,STRINGS,AND FUNCTIONS 

360 DIMENSION X(101),Y(101) ,YH(1@1) 

370 COMMON X,Y,YH 

388 CHARACTER#*1 L 

398 CHARACTER*46@8 ST(20) 

400 G(X) =We(L-EXP(V2%(X-V1))) 

410 ST (1) = RRR KE RHEE HERE HERR EHR HHH R KKH RHEE HEHEHE EHH EEE EHH H ER HE 
420 ST(2)='* "BWLOSS1" IS A PROGRAM WHICH FITS AN EXPONENTIAL * 
430 ST(3)='* CURVE TO HISTORICAL INFORMATION ON ECONOMIC LOS- * 
440 ST(4)=’* SES CAUSED BY WAVE ATTACK. EACH STORM CAUSING * 
458 ST(5)='* LOSSES IS ASSUMED TO BE CHARACTERIZED BY A SIN- * 
460 ST(6)='# GLE SIGNIFICANT WAVE HEIGHT. THE PROGRAM RE- + 
4708 ST(7)='* QUIRES ESTIMATES OF THE MAXIMUM LOSS SUSTAINABLE * 
482 ST(8)='* FROM WAVE ATTACK,THE MAXIMUM SIGNIFICANT WAVE * 
499 ST(9)='* HEIGHT FOR WHICH LOSSES CAN BE NEGLECTED AND AT * 
588 ST(10)='* LEAST ONE HISTORICAL LOSS WITH ASSOCIATED WAVE 

518 ST(11)='* HEIGHT. THE EXPONENTIAL CURVE IS COMPUTED BY * 


A5 


520 ST(12)='* THE LEAST SQUARES METHOD. ITS PARAMETERS AND * 
530 ST(13)='* NON-LINEAR CORRELATION ARE PRINTED. P 
540 ST(14)='4 THE PROGRAM WILL ALSO ESTIMATE EXPECTED ANNUAL # 
550 ST(15)='* LOSSES GIVEN THE PARAMETERS FOR THE LONG-TERM r 
560 ST(16)='* CUMULATIVE PROBABILITY DISTRIBUTION OF SIGNIFICANT +# 
570 ST(17)='* WAVE HEIGHTS AT THE SITE. THE PROGRAM WILL AC- * 
580 ST(18)='* CEPT THREE DIFFERENT DISTRIBUTIONS: (1) EXTREMAL * 
590 ST(19)='* TYPE I; (2) WEIBULL; AND (3) LOG EXTREMAL. P 
600 $1 (20) =ST(1) 

610 DO 5@ 1=1,20 

620 WRITE(6,407) ST(I) 

630 407 FORMAT(1X,A6Q) 

640 50 CONTINUE 

650 

66@C GET THE FACTS 

670 WRITE (6,408) 

688 408 FORMAT(///) 

690 1 WRITE (6,181) 

700 101  FORMAT(1X,"INPUT THE MAXIMUM CONCEIVABLE LOSS IN MILLIONS OF DOLLARS") 
710 READ ,W 

720 IF(W .LE. @) GO TO 1 

730 2 WRITE (6,201) 

740 201 FORMAT(1X,"INPUT THE MAXIMUM SIGNIFICANT WAVE HEIGHT FOR WHICH",/,1X, 
750 & "LOSSES ARE NEGLIGIBLE - USE CONSISTENT UNITS") 

760 READ, X(1) 

778 IF(X(1) .LT. @) GO TO 2 

780 YH(1)=0 

790 ¥(1)= 

Baa 4 WRITE (6,102) 

810 182 FORMAT(1X,"HOW MANY SIGNIFICANT WAVE HEIGHT VS LOSS",/,1X, 

820 & "DATA POINTS DO YOU HAVE?") 

830 READ,N 

840 IF(N .LE. 1) GO TO 4 

850 IF( N .GT. 108) PRINT, ‘10@ POINTS IS MAXIMUM-REINPUT’ 

868 IF( N .GT. 108) GO TO 4 

870 8 WRITE (6,104) 

888 104 FORMAT(/,1X,"ENTER SIGNIFICANT WAVE HT. ,COMMA,LOSS IN 

890 & MILLIONS OF DOLLARS",/,1X,"AND RETURN FOR EACH POINT") 

900 1=2 

910 15 READ, X(1),VH(T) 

920 IF( YH(I) .GT. W) GO TO 17 

930 ¥(1)=ALOG(1-YH(I) /W) 

940 IF(I .EQ. (N+1)) GO TO 18 

950 l=I+1 

960 60 TO 15 


978 17 WRITE (6,185) 
98@ 105 FORMAT(/,1X,"ERROR-YOUR INPUT LOSS IS MORE THAN YOUR MAXIMUM" 


990 & ,/,1X,"LOSSES. RE-INPUT POINT") 
1000 60 TO 15 

1010 18 N=N+li 

1020 

1@30C FIT CURVE TO INPUT DATA 

1040 CALL LOG(N,W,V1,V2) 


A6 


1050 
1068 30 WRITE (6,110) 
1070 110 FORMAT(//,1X,"DO YOU WANT TO MAKE SOME LOSS PREDICTIONS ",/,1X, 


1080 & "FROM SIGNIFICANT WAVE HEIGHT DATA(Y/N) ?") 
1090 CALL ANS(L) 
1100 TRC SEG. ON) GO 10) 775 


11108 45 PRINT, INPUT SIGNIFICANT WAVE HEIGHT’ 
1128 46 READ,H 


1138 IF( H .LT. @) GO TO 45 
1140 SLOG=G(H) 

1158 IF( H .LE. X(1) ) SLOG=0 
1160 WRITE(6,114) SLOG 


1178 114 FORMAT(1X,"PREDICTED LOSS IN MILLIONS OF DOLLARS IS ",F7.2) 
118@ 115 PRINT, ‘DO YOU WISH TO MAKE ANOTHER PREDICTION(Y/N) ?' 

1198 CALL ANS(L) 

1200 IF(L .E@. ‘¥‘) GO TO 45 

1210 

1220 

1230C FIND SIGNIFICANT WAVE HEIGHT GIVEN DAMAGE 

124@ 75 WRITE (6,120) 

125@ 128 FORMAT(/,1X,"DO YOU WANT TO PREDICT SIGNIFICANT WAVE HEIGHTS", 


1268 & /,1X,"FROM LOSS DATA(Y/N)?") 

1278 CALL ANS(L) 

1288 IF(L .E@. ‘N’) GO TO 380 

1298 80 PRINT, ‘INPUT LOSS IN MILLIONS OF DOLLARS' 
1300 READ, SAB 

1318 IF(SAB .GT. W) GO TO 88 

1320 WHT=ALOG(1.Q-SAB/W) /V2+V1 

1330 WRITE(6,133) WHT 


134@ 133 FORMAT(1X,"PREDICTED SIGNIFICANT WAVE HEIGHT IS ",F7.2) 
1358 98 PRINT, DO YOU WISH TO MAKE ANOTHER PREDICTION(Y/N)?° 
1360 CALL ANS(L) 


1370 IF(L .EQ. ‘¥’) GO TO 82 

1380 

1392 300 PRINT, ’DO YOU WANT TO PREDICT EXPECTED ANNUAL LOSSES(Y/N)?’ 
1400 CALL ANS(L) 

1410 IF(L .EQ. 'N’) GO TO 400 

1420 CALL EXPCT(W,V2,X(1)) 

1438 400 STOP 

1440 END 

1450 

1460 

1470 

1480 

1490 

150@C SUBROUTINE LOG TO FIT EXPONENTIAL CURVE TO INPUT DATA 
1510 SUBROUTINE LOG(N,W,V1,V2) 

1520 CHARACTER#1 L 

1530 COMMON X,Y, YH 

1540 DIMENSION X(1O1),Y(1O1),YH(101),D1(101),Z1(101) 
1558 G(X) =We (1. O-EXP(V2*(X-V1))) 


15608C CALCULATE PARAMETERS V1 AND V2 BY THE LEAST SQUARES METHOD 
1578 YSUM=0 


A7 


1588 DT=0 


1590 DB=0 

1600 DO 20 K=2,N 

1610 YSUM=YSUM+YH(K) 

1620 DT=DT+Y (K) #(X(K)-X(1)) 

1630 20 DB=DB+(X(K)-X(1)) ##2 

1640 Vi=X (1) 

1650 V2=DT/DB 

1660C CALCULATE CORRELATION COEFFICIENT (R#*2) 
167 YAVG=YSUM/N 

1680 ST=0 

1690 SB=0 

1708 DO 78 I=1,N 

1710 11(1)=6(X(1)) 

1720 D1(1)=(YH(1)-24 (1) ) ##2 

1730 D2=(YH(1)-YAVG) #2 

1740 ST=ST+D1 (1) 

1750 70 SB=SB+D2 

1760 IF( (1-ST/SB) .LT. @) CORR=1-ST/SB 
1770 IF( (1-ST/SB) .LT. @) GO TO 157 
1780 CORR=SORT(1.0-ST/SB) 


1798C PRINT OUT PARAMETERS AND OTHER DATA 
1800 157 WRITE (6,80) 


1810 8 FORMAT(//,1X,"DATA ON EXPONENTIAL CURVE...",/,1X, 

1820 & "CURVE HAS FORM: $L(Hs)=$Lmax#{(1-explA®(HS-HLo) 1}") 

1838 WRITE(6,82) W,V1,V2 

1848 82 FORMAT (1X, "$Lmax=",F14.7,/,1X,"HLO=",3X,F10.4,/,1X, "A=", 4X F147 
1850 & ,/,1X,"$L (Hs) =",4X,"LOSSES",/,1X,"Hs=",8X, "SIGNIFICANT WAVE HEIGHT") 
1860 WRITE (6,84) CORR 

1870 84 FORMAT(1X, "NON-LINEAR CORRELATION IS",3X,F9.7) 

1888 WRITE (6,86) ST 

1890 86 FORMAT(1X,"SUM SQR RESIDUALS.....",4K,F11.7) 

1900 35 WRITE (4,37) 

1910 37 FORMAT(////,1X,"PRINT RESIDUAL TABLE(Y/N)?") 

1920 CALL ANS(L) 

1930 TFC eb. “Ys Kot 

1948 IF( L .NE. ‘N’ .AND. L NE. ‘Y‘) GO TO 35 

1958 IF( K .EQ. 1) WRITE(6,45) 

1960 45 FORMAT(///,1X," XVALUE YVALUE YEST DIFF") 
1978 DO 6@ I=1,N 

1980 IF( K.EQ. 1) WRITE(6,51) X(1),YH(1),Z1(1) ,SQRT(D1(1)) 

1990 51 FORMAT(F11.4,F11.4,F11.4,F11.4) 

2000 68 CONTINUE 

2010 RETURN 

2020 END 

2030 

2040 

2050 SUBROUTINE EXPCT(W,V2,CUT) 

2068 DOUBLE PRECISION BU 

2078 REAL LAMBDA 

2080 FD1(X)=-(ALOG(-ALOG(X)) #PHI)+EPSI 


A8 


2098 FD2(X)=((-ALOG(1-X)) ##(1/A1)) #BI 


2100 FD3(X)=B2/((-ALOG(X)) ##(1/A2)) 

2110 PDF1(X) SEXP (-EXP(-(X-EPSI) /PHI)) #EXP (-(X-EPSI)/PHI) /PHI 
2120 PDF2(X)=AL#(X¥##(AL-1)) FEXP(-(X/B1) ##AL) /(BLERAL) 
2130 PDF3(X) =A2#(B2##AZ) €EXP (-(B2/X) ##AZ) /(X#* (A241) ) 
2140 CDF1(X) =EXP(-EXP( (EPSI-X) /PHI) ) 

2150 CDF2(X)=1.0-EXP(-(X/B1) ##A1) 

2160 CDF3(X) =EXP (-(B2/X) ##A2) 

2170 G(X) =W#(1. O-EXP(V2*(X-CUT))) 

2180 78 PRINT, ‘SELECT A DISTRIBUTION... ’ 

2190 PRINT, EXTREMAL TYPE I...1' 

2200 PRINT, “WETBULLe ans cecneced” 

2210 PRINT, LOG-EXTREMAL......3' 

2220 PRINT, SELECT 1, 2, OR 3° 

2230 READ, ID 

2240 IF(ID .LT. 1 .0R. ID .6T. 3) GO TO 72 

2250 IF(ID .EQ. 1) WRITE(6,104) 

2260 104 FORMAT(/,1X,"INPUT EXTREMAL TYPE I EPSILON, AND PHI ") 
2278 IF(ID .EQ. 1) READ,EPSI,PHI 

2280 IF( ID .EQ. 2) WRITE(6,114) 

2290 114 FORMAT(/,1X,"INPUT WEIBULL ALPHA AND BETA") 

2300 IF( ID .EQ@. 2) READ,AL,B1 

2310 IF( ID .E€Q. 3) WRITE(6,124) 

2328 124 FORMAT(/,1X,"INPUT LOG-EXTREMAL ALPHA AND BETA") 

2330 IF( ID .EQ. 3) READ,A2,B2 

2340 WRITE (6,5) 

2350 5 FORMAT(/,1X,"INPUT AVERAGE NUMBER OF EXTREMAL EVENTS PER YEAR,", 
2360 & /,1X,"THE POISSON ‘LAMBDA’ PARAMETER") 

2370 READ, LAMBDA 

2380 BL=CUT 

2390 IF(ID .£Q. 1) BU=FD1(1-.0000001) 

2400 IF(ID .€Q. 2) BU=FD2(1-.0000001) 

2418 IF(ID .EQ. 3) BU=FD3(1-.0000001) 

2420 SUM=8 

2430 D=BU-BL 

2440 K=-1 

2450 DO 10 I=1,101 

2460 K=-K 

2470 ip Ki alte @) €=4 

2480 1F( K .GT. @ )C=2 

2490 IPC 1 SEOs i ORs lesEG. tebe! 

2500 ADD=FLOAT(I-1) *D/100.® 

2510 XV=BL+ADD 

2520 IF(ID .EQ. 1 . AND. EXP(-(XV-EPSI)/PHI) .6T. 82.0) GO TO 10 
2538 IF(TD .EQ. 2 .AND. ((XV/B1)*#AL) .6T. 82.0) GO TO 10 
2540 IF(ID .EQ. 3 .AND. ((B2/XV)##A2) .GT. 82.0) GO TO 10 
2550 IF( ID .€Q. 1) FACL=PDF1(XV) 

2560 IF( ID .EQ. 2) FAC1=PDF2(XV) 

2570 IF( ID .E@. 3) FACL=PDF3(XV) 

2580 FAC2=G(XV) 

2598 IF(XV «LT. CUT) FAC2=0 

2600 SUM=SUM+C#FACI*FAC2 

2610 10 CONTINUE 


Ag 


2628 SUM=SUM/308.0*D 


2630 WRITE(6,14) SUM#LAMBDA 

2648 14 FORMAT(/,1X,"EXPECTED ANNUAL LOSS IN MILLIONS OF DOLLARS 1S",F14.7) 
2658 RETURN 

2668 END 

2678 

2688 

2698C SUBROUTINE TO ANSWER YES/NO QUESTIONS 

2708 SUBROUTINE ANS(L) 

2710 CHARACTER#1 L 

2720 60 TO 25 

2738 30 PRINT, ‘REINPUT RESPONSE’ 

2748 25 READ(5,1@) L 

2758 18 FORMAT(A1) 

2768 IF( L .NE. ‘Y’ .AND. L .NE. 'N’ ) BO TO 30 
2770 RETURN 

2788 END 
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APPENDIX B: COMPUTER PROGRAM BWCOMP 


Comparison of Breakwater Volumes and Costs 


Program purpose 
1. The program BWCOMP calculates breakwater volumes and costs demon- 


strating the effect of varying breakwater slopes on wave transmission, the 
choice of armor size and shape, and overall volume and cost. 
Background 

2. Systematic comparison of the relative cost of rubble-mound break- 
waters designed with varying combinations of slope is tedious and awkward to 
present in project reports, yet it must be accomplished to assure that a cost- 
effective cross section is chosen. Wave transmission by overtopping during a 
given wave condition is a function of a breakwater's seaward slope, its crest 
elevation and width, and its surface roughness. A breakwater's stability is a 
function of its seaward slope, its leeward slope, the size and shape of its 
armor units, and other factors which affect the overall cost to a lesser de- 
gree. The cost of armor units varies with size and shape. The problem of 
cost comparison is further complicated by the interrelation of most of these 
factors. 

3. The program BWCOMP was designed to make the task of comparing varia- 
tions of these factors easier to accomplish and present. The following sim- 
plifying assumptions make the program economical to use: 


a. Wave transmission by permeation through the structure is typi- 
cally much smaller than transmission by overtopping and can be 
neglected for this comparative analysis. 


The unit price ($ per unit volume or unit weight) and availabil- 
ity of primary armor units tend to be most critical to the over- 
all breakwater cost and constructibility, as compared to unit 
price variations and availability of secondary armor, under- 
layers, filter material, or core material. A single average 
unit price can therefore be derived, for the purposes of this 
comparative analysis, to include all materials except the pri- 
mary armor. 


Io 


Most rubble-mound breakwaters intended primarily as wave bar- 
riers for harbors or ports must be designed for some overtopping 
during extreme events, with primary armor extending down the 
leeward slope below the water line. Final designs of rubble- 
mound breakwaters may have complex features in detail, but the 
above assumptions allow adoption of a standard parameterized 
cross section, as shown in Figure 17 of the main text. A 


Ke) 


Bl 


modification to this cross section is required for jetties 
oriented straight into oncoming waves or breakwaters with 
monolithic superstructures. 


4. The program BWCOMP uses the above assumptions to make cost compari- 
sons of alternate armor material and slope combinations with accuracy appro- 
priate for the earliest stages of planning or for Step 4 of the optimization 
procedure proposed in the main text. Final design should involve all the de- 
tailed considerations recommended in the Shore Protection Manual (SPM) (1984)# 


and other guidance available. 


Program Input 


5. The program is fully interactive in its present form and accommo- 
dates either English or metric units. It is written in FORTRAN as implemented 
on the US Army Engineer Waterways Experiment Station (WES) Honeywell DPS-8 
mainframe computer system. A BASIC version is also available for use on mi- 
crocomputers. The interactive input required is demonstrated by the example 
interactive session included in this appendix. The associated output is shown 
(in part) in the example output included in this appendix. Two wave condi- 
tions must be specified: one for determining armor size and the other for 
determining crest elevation as a function of a specified maximum transmitted 
wave height. A percentage exceedance must be associated with the specified 
maximum transmitted height such that x% of the transmitted waves can exceed 
the maximum height during the sea state represented by the second specified 
wave condition. The input unit prices (cost per volume) should be average 
values for the materials (rock or concrete) and armor unit types. This is a 
comparative analysis, so fine precision in these estimates is not necessary, 
but consistency is important. Prices that vary with the weight or volume of 
the individual armor units may require successive runs of BWCOMP since an 
estimate of these individual unit weights is necessary to input the appropri- 


ate unit price. 


* References cited in this appendix are included in the References at the end 
of the main text. 
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Computations 


6. The sequence of computations is summarized in the BWCOMP flowchart 
in Figure B1. The narrative below describes the assumptions and equations 
applied in this sequence. 

7. The program performs all computations for each of seven pairs 
of seaward and leeward slopes: 1:1.5/1:1.5, 1:2.0/1:1.5, 1:2.0/1:2.0, 
152.5/151.5% W257 1220, 1320/1215, cand 1:3;,0/ 12.0. Identical computa- 
tions are performed for each of 10 armor units for each of these slope com- 
binations. The stability, geometry, and runup coefficients which are assumed 
for each armor unit are specified in DATA statements at the beginning of the 
program listing, as summarized in Table B1. The crest elevation is first as- 
sumed as 0.3 m then increased in 0.3-m increments until the estimated trans- 
mitted wave height is less than the specified maximum. The computed dimen- 
sions and costs for all 10 armor units are then printed in a table for each 
slope combination (i.e. in seven tables). 

8. The wave conditions are checked for breaking or nonbreaking condi- 
tions by Goda's breaker index formula (Goda 1975) assuming a horizontal bot- 
tom. The stability or transmission incident heights are set equal to the 
breaker height at the specified depth if the breaker height is smaller. The 
Stability coefficient K, for Hudson's formula (Equation 1 in the main text) 
is chosen accordingly for each armor unit type. The weight computed by Hud- 
son's formula is then applied to compute the armor thickness and minimum crest 
width by Equation 9 (main text) by assuming "n" values of 2 and 3, respec- 
tively. The crest elevation derived from wave transmission computations then 
allows all dimensions of the parameterized cross section (Figure 17, main 
text) and the corresponding volumes and costs per unit trunk length to be 
estimated. Specifications from Figures 7-109 through 7-115 in the SPM (1984) 
are applied to estimate the armor thickness and number of individual armor 
units per unit trunk length. 

9. The crest elevation is determined by first assuming a crest eleva- 
tion (initially 1 ft* above the still-water level) and then estimating the 
x% transmitted wave for the specified incident wave condition and the current 


breakwater geometry. The estimated transmitted wave height is compared to the 


* To convert feet to metres, use a conversion factor of 0.3048. 
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Figure B1. 
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specified maximum, and the crest is increased 1 ft for another round of trans- 
mission computations if the condition is not satisfied. The computations ap- 


ply to Equations 26-29 and the armor unit data from Table 5 in the main text. 


Sample Interactive Session 


10. Below is a sample interactive session for program "BWCOMP" (FORTRAN 


version). 


BWCOMP IS AN INTERACTIVE PROGRAM WHICH COMPUTES BREAKWATER 
VOLUMES AND COSTS FOR A SIMPLE PARAMETERIZED CROSS SECTION FOR 
THE PURPOSE OF COMPARING THE EFFECT OF CHANGING THE SEAWARD 
AND LEEWARD SLOPES ON THE SIZE AND RELATIVE COST OF A RANGE OF 
ARMOR UNIT TYPES. THE QUTPUT OF BWCOMP SHOULD NOT BE USED 
AS A COST ESTIMATING TOOL IN ANY STAGE OF A PLANNING OR DESIGN 
PROJECT. THE ASSUMPTIONS APPLIED ARE INTENDED TO TENTATIVELY 
IDENTIFY THE OPTIMUM COMBINATION OF SLOPES AND ARMOR UNIT 
FOR A GIVEN SET OF DESIGN CRITERIA. 


TITLE OF THIS RUN? (UP TO 78 CHARACTERS) 
=EXAMPLE EXECUTION OF PROGRAM BWCOMP 


ENGLISH(@) OR METRIC(1) UNITS? 


=0 

SIG. WAVE HEIGHT , FT(F4.1)? 
=15. 

PEAK WAVE PERIOD, SEC(F4.1)? 
=12. 


DO YOU WANT TO SPECIFY A SEPARATE WAVE HEIGHT 
AND PERIOD FOR WAVE TRANSMISSION COMPUTATIONS(Y OR N)? 


=y 

H SIG. (FOR TRANSMISSION) IN FT =? 
=12. 

PEAK PERIOD(FOR TRANSMISSION) IN SECONDS =? 
=10. 


MAX. ALLOWABLE TRANSMITTED WAVE HEIGHT, FT(F4.1)? 
=2. 


ZEXCEEDANCE OF TRANSMITTED HEIGHT IN SEA STATE? 


DO YOU WANT TO ALLOW FOR INCREASED RUNUP AND OVERTOPPING 
CAUSED BY ONSHORE WIND(Y OR N)? 
=Y 


ONSHORE WIND VELOCITY, MPH (F5S.1)? 
=35, 

WATER DEPTH AT TOE OF STRUCTURE, FT(F4.1)? 
=19. 

FRESH WATER(@) OR SALT WATER(1)? 
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UNIT WEIGHT OF ROCK IS ASSUMED TO BE 
165 LBS/CUFT(2643 KG/CUM). DO YOU WISH TO ENTER AN 
ALTERNATE VALUE? (Y OR N) 


UNIT WEIGHT OF CONCRETE IS ASSUMED TO BE 
149.5 LBS/CUFT(2423 KG/CUM). DO YOU WISH TO ENTER AN 


ALTERNATE VALUE? (Y OR N) 
=N 


INPUT UNIT PRICE OF UNIFORM ARMOR ROCK,IN PLACE ($/TONS) 
=50. 


INPUT UNIT PRICE OF GRADED ARMOR ROCK,IN PLACE($/TONS) 
345. 


INPUT UNIT PRICES($/CY, IN PLACE) FOR THESE ARMOR UNITS: 
PLAIN CUBE,MOD.CUBE,TETRAPOD,QUADRIPOD,HEXAPOD,TRIBAR,TOSKANE,& DOLOS 


FORMAT(8F7.2);SEPARATE PRICES BY COMMAS 
=80.,95.,100.,100.,100.,100.,105.,115. 


INPUT UNIT PRICE FOR CORE MATERIAL(GRADED ROCK) ,$/TONS 
=35. 
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Sample Output 


11. Below is a sample output for computer program BWCOMP. 


# EXAMPLE EXECUTION OF PROGRAM BWCOMP # 
TABLE OF COMPARATIVE BREAKWATER QUANTITIES AND COSTS 


SEAWARD SLOPE = 1 : 1.5 LEEWARD SLOPE = 1: 1.5 WATER DEPTH = 19 FT 
INCIDENT WAVE : SIG HEIGHT = 15 FT PEAK PERIOD = 12 SEC 

(DEPTH LIMITED WAVE HEIGHT = 14.3331 FT) 

INCIDENT WAVE (FOR TRANSMISSION COMPUTATIONS): 

$1G. HEIGHT = 12 FT PEAK PERIOD = 1@ SEC 

MAX TRANSMITTED WAVE HEIGHT = 2 FT({ 1 % EXCEEDANCE) 


RELATIVE VALUES PER UNIT TRUNK LENGTH 


ARMOR UNIT SIZE CREST CREST NO.ARMOR ARMOR CORE ARMOR CORE 


HEIGHT WIDTH UNITS VOL VOL COST COST 
TONS Fan FT CY CY ($) ($) 
ROCK, UNIF 20.6 25.0 18.5 4.9 71.4 66.4 9808 R21 
ROCK, GRAD 18.7 29.0 18.8 5.64 75.3 84.6 4756 154 
CUBE 17.6 25.0 20.9 3.3 46.6 SI aay 3726 109 
MOD. CUBE aS 20.0 16.3 6.2 28.9 88.3 2749 146 
TETRAPOD 8.8 18.8 15.08 5.4 235.6 49.6 2361 98 
QUADRIPOD 8.8 19.0 13.7 Del 22.4 59.8 2240 101 
HEXAPOD Tend 19,8 15.9 Unt 27.4 91.2 2742 93 
TRIBAR 6.8 28.8 13.6 6.8 28.5 68.7 2050 110 
TOSKANE 5.6 21.0 12.8 7.4 20.6 65.7 2168 120 
DOLOS 4.1 21.8 10.6 Tas 15.3 69.9 1763 127 
NOTES: CREST HEIGHT IS ABOVE STILL WATER LEVEL 


1. 
2. ARMOR VOLUME FOR ROCK IS THE TOTAL CROSS-SECTION VOLUME 
3. ARMOR VOLUME FOR ARTIFICIAL UNITS IS THE CONCRETE VOLUME 
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SEAWARD SL 


% * 


* *& 


EXAMPLE EXECUTION OF PROGRAM BWCOMP 


* * 


TABLE OF COMPARATIVE BREAKWATER QUANTITIES AND COSTS 


OPE = 


INCIDENT WAVE : 
(DEPTH LIMITED WAVE HEIGHT = 14.3331 FT) 


INCIDENT WAVE 


SIG. HEIGH 


ARMOR UNIT 


ROCK, UNIF 
ROCK, GRAD 
CUBE 

NOD. CUBE 
TETRAPOD 
QUADRIPOD 
HEXAPOD 
TRIBAR 
TOSKANE 
DOLOS 


NOTES: 1. 


T = 


AGTH ES 


LEEWARD SLOPE = 13; 1.5 WATER DEPTH 
SIG HEIGHT = 15 FT 


(FOR TRANSMISSION COMPUTATIONS): 
12 SE 
MAX TRANSMITTED WAVE HEIGHT = 2 FT( 1 % EXCEEDANCE) 


PEAK PERIOD = 10 SEC 


PEAK PERIOD 


% 


=) 419) Fal 
= 12, SEC 


RELATIVE VALUES PER UNIT TRUNK LENGTH 


SIZE CREST CREST NO.ARMOR ARMOR CORE ARMOR 
HEIGHT WIDTH UNITS VOL VOL COST 

TONS FT FT CY CY ($) 
18.3 15.0 14.8 8.1 59.1 595.8 4147 
9.4 17.0 14.3 9.1 60.6 66.4 3828 

8.8 15.8 16.6 8.8 38.4 49.9 3074 

4.7 13.0 13.8 10.9 25.6 76.2 2428 

4.4 12.0 12.0 9.8 21.4 31.3 2136 

4.4 11.8 10.9 8.8 19,1 49.9 1998 

3.8 12.0 L257 12.6 24.1 49.4 2407 

3.4 15.8 10.8 11.5 19.4 68.8 1942 

2.8 14.0 10.2 13.3 18.5 65.4 1941 

2.1 14.0 8.4 13.6 13.8 69.9 1586 

CREST HEIGHT IS ABOVE STILL WATER LEVEL 

ARMOR VOLUME FOR ROCK IS THE TOTAL CROSS-SECTION VOLUME 
ARMOR VOLUME FOR ARTIFICIAL UNITS IS THE CONCRETE VOLUME 
* * #* #* EXAMPLE EXECUTION OF PROGRAM BWCOMP *  # 


TABLE OF COMPARATIVE BREAKWATER QUANTITIES AND COSTS 


SEAWARD SLOPE = 


INCIDENT W 


AVE : 


Nr ees 


LEEWARD SLOPE = 1: 2 WATER DEPTH 
SIG KETGHT = lelS.FT 


(DEPTH LIMITED WAVE HEIGHT = 14.3331 FT) 
AVE (FOR TRANSMISSION COMPUTATIONS): 
SIG. HEIGHT = 


INCIDENT W 


ARMOR UNIT 


ROCK, UNIF 
ROCK, GRAD 
CUBE 

MOD. CUBE 
TETRAPOD 
QUADRIPOD 
HEXAPOD 
TRIBAR 
TOSKANE 
DOLOS 


NOTES: 


12 FT 
MAX TRANSMITTED WAVE HEIGHT = 2 FT( 1 % EXCEEDANCE) 


PEAK PERIOD = 1@ SEC 


PEAK PERIOD 


=. 12 


LO aRiT 
SEC 


* # 


RELATIVE VALUES PER UNIT TRUNK LENGTH 


SIZE CREST CREST NO.ARMOR ARMOR CORE ARMOR 
HEIGHT WIDTH UNITS VOL VOL COST 

TONS ais al CY CY ($) 
15.9 20.8 16.9 6.2 68.6 68.4 4817 
14.0 24.0 16.3 TES 73.5 89.5 4643 
13.2 20.0 19.8 6.8 44,7 at 7 3576 
ese 17.0 14.8 8.2 29.08 88.0 2754 

6.6 15.0 13.7 7.2 23.9 56.08 2353 

6.6 16.0 W2ata 6.9 22.4 63.1 2244 

9.8 16.0 14.5 9.6 27.4 58.1 2739 

Del 18.0 12.4 8.4 2a 73.4 2123 

4.2 18.@ 11.7 10.2 20.8 74.8 2189 

Send 18.0 9.6 10.2 15.9 79.95 1786 

CREST HEIGHT IS ABOVE STILL WATER LEVEL 

ARMOR VOLUME FOR ROCK IS THE TOTAL CROSS-SECTION VOLUME 
ARMOR VOLUME FOR ARTIFICIAL UNITS IS THE CONCRETE VOLUME 
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# # EXAMPLE EXECUTION OF PROGRAM BWCOMP * # 
TABLE OF COMPARATIVE BREAKWATER QUANTITIES AND COSTS 


SEAWARD SLOPE = 1 : 2 LEEWARD SLOPE = 1: 1.5 WATER DEPTH = 19 FT 
INCIDENT WAVE : SIG HEIGHT = 15 FT PEAK PERIOD = 12 SEC 
(DEPTH LIMITED WAVE HEIGHT = 14.3331 FT) 

INCIDENT WAVE (FOR TRANSMISSION COMPUTATIONS) : 

SIG. HEIGHT = 12 FT PEAK PERIOD = 18 SEC 

MAX TRANSMITTED WAVE HEIGHT = 2 FT( 1 % EXCEEDANCE) 


RELATIVE VALUES PER UNIT TRUNK LENGTH 


ARMOR UNIT SIZE CREST CREST NO.ARMOR ARMOR CORE ARMOR CORE 


HEIGHT WIDTH UNITS VOL VOL COST COST 

TONS FT Fal CY CY ($) ($) 

ROCK, UNIF 15.5 20.0 16.9 5.8 63.7 99.2 4473 108 
ROCK, GRAD 14.0 24.0 16.3 6.8 68.1 77.8 4299 142 
CUBE 13.2 20.0 19.8 6.4 41.5 93.2 S323 97 
MOD. CUBE ea! 17.0 14.8 Tat 27.0 77.9 2561 142 
TETRAPOD 6.6 15.8 13.7 Gia7. 219 48.5 2190 88 
QUADRIFOD 6.6 16.0 12.9 6.4 20.9 54.8 2086 108 
HEXAPOD 5.8 16.8 14.5 8.9 25.9 58.3 2548 92 
TRIBAR ee 18.0 12.4 7.8 HON. 64.8 1978 116 
TOSKANE 4.2 18.8 WET 9.3 LOS 65.2 2031 119 
DOLOS 3.1 18.0 9.6 9.4 14.4 69.4 1655 126 


NOTES: 1. CREST HEIGHT IS ABOVE STILL WATER LEVEL 
2. ARMOR VOLUME FOR ROCK IS THE TOTAL CROSS-SECTION VOLUME 
3. ARMOR VOLUME FOR ARTIFICIAL UNITS IS THE CONCRETE VOLUME 


* #* # EXAMPLE EXECUTION OF PROGRAM BWCOMP * * # 


TABLE OF COMPARATIVE BREAKWATER QUANTITIES AND COSTS 


SEAWARD SLOPE 1: 2.5 LEEWARD SLOPE = 1: 1.5 WATER DEPTH = 19 FT 
INCIDENT WAVE : SIG HEIGHT = 15 FT PEAK PERIOD = 12 SEC 

(DEPTH LIMITED WAVE HEIGHT = 14,3331 FT) 

INCIDENT WAVE (FOR TRANSMISSION COMPUTATIONS): 

SiG. HEIGHT = %2 Fi PEAK PERIOD = 18 SEC 

MAX TRANSMITTED WAVE HEIGHT = 2 FT( 1 % EXCEEDANCE) 


RELATIVE VALUES FER UNIT TRUNK LENGTH 


ARMOR UNIT SEZE CREST CREST NO.ARMOR ARMOR CORE ARMOR CORE 


HEIGHT WIDTH UNITS VOL VOL COST COST 
TONS Ean Fall CY GN, ($) ($) 
ROCK, UNIF 12.4 17.0 at 6.9 68.4 94.9 4240 103 
ROCK, GRAD 1i.2 20.0 USi02 Unie) 63.4 iliee 4806 129 
CUBE 18.6 17.8 Veal ss 39.3 98.6 3147 92 
MOD. CUBE De7 15.0 13.8 ices 26.2 Wks) 2489 142 
TETRAPOD 5.3 13.0 V2 7 Chae rash, ve 48.6 2122 88 
QUADRIPOD 3.9 13.8 THIEN Tis 1956 o1.3 1961 93 
HEXAPOD 4.6 14.0 13.5 10.8 PaCS 5Q.7 2472 92 
TRIBAR 4.1 16.0 Ap ores etl LSS 64.9 922 118 
TOSKANE 3.4 15.08 10.8 11.0 18.3 61.9 1921 113 
DOLOS Psi) 16.0 Bin? PLES Lara 70.6 1618 128 


NOTES: 1. CREST HEIGHT IS ABOVE STILL WATER LEVEL 
2. ARMOR VOLUME FOR ROCK IS THE TOTAL CROSS-SECTION VOLUME 
3. ARMOR VOLUME FOR ARTIFICIAL UNITS IS THE CONCRETE VOLUME 
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*# # * % % #% EXAMPLE EXECUTION OF PROGRAM BWCOMP * * #* #* ¥# 
TABLE OF COMPARATIVE BREAKWATER QUANTITIES AND COSTS 


SEAWARD SLOPE = 1 : 2.5 LEEWARD SLOPE = 1 : 2 WATER DEPTH = 19 FT 
INCIDENT WAVE : SIG HEIGHT = 15 FT PEAK PERIOD = 12 SEC 
(DEPTH LIMITED WAVE HEIGHT = 14.3331 FT) 

INCIDENT WAVE (FOR TRANSMISSION COMPUTATIONS): 

SIG. HEIGHT = 12 FT PEAK PERIOD = 1@ SEC 

MAX TRANSMITTED WAVE HEIGHT = 2 FT( 1 % EXCEEDANCE) 


RELATIVE VALUES PER UNIT TRUNK LENGTH 


ARMOR UNIT SIZE CREST CREST NO.ARMOR ARMOR CORE ARMOR CORE 


HEIGHT WIDTH UNITS VOL VOL COST COST 
TONS FT FT Cy CY. ($) ($) 
ROCK, UNIF 12.4 17.8 Need 7.3 64.5 64.4 4324 117 
ROCK, GRAD Vt.:2 20.0 Nei 8.5 67.8 80.8 4284 147 
CUBE 18.6 17.0 Lia7 8.0 41.9 58.8 3355 105 
MOD. CUBE 5.7 15.0 {S58 ios 27.9 84.9 2653 158 
TETRAPOD 343 13.0 12.7 8.6 220 55.4 2259 181 
QUADRIPOD 5.3 13.0 1G 8.0 20.9 58.3 2089 106 
HEXAPOD 4.6 14.8 13%05 11.5 26.3 iol 2634 105 
TRIBAR 4.1 16.0 Oa) 18.1 20.5 T3330 2053 13 
TOSKANE 3.4 15.0 10.8 Die? 19.5 70.1 2051 127 
DOLOS 2.9 16.0 8.9 12.3 15.8 Ub kel 1730 145 


NOTES: 1. CREST HEIGHT IS ABOVE STILL WATER LEVEL 
2. ARMOR VOLUME FOR ROCK IS THE TOTAL CROSS-SECTION VOLUME 
3. ARMOR VOLUME FOR ARTIFICIAL UNITS IS THE CONCRETE VOLUME 


* * * # # *% # EXAMPLE EXECUTION OF PROGRAM BWCOMP * * # #8 # 
TABLE OF COMPARATIVE BREAKWATER QUANTITIES AND COSTS 


SEAWARD SLOPE 1: 3 LEEWARD SLOPE = 1: 2 WATER DEPTH = 19 FT 
INCIDENT WAVE : SIG HEIGHT = 15 FT PEAK PERIOD = 12 SEC 
(DEPTH LIMITED WAVE HEIGHT = 14.3331 FT) 

INCIDENT WAVE (FOR TRANSMISSION COMPUTATIONS): 

SIG. HEIGHT = 12 FT PEAK PERIOD = 18 SEC 

MAX TRANSMITTED WAVE HEIGHT = 2 FT( 1 % EXCEEDANCE) 


RELATIVE VALUES PER UNIT TRUNK LENGTH 


ARMOR UNIT SIZE »sCREST CREST NO.ARMOR ARMOR CORE ARMOR CORE 


HEIGHT WIDTH UNITS VOL VOL COST COST 
TONS Pu FT CY Gy ($) ($) 
ROCK, UNIF 10.3 15.0 14.8 8.95 62.6 63.0 4393 1S) 
ROCK, GRAD 9.4 17.0 14.3 9.6 64.3 Teac 4062 136 
CUBE 8.8 15.0 16.6 9.3 40.7 536.6 3254 103 
MOD. CUBE 4.7 13.0 13.d 11.5 27.8 B4.3 2569 153 
TETRAPOD 4.4 12.0 12.8 18.4 22.6 57.7 2259 185 
QUADRIPOD 4.4 11.@ 10.9 a) 28.2 36.1 2016 1082 
HEXAPOD S08 12.0 2.7 13.4 25.4 55.7 2945 101 
TRIBAR 3.4 15.0 10.8 12.2 20.6 77.8 2068 148 
TOSKANE 2.8 14.0 10.2 14.1 19.6 73.3 2057 LSS 
DOLOS Fare 14.80 B.4 14.4 14.6 78.1 1683 142 
NOTES: 1. CREST HEIGHT IS ABOVE STILL WATER LEVEL 


1 
2. ARMOR VOLUME FOR ROCK IS THE TOTAL CROSS-SECTION VOLUME 
3. ARMOR VOLUME FOR ARTIFICIAL UNITS IS THE CONCRETE VOLUME 
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Sample Program Listing 


12. The program listing for BWCOMP (FORTRAN version) is as follows: 


1@C SEH HEA AE SE HE AE AE HE HEE TEE AE TE SE HE GE TEA AE TE THE HE SESE HE SE HEE ATE ATER HEHE HE HE 

20C % PROGRAM BWCOMP - VERSION 9/85 % 

30C % CONTACT FOR QUESTIONS ON PROGRAM USE: * 

40C * ORSON SMITH * 

5@C % COASTAL DESIGN BRANCH * 

60C * COASTAL ENGINEERING RESEARCH CENTER * 

70C * U. S. ARMY ENGINEER WATERWAYS EXPERIMENT ¥# 

Bac * STATION, P. O. BOX 431 * 

98C # VICKSBURG, MS 39180 * 

10@C 3 He aE AE SE AE SE ETE AE EAE TE JE EE HE JE EE IE HE aE SE GE ae ae SE HE EE OE JER ab OE SE aE ea dE RE SE aE aE a a HE 

118 CHARACTER ANS#2,ANSW#2,WINDU#3, ANSWIND#2 

120 CHARACTER ANW#2,LU#2,RHOU*B8, VOLU#2,UNIT#10,WU%4 

130 CHARACTER TITLE*8O 

148 REAL KT,NO,KD,LO 

158 DIMENSION XKDELTA(1G) ,H(2) ,T(2) ,RA(18) ,RB(10) 

168 DIMENSION XKDBR(10) ,XKDNBR(1@) ,P(10) ,UNITC(11) 

1780C * DATA 1S FOR ARMOR UNITS: (1) UNIFORM QUARRYSTONE, (2) GRADED * 
180C # RIPRAP, (3) PLAIN CUBES, (4) MODIFIED CUBES, (5) TETRAPODS % 
198C * (6) QUADRIPODS, (7) HEXAPODS, (8) TRIBARS, (9) TOSKANES, * 
208C * (18) DOLOSSE - ALL FOR RANDOM PLACEMENT ON TRUNK + 
210C x 

220C * * STABILITY COEFFICIENTS - BR: BREAKING; NBR: NON-BREAKING * *# 
2308 DATA XKDBR/2.0,2.2,3.5,6.5,7.0,7.0,8.0,9.0,11.0,15.8/ 

248 DATA XKDNBR/4.0,2.5,4.0,7.5,8.8,8.8,9.5,10.0,22.0,31.8/ 

250C * * LAYER COEFFICIENTS * * 

260 DATA XKDELTA/1.@,1.0,1.15,1.1,1.04,.95,1.15,1.82,1.03,.94/ 

270C * * POROSITIES # # 

280 DATA P/. 4423752434247 526546494147 ,.54,002,.06/ 

298C # * RUNUP COEFFICIENTS A & B(AHRENS & McCARTNEY, 1975) * * 

308 DATA RA/.775,.956,.775,.95,1.01,.59,.92,1.81,.986,.988/ 

318 DAMA RE SGL, CSI Ba SO OF ga lig Oye Og lie il yrs OS LO S/ 

328C HEKKKKEAHKEKHEERKEHRHAHHE RHEE EEH EERE ERE REHHEREKE 

330C 

34QC ***##BEGIN INTERACTIVE INPUT OF DATA##** 

350C 

360 WRITE (6,18) 

378 10 FORMAT(1H1,///,5X,"BWCOMP IS AN INTERACTIVE PROGRAM WHICH 

388 & COMPUTES BREAKWATER ",/,"VOLUMES AND COSTS FOR A SIMPLE 

3968 & PARAMETERIZED CROSS SECTION FOR ",/,"THE PURPOSE OF COM 

400 &PARING THE EFFECT OF CHANGING THE SEAWARD ",/,"AND LEEWARD 

418 & SLOPES ON THE SIZE AND RELATIVE COST OF A RANGE OF 

4208 & ",/,"ARMOR UNIT TYPES. THE OUTPUT OF BWCOMP SHOULD NOT 

430 & BE USED ",/,"AS A COST ESTIMATING TOOL IN ANY STAGE OF A 

448 & PLANNING OR DESIGN ",/,"PROJECT. THE ASSUMPTIONS APPLIED 

458 & ARE INTENDED TO TENTATIVELY ",/,"IDENTIFY THE OPTIMUM 

468 & COMBINATION OF SLOPES AND ARMOR UNIT ",/,"FOR A GIVEN 


478 &SET OF DESIGN CRITERIA.") 

480 15 WRITE(6,20) 

490 28 FORMAT(/,5X,"TITLE OF THIS RUN? (UP TO 78 CHARACTERS} ") 
508 READ(S,21) TITLE 
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518 21 FORMAT (ABQ) 


520 WRITE (6,300) 
538 300 FORMAT(/,5X,"ENGLISH(®) OR METRIC(1) UNITS?") 
348 READ(S, 301) IUNITS 


55@ 381 FORMAT(12) 
568 23 FORMAT (AZ) 
578 24 FORMAT (12) 
58@ 25 FORMAT(V) 


598 TF (IUNITS) 302,382,303 
600 302 LU="FT" 
618 60 TO 304 


628 303 LU=" M" 
638 304 WRITE(6,35) LU 
640 35 FORMAT(/,3X,"SIG. WAVE HEIGHT , ",A2,"(F4.1)?") 


658 READ(S,25)H(1) 

668 WRITE (6,45) 

678 45 FORMAT(3X,"PEAK WAVE PERIOD, SEC(F4.1)?") 

688 READ(5,25)T(1) 

690 WRITE (6,46) 

708 46 FORMAT(3X,"DO YOU WANT TO SPECIFY A SEPARATE WAVE HEIGHT", 
710 &/,"AND PERIOD FOR WAVE TRANSMISSION COMPUTATIONS(Y OR N)?") 
728 READ(S,48) ANSW 

738 48 FORMAT (AL) 

740 IF(ANSW.EQ."N") GO TO 119 

730 WRITE(6,49) LU 


768 49 FORMAT(3X,"H SIG. (FOR TRANSMISSION) IN “,A2," 2?") 

770 READ(S,51) H(2) 

788 Si FORMAT (F4.1) 

798 WRITE (4,52) 

BOO 52 FORMAT(3X,"PEAK PERIOD(FOR TRANSMISSION) IN SECONDS =?") 
818 READ(5,53) T(2) 

820 53 FORMAT(F4.1) 


838 GO TO 121 
B40 119 H(2)=H(1) 
850 T(2)=T(1) 


B68 121 WRITE(6,47) LU 
870 47 FORMAT(/,3X,"MAX. ALLOWABLE TRANSMITTED WAVE HEIGHT, ", 
880 &A2,"(F4.1)2") 


890 READ(S,25) HTMAX 

908 WRITE (6,122) 

91@ 122 FORMAT(/,3X,"ZEXCEEDANCE OF TRANSMITTED HEIGHT IN SEA STATE?") 
920 READ(S,91) EXC 

930 EXCP=EXC 

940 EXC=EXC/100. 

958 WRITE (6,123) 

960 123 FORMAT(/,3X,"DO YOU WANT TO ALLOW FOR INCREASED RUNUP 

970 & AND OVERTOPPING",/,3X,"CAUSED BY ONSHORE WIND(Y OR N)?") 
980 READ(S,48) ANSWIND 

998 IF(ANSWIND.EQ@."N") GO TO 124 

1000 IF(ITUNITS.EQ.0) WINDU="MPH" 

1010 IF(IUNITS.EQ@.1) WINDU="M/S" 

1828 WRITE (6,125) WINDU 


1038 125 FORMAT(/,3X,"ONSHORE WIND VELOCITY, ",A3," (F5.1)?") 
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1048 
1858 
1060 
10878 
1088 
1098 
1188 
1110 
1120C 
1138C 
1140C 
1158 
1160 
1178 
1188 
1198 
1200 
1210 
1220 
1238 
1248 
1250 
1268 
1278 
1280 
1298 
1308 
1318 
1328 
1338 
1348 
1350 
1368 
1378 
1380 
1398 
1488 
1418 
1428 
1430 
1448 
1458 
1468 
1478 
1480 
1490 
1508 
1518 
1528 
1538 
1548 
1558 
1560 


READ(5,126) WIND 
126 FORMAT(FS.1) 
124 WRITE(6,120) LU 
120 FORMAT(3X,"WATER DEPTH AT TOE OF STRUCTURE, ",A2,"(F4.1)?") 
READ (5,25) D 
WRITE (6,130) 
130 FORMAT(3X,"FRESH WATER(@) OR SALT WATER(1)?") 
READ(5,24) WATER 


* CHECK FOR BREAKING BY GODA’S FORMULA FOR HORIZONTAL, BOTTOM 


DO 6@ I=1,2 
HI=H(1) 
IF(IUNITS.EQ.@) LO=5.12#T(1) *T (1) 
IF(IUNITS.EQ.1) LO=1.56#T(1)#T(1) 
HB=8.17*LO#(1.-EXP(-4.71239*D/LO)) 
IF(HB.LT.H(I)) HI=HB 
IF(I.EQ.1) HI=HI 
IF(1.EQ.2) H2=HI 

68 CONTINUE 
H2P=H2*SQRT ( (-ALOG(EXC)) /2.) 
WRITE (6,135) 

135 FORMAT(3X,"UNIT WEIGHT OF ROCK IS ASSUMED TO BE",/, 
&3X,"165 LBS/CUFT(2643 KG/CUM). DO YOU WISH TO ENTER AN ",/, 
&3X, "ALTERNATE VALUE? (Y OR N) ") 

READ (5,23) ANS 

IF(ANS.EQ."¥") GO TO 55 

IF(ANS.EQ."N") GO TO 77 
55 IF(IUNITS) 305,305,306 
305 RHOU="LBS/CUFT" 

GO TO 307 

306 RHOU=" KG/CUM " 

307 WRITE(6,145) RHOU 

145 FORMAT(3X,"ENTER UNIT WEIGHT OF ROCK, ",A8,"(F6.1)") 

READ (5,25) RHOR 
GO TO 88 
77 IF(IUNITS) 78,78,79 
78 RHOR=165. 
GO TO 88 
79 RHOR=2643. 
88 CONTINUE 
WRITE (6,335) 

335 FORMAT(3X,"UNIT WEIGHT OF CONCRETE IS ASSUMED TO BE",/, 
&3X,"149.5 LBS/CUFT(2423 KG/CUM). DO YOU WISH TO ENTER AN ",/, 
&3X, "ALTERNATE VALUE? (Y OR N) ") 

READ(5,23) ANS 
IF(ANS.EQ."Y") GO TO 555 
IF(ANS.EQ."N") GO TO 777 
555 WRITE(6,345) RHOU 
345 FORMAT(3X,"ENTER UNIT WEIGHT OF CONCRETE, ",AB,"(F6.1)") 
READ(5,25) RHOC 
GO TO 888 
777 IF(IUNITS) 778,778,779 
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1570 778 RHOC=149.5 
1580 60 TO 888 

1590 779 RHOC=2423. 

1608 888 CONTINUE 

1610 IF(IUNITS) 88,80,93 
1620 8@ IF(WATER) 81,81,82 
1630 81 RHOW=62.4 


1640 GO TO 86 
1658 82 RHOW=464. 
1660 GO TO 86 


1670 83 IF(WATER) 84,84,85 
1680 84 RHOW=1000. 

1698 GO TO 86 

1700 85 RHOW=1025.6 

1718 86 CONTINUE 


1728 IF(TUNITS.EQ@.0) VOLU="CY" 
1730 IF(IUNITS.EQ.1) VOLU="CM" 
1748 IF(IUNITS.EQ@.8) WU="TONS" 
1758 IF(IUNITS.EQ.1) WUs" MT " 
1768 WRITE(6,700) WU 


1778 708 FORMAT(/,3X,"INPUT UNIT PRICE OF UNIFORM ARMOR ROCK, 
1788 &IN PLACE($/",A4,")") 


1798 READ(S,71@) UNITC(1) 

1888 710 FORMAT(F7.2) 

1818 WRITE (6,728) WU 

1828 728 FORMAT(/,3X,"INPUT UNIT PRICE OF GRADED ARMOR ROCK, 
1838 &IN PLACE ($/",A4,")") 

18480 READ(5,718) UNITC(2) 

1850 WRITE (6,808) VOLU 

1868 880 FORMAT(/,3X,"INPUT UNIT PRICES($/",A2,", IN PLACE) 
1878 & FOR THESE ARMOR UNITS:") 

1880 WRITE (6,810) 

1898 8108 FORMAT(/,3X,"PLAIN CUBE,MOD.CUBE, 

1908 &TETRAPOD,QUADRIPOD,HEXAPOD,TRIBAR,TOSKANE,& DOLOS") 
1918 WRITE (6,820) 

1928 9828 FORMAT(/,3X,"FORMAT(8F7.2);SEPARATE PRICES BY COMMAS") 
1938 READ(S5,830) (UNITC(I) ,1=3,10) 

1948 838 FORMAT (8F7.2) 

1950 WRITE (6,848) WU 

1968 840 FORMAT(/,3X,"INPUT UNIT PRICE FOR CORE MATERIAL 
1970 & (GRADED ROCK) ,$/",A4) 

1980 READ(5,850) UNITC(11) 

1998 850 FORMAT (F7.2) 

2000C 

2018C * * SET COT LEEWARD SLOPE * * 

2028C 

2030 DO 999 I1=1,2 

2040 IF(II.€Q@.1) COTP=1.5 

2050 IF(IT.EQ.2) COTP=2.8 

2068C 

2078C * * SET COT SEAWARD SLOPE * 

2088C 

2090 DO 998 JJ=1,4 
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2100 IF(II.E@.2 .AND. JJ.EQ.1) GO TO 998 


2110 IF(JJ.EQ.1) COT#1.5 
2120 IF(JJ.EQ@.2) COT=2.8 
2130 IF(JJ.E@.3) COT=2.5 
2148 IF(JJ.E@.4) COT=3.0 
2158 THETA=ATAN(1./COT) 
2168 THETAP=ATAN(1./COTP) 
21706C 


2180C * * PRINT HEADINGS FOR OUTPUT TABLE * * 
2190C 


2200 WRITE(6,851) TITLE 

2210 851 FORMAT(1H1,11X,"*  ",ABO," #") 

2220 WRITE (6,920) 

2230 920 FORMAT(//,30X,"TABLE OF COMPARATIVE BREAKWATER 

2248 & QUANTITIES AND COSTS") 

22508 WRITE(6,930) COT,COTP 

2268 930 FORMAT(/,32X,"SEAWARD SLOPE = 1:",F3.1, 

2270 &” LEEWARD SLOPE = 1:",F3.1) 

2280 WRITE(6,931) D,LU 

2298 931 FORMAT(32X,"WATER DEPTH = ",F4.1,1X,A2) 

2300 WRITE(6,940) H(1),LU,T(1) 

2310 940 FORMAT(32X,"INCIDENT WAVE: SIG. HEIGHT = ",F4.1,1X,A2, 
2320 &" PEAK PERIOD = ",F4.1,"SEC") 

2330 IF(H1.LT.H(1)) WRITE(6,945) H1,LU 

2340 945 FORMAT(32X,"(DEPTH LIMITED WAVE HEIGHT = ",F4.1,1X,A2,")") 
2350 IF(ANSW.EQ."Y") WRITE(6,951) H(2),LU,T(2) 

2360 951 FORMAT(32X,"INCIDENT WAVE(FOR TRANSMISSION COMPUTATIONS):", 
2370 &/,32X%,"S1G. HEIGHT = ",F4.1,1X,A2," PEAK PERIOD = ",F4.1, 
2380 &"SEC") 

2390 IF (ANSW.EQ."Y".AND.H2.LT.H(2)) WRITE (6,945) H2,LU 

2400 WRITE(6,950) HTMAX,LU 

2410 950 FORMAT(32X,"MAX TRANSMITTED WAVE HEIGHT = ",F4.1,1X,A2) 
2420 WRITE(6,952) EXCP 

2430 952 FORMAT(32X,"(",F4.1,"% EXCEEDANCE)") 

2440 WRITE (6,960) 

2450 960 FORMAT(/,53X,"* * #* RELATIVE VALUES PER UNIT 

2460 & TRUNK LENGTH # * #") 

2478 WRITE (6,972) 

2480 970 FORMAT(/,11X,"ARMOR UNIT",2X," SIZE ",2X," CREST ", 
2490 &2X," CREST ",2X,"NO.ARMOR",2X," ARMOR ",2X," CORE ", 
2500 &2X," ARMOR ",2X," CORE ",2X," TOTAL ") 

2510 WRITE(6,980) WU,LU,LU,VOLU, VOLU 

2520 980 FORMAT(23X," (",A4,") ",2X," HT. (",A2,")",2X,"WDTH(",A2Z,")", 2X, 
2530 &" UNITS ",3X,"VOL(",A2,")", 3%, "VOL6",A2,")", 2X," COST(#)", 
2540 &2X," COST($)",2X," COST($)") 

2550C 

2560C * «* SET ARMOR UNIT TYPE * «& 

2570C 

2580 DO 997 NN=1,12 

2598 ICHECK=0 

2600 RHO=RHOC 

2610 IF(NN.EQ.1) RHO=RHOR 

2620 IF(NN.EQ.2) RHO=RHOR 
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2630 IF(ABS(H1-H(1)) «LT. .888000001) KD=XKDNBR (NN) 


2648 IF(H1.LT.H(1)) KD=XKDBR(NN) 

2658 W=RHO#(H1#*3) / (KD¥COT#( ( (RHO/RHOW) ~1.) #3) ) 
2668 AT=2.*((W/RHO) ##0. 33) 

2670 IF(NN.EQ.3) AT22.3#((W/RHO) #*O. 33) 

2688 IF(NN.EQ.4) AT=1.03%((W/(Q.78#RHO) ) ¥#O. 33) 
2690 IF(NN.EQ@.5) AT=1.361#((W/(@.28*RHO) ) *#O. 33) 
2700 IF(NN.EQ.6) AT=1.502%((W/ (8. 495#RHO) ) #*O. 33) 
2710 IF(NN.EQ.7) AT=1.29%((W/(0.176*RHO) ) ##O. 53) 
2720 IF(NN.EQ.8) AT=3.68%((W/(6.4B8*RHO) ) #*O.33) 
2730 IF(NN.EQ.9) AT=0.889#( (W/ (0. @83*RHO) ) **O. 33) 
2748 IF(NN.EQ@.1@) AT=1.02%((W/ (0. 16*RHO) ) ##O. 33) 
2758 TOPW=3.Q*#XKDELTA(NN) *( (W/RHO) ##8. 33) 

2768 IF(IUNITS.EQ@.8) C=1. 

2778 IF(IUNITS.EQ@.1) C=.3 

2780C 

2790C ¥ * INCREASE CREST ELEVATION (C) TO SATISFY HTMAX * 
2800C 

2818 996 IF(ICHECK.EQ@.8) GO TO 588 

2820 IF(IUNITS.EQ.@) C=C+1. 

2830 IF(IUNITS.EQ.1) C=C+.3 

2848 S@@ HS=C+D 

2858 TANB=1./COT 

2868 IF(IUNITS.EQ@.8) SURF=TANB/SQRT(H2P/(5.12#T(2)#T(2))) 
2870 IF(IUNITS.EQ.1) SURF=TANB/SQRT(H2P/(1.56*T(2)*#T(2))) 
2880 55@ RH=RA(NN) *SURF/(1.+RB(NN) #SURF) 

2890 R=H2P#RH 

2908 TF(ANSWIND.EQ@."N") GO TO 551 

2918 IF(ITUNITS.EQ.1) WIND=WIND/. 447 

2920 WF=(WIND##2)/1800. 

2938 IF(WIND.GT.60.) WF=2. 

2940 WK=1.+(WF#(C/R+.1) *SIN( THETA) ) 

2950 R=R*¥WK 

2968 S51 FR=C/K 

2978 CR=0.51 -O.11#TOPW/HS 

2980 KT=CR#(1.-FR) 

2990 HT=H2P*#KT 

3008C 

SO1@C * * CHECK IF HTMAX CRITERIA MET * * 
3820C 

3830 604 IF(HT.LE.HTMAX) GO TO 605 

3040 ICHECK=1 

3050 GO T0 996 

3060C 

3878C * * VOLUME AND COST COMPUTATIONS * * 
3080C 

3090 685 CONTINUE 

3100 TVOL=(TOPW#HS) +( (HS#*#2) #COT/2.)+( (HS**#2) #(COTP/2.)) 
3118 AH=C+(1.5#H1) 

3120 IF (AH. GT.HS) AH=HS 

3138 SAH=C+(@.5#H1) 

3140 AVOL=(TOPWHAT) +(AT*(AH/SIN(THETA) )) +(AT#(SAH/SIN(THETAP) )) 
3158 CVOL=TVOL-AVOL 
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988 


Wt 


982 


598 
999 


985 


986 


14 


AA=TOPW+(AH/SIN( THETA) )+(SAH/SIN(THETAP) ) 
IF(TUNITS.EQ.8) CVOL=CVOL/27. 

P(1)=0.37 

NO=2. #AAXXKDELTA(NN) #(1.-P (NN) )#((RHO/W) ##8. 67) 
IF(IUNITS.E@.@) CC=(UNITC(11)#. 63#CVOL#RHOR) /2000. 
IF(IUNITS.EQ@.1) CO=(UNITC(11)*. 63*CVOL#RHOR) /1000. 
IF (NN.GE.3) AVOL=(W/RHOC) #NO 

IF (JUNITS.EQ.@) AVOL=AVOL/27. 

IF(IUNITS.E@.8) W=W/2008. 

IF(IUNITS.EQ.1) W=W/1000. 

AC=AVOL*UNITC(NN) 

IF(NN.LE.2) AC=NO#W*HUNITC (NN) 

TC=AC+CC 

IF(NN-EQ.1) UNIT="ROCK,UNIF." 

IF(NN.E@.2) UNIT="ROCK,GRAD.” 

IF(NN.EQ.3) UNIT="CUBE sy 

IF(NN.EQ@.4) UNIT="MOD. CUBE " 

IF(NN.EQ@.5) UNIT="TETRAPOD " 

IF(NN.EQ.6) UNIT="QUADRIPOD " 

IF(NN.EQ.7) UNIT="HEXAPOD a 

IF(NN.EQ.8) UNIT="TRIBAR * 

IF(NN.EQ.9) UNIT="TOSKANE z 

IF(NN.EQ@.10) UNIT="DOLOS ; 


* TABLE OUTPUT * # 


WRITE (6,900) UNIT,W,C,TOPW,NO,AVOL,CVOL,AC,CC,TC 
FORMAT(/,11X,A10,F8.1,2X,8(2X,F8.1)) 

CONTINUE 

WRITE (6,982) 


FORMAT(//,11X,"NOTES: 14. CREST HEIGHT IS ABOVE STILL 
& WATER LEVEL",/,19X,"2. ARMOR VOLUME FOR ROCK IS THE 
& TOTAL CROSS-SECTION VOLUME",/,19X,"3. ARMOR VOLUME 


& FOR ARTIFICIAL UNITS IS THE CONCRETE VOLUME") 
CONTINUE 

CONTINUE 

WRITE (6,985) 

FORMAT(//,3X,"D0 YOU WISH TO MAKE ANOTHER RUN?") 
READ (5,986) ANW 

FORMAT (A2) 

IF(ANW.EQ."Y")GO TO 15 

IF(ANW.EQ."N")GO TO 14 

STOP 

END 
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APPENDIX C: COMPUTER PROGRAM BWLOSS2 


Estimation of Economic Losses from Transmitted Waves 


Program purpose 
1. BWLOSS2 fits an Extremal Type I long-term probability distribution 


to transmitted wave height data and estimates expected annual economic losses 
due to wave attack after a protective breakwater has been built. 
Program capabilities 

2. Estimation of incremental economic benefits directly related to a 
rubble-mound breakwater built for wave protection requires that the costs to 
the beneficiaries with the breakwater in place be determined. A rubble-mound 
breakwater built so high that no waves are transmitted during the worst con- 
ceivable conditions is seldom affordable. It is often more cost effective to 
accept a small amount of risk that waves from a very severe storm will be 
transmitted and cause an estimable degree of economic loss. The damage caused 
by these transmitted waves is assumed to follow a previously derived economic 


loss function of wave height. 


$L(H_) = $L I z gs te | C1) 
s max 
where 
$L(H,) = the economic losses caused by a storm of significant wave 
height He 
Pmax = the maximum conceivable loss due to wave attack 
Ho = the maximum wave height for which losses can be neglected 


A = a coefficient determined by regression of historical H, , 
$L(H,) information 


The computer program BWLOSS1 (Appendix A) is available to derive this economic 
loss function from property valuations, coastal engineering data on wave cli- 


mate, and historical economic loss data. The program BWLOSS2 requires that 


$Linax » HLo , and A be input, along with the e« and 6 parameters of the 
Extremal Type I cumulative probability distribution of significant wave 
heights H, which defines the wave climate incident on the seaward side of 
the breakwater as follows: 


C1 


[(<-4,)/¢] 
F(H,) = ale (C2) 


3. At least two data points of transmitted wave height versus a return 
period of the associated incident wave height are also required to transform 
the incident probability distribution to an Extremal Type I cumulative proba- 
bility distribution of transmitted waves F(H,) . The transformation is ac- 
complished by a least squares fit of the Extremal Type I function above to the 
Ay points, given the cumulative probability of the corresponding incident Ae 
value as represented by the traditional return period. The nonlinear coeffi- 
cient of correlation and sum of least squares are computed to indicate the 
goodness of the fit. A table of residuals is optionally provided. 

4. The transmitted wave heights during any storm represented by He 
are probably not Rayleigh distributed, but the transmitted wave height associ- 
ated with a given-incident significant wave height is assumed to be at the 
13.5 percent exceedance level among all transmitted waves (including those of 
zero height). This is the same exceedance level as the significant wave 
height in a Rayleigh distributed sea state. The methods of Andrew and Smith 
(in preparation) can be applied to estimate the transmitted wave heights at 
other exceedance levels. 

5. The program BWLOSS2, in a manner similar to its sister program 
BWLOSS1, computes an expected, or long-term, average annual economic loss due 


to transmitted waves by the following formulation: 


H 
SY dF(H, ) 
$L'\ _ mca Sa" 
<tc zi) see | TTT ah (C3) 
Lo 


where [dF (H, )/dH, ] is the probability density function f(H,) associated 
with F(H,) . % is the Poisson parameter, or average number of extreme 
events per year, as defined to derive E(He) . The assumption of a random 
number of extreme events per year which can be represented by a mean value in- 
dependent of individual H, (or Hy) values, is critical to the above defini- 
tion of E{$L'/yr} . $L(H,) is taken to describe also $L(H,) since the 
transmitted waves are now the incident waves to the facilities and operations 


incurring losses. H, thus represents the intensity of the transmitted wave 
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climate associated with an incident wave climate of intensity H, and return 
period defined as 


RT(H.) = (C4) 


i 
vt - PCH) ] 

6. The upper limit of integration H 0 is taken as the H, ‘ value 
corresponding to a probability of exceedance in any year of 0.0000001. The 
lower limit of integration is the Die value below which losses are assumed 
as zero. The extrapolation of F(H,) to values of H, below that originally 
used as a threshold for data to derive the e« and @ parameters is probably 
conservative, but this question will be the subject of further study. The 
choice of the threshold H, as equal to Hy, would resolve any problems, 
however. The integration between H,, and ae is approximated numerically 
by an application of Simpson's rule with 100 intervals. 

Program input 

7. The program BWLOSS2 is completely interactive in its present form. 
It is written in FORTRAN IV as implemented on the US Army Engineer Waterways 
Experiment Station (WES) Honeywell DPS-8 mainframe system, but it is easily 


adaptable to microcomputer systems. 


Sample Interactive Session 


8. The following sample interactive session demonstrates the required 


user input: 
INPUT $Lmax FOR LOSS CURVE: $L(Hs)=$Lmax#(1-EXP(A(Hs-HLo) )) 


=20.8 

INPUT HLo, THE MAXIMUM WAVE HEIGHT FOR WHICH LOSSES ARE NEGLIGIBLE 
=2.8 

INPUT A, THE REGRESSION COEFFICIENT 

=-, 086 


INPUT AVERAGE NUMBER OF EXTREMAL EVENTS PER YEAR, 
THE POISSON "LAMBDA® PARAMETER 

=4.0 

INPUT THE NUMBER OF TRANSMITTED WAVE HEIGHT, 

RETURN PERIOD DATA POINTS YOU HAVE. 

23 

INPUT THE DATA POINTS-ONE AT A TIME. 

INPUT TRANSHITTED WAVE HEIGHT,COMMA, THEN RETURN PERIOD IN YEARS 
22,0,20 

=2,7,50 

24,0, 100 

PRINT RESIDUAL TABLES (Y/N)? 

=y 
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LEAST SBUARES RESULTS 


EXTREMAL TYPE I 


Fihs)=Pr(Hs<hs)= EXP(-EXP(-(hs-EPSI) /PHI)) 
EPSI= 3,856 


PHI= 1.294 
HEAN= 3.189 
VARIANCE= 2.755 
XVALUE YVALUE YEST DIFF 


2.0008 8.9875 0.9892 0.8017 
2.7008 8.9958 0.9937 0.8013 
4, 0800 8.9975 8.9977 0.8002 


NON-LINEAR CORRELATION 1S 8. 91378593 


SUM S@R RESIDUALS IS 0. 08808 

RETURN PERIOD TABLE 

YEAR Hs 
5.00 8.01 
10.68 8.98 
25.88 2.18 
50. 88 3.08 
180. 88 3.98 

WEIBULL 


Fihs)=Pr(Hs<hs)= — 1-EXP((-(hs/BETA) ##ALPHA) | 
ALPHA= 0.4443 


BETA= @.869 
NEAN= 8.179 
VARIANCE= 8.217 
XVALUE YVALUE YEST DIFF 


2. 0008 Q.9875 8.9886 8.881! 
2. 7008 8.9958 8.9948 8.0018 
4, 0008 8.9975 8.9977 8.0882 


NON-LINEAR CORRELATION IS 8. 95775742 


SUM SAR RESIDUALS iS 8. 88088 

RETURN PERIOD TABLE 
YEAR Hs 
5.08 8.81 

10.88 1.38 

25. 80 2.14 

38. 08 2.93 

108. 88 3.86 


ch 


LOG EXTREMAL 


Fihs)=Pr(Hs<hs)= EXP(-((BETA/hs) ##ALPHA) ) 


ALPHA= 2.3089 
BETA= @. 288 
MEAN= 8.453 
VARIANCE= 8.391 
XVALUE YVALUE YEST DIFF 


2, 0088 8.9875 8. 9885 8.0018 
2. 7008 0.9998 8.9942 0. 0088 
4, 0008 8.9975 0.9977 0.0002 


NON-LINEAR CORRELATION 15 8.97882341 


SUM S@R RESIDUALS IS 8. 80000 

RETURN PERIOD TABLE 
YEAR Hs 

5.88 1.85 

18. 88 1.42 

25, 88 2.13 

50. 08 2.88 

108.80 3.89 


SELECT A DISTRIBUTION... 
EXTREMAL TYPE I..01 
WEIBULL. .... secened 
LOB-EXTREHAL...0003 
SELECT 1, 2, OR 3 

=! 


EXPECTED ANNUAL LOSS IN MILLIONS OF DOLLARS IS 


t 


0. 0865485 
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Sample Program Listing 


9. A sample program listing for BWLOSS2 is given below: 


18C PROSRAM “BWLOSS2". 9/85 VERSION 

20C DESIGN BRANCH-COASTAL ENGINEERING RESEARCH CENTER 
JOC U.S. ARMY ENGINEERS WATERWAY EXPERIMENT STATION 
4OC VICKSBURG, HS 39188 

5QC FOR FUTHER INFORMATION CONCERNING THE APPLICATION 
6@C OF *BWLOSS2", CALL 

70C ORSON P. SMITH (681)-634-2813 FTS:542-2813 OR 
BOC ROBERT B. LUND (681)-634-2868 FTS:542-2868 


50c 

180C INITIALIZE STRINGS, VARIABLES, AND FUNCTIONS 

118 

128 DIMENSION X (208) ,Y¥ (208) ,D4(20@) ,¥2(20@) , 21 (288) 

138 DIMENSION ALPHA(3) ,BETA(3) 

148 REAL LAMBDA 

158 CHARACTER#6@ WORD(14) 

168 CHARACTER#1 L6 

178 WORD (1) =‘ SESGSESEEEEESESEEEESESEELERE ESLER ESE SESE ETRE ESET ESTE SESE 
188 WORD(2)='# "BWLOSS2" IS A PROGRAM WHICH FITS LONG-TERM CURULA- # 
198 WORD(3)='# TIVE PROBABILITY DISTRIBUTIONS TO TRANSMITTED SI6- # 
208 WORD(4)='# NIFICANT WAVE HEIGHT DATA TO ESTIMATE EXPECTED $ 
218 WORD(5)=' ANNUAL ECONOMIC LOSSES DUE TO WAVE ATTACK AFTER A a’ 
228 WORD(6)='# BREAKWATER HAS BEEN BUILT. THE PROGRAM REQUIRES &' 
238 WORD(7)='# DATA ON A PREVIOUSLY DEFINED LOSS .VS. SIGNIFICANT #° 
240 WORD(8)='# WAVE HEIGHT FUNCTION (REF. PROGRAM BWLOSSi) AND AT #° 
258 WORD(9)='# LEAST 2 POINTS OF TRANSMITTED WAVE HEIGHT, RETURN #' 
268 WORD(18)='# PERIOD DATA. THE DISTRIBUTIONS OF TRANSHITTED WAVE #' 
278 WORD(14)='# HEIGHTS ARE CALCULATED BY THE HETHOD OF LEAST #" 
288 WORD(12)='# SQUARES. A RESIDUAL TABLE CAN BE OPTIONALLY #' 
298 WORD(13)='# PRINTED. #! 
300 WORD (14) =WORD(1) 

318 DO 3 T=1,14 

328 WRITE(6,4) WORD(T) 


338 4 FORMAT (A68) 
348 3 CONTINUE 


350 PRINT, ‘ 
368 PRINT, ‘ 
370 


388C GET THE FACTS 
398 5 PRINT, INPUT $Laax FOR LOSS CURVE: $L(Hs)=$Laax#(1-EXP(A(HS-HLo))) ’ 


408 READ, C1 

410 IFIC] .LE. @) 60 10 5 

428 18 PRINT, ‘INPUT HLo, THE MAXIMUM WAVE HEIGHT FOR WHICH LOSSES ARE NEGLIGIBLE’ 
430 READ, C3 

44a IF(C3 .LE. @) 60 70 18 

450 15 PRINT, INPUT A, THE REGRESSION COEFFICIENT’ 

460 READ,C2 

478 IF{ C2 .6T. 8) 60 TO 15 

488 17 PRINT, “INPUT AVERAGE NUMBER OF EXTREMAL EVENTS PER YEAR, ‘ 
490 PRINT, ‘THE POISSON "LAMBDA" PARAMETER’ 

500 READ, LAMBDA 
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O10 
528 20 
038 
540 
308 
568 
378 
588 
998 
600 
618 
628 
638 38 
648 


IF( LAMBDA .LE. @) 60 TO 17 

PRINT, INPUT THE NUMBER OF TRANSMITTED WAVE HEIGHT, ' 
PRINT, ‘RETURN PERIOD DATA POINTS YOU HAVE. ‘ 

READ ,N2 

IF( N2 .LT. 2) PRINT, YOU NEED AT LEAST TWO POINTS’ 

IF( N2 .LT. 2) 60 10 20 

IF( N2 .6T. 188) PRINT, ‘TOO MANY POINTS-100 1S MAXINUN’ 
IF( N2 .6T. 108) 60 TO 20 

PRINT, ‘INPUT THE DATA POINTS-ONE AT A TIME.’ 

PRINT, ‘INPUT TRANSMITTED WAVE HEIGHT, COMMA, THEN RETURN PERIOD IN YEARS ‘ 
DO 30 1=1,N2 

READ, X(1) ,¥(1) 

V(1)=1,-1.0/ (LAMBDAY (1) ) 


650C PUT TRANSMITTED WAVE HEIGHTS IN ORDER 


668 
678 
688 
698 
700 
718 


CALL ORDER (N2,X,Y) 

CALL WAVTRNS(X,Y,ALPHA, BETA, LAMBDA ,N2) 
CALL EXPCT(C1,C2,C3,ALPHA, BETA, LANBDA) 
STOP 

END 


720C SUBROUTINE TO PUT DATA IN ORDER-LOWEST TO HIGHEST 


730 
740 
750 
768 
778 
780 
798 
800 
818 
828 
830 
840 
850 10 
868 28 
870 
888 
898 
900 
910 
920 
930 
940 
958 
960 
970 
980 
990 
1080 
1018 
1020 
1038 


SUBROUTINE ORDER(N2, X,Y) 
DINENSION ¥ (180) ,¥ (108) 
DO 28 K=2,N2 

J=N2-K#2 

DO 18 I=1,J-1 

IF( X(1) LT. X(1#1)) 60 TO 18 
TL=X(I) 

T2=¥(1) 

K(I)=X(141) 

V«T)=V(I+1) 

X(1#1) 271 

¥(141)212 

CONTINUE 

CONTINUE 

RETURN 

END 


SUBROUTINE WAVTRNS(HS,Y,ALPHA, BETA, LAMBDA ,N) 

DIMENSION YACT(200,3) , YEST(200,3) ,DUML (208) ,DUM2(20@) ,HS (208) ,¥ (208) 
DIMENSION YAVG(3) ,CORR(3) ,A(3) ,B(3) ,ST(3) ,SB(3) ,ALPHA(3) ,BETA(3) , VAR(S) 
DIMENSION RET(S) ,CHS(5,3) 

REAL NEAN(3) 
REAL LANBDA 

FL(X)=EXP(-EXP (-(X-EPS1) /PHI)) 
F2(X)=1.O-EXP((-(X/SIGMA) ##C) ) 

F3(X)=EXP(-((SIGMA2/X) ##U) ) 
CHARACTER#2@ IFLAG(3) 

CHARACTER#17 DEF 

CHARACTER#34 FORM(3) 

CHARACTER#26 TITLE(1) 

CHARACTER#! LOGIC 
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1048 CHARACTER#6@ BOX (15) 


1058 

1860C INITIALIZATION OF STRINGS AND CONSTANTS 
1878 IFLAG(1)="EXTREMAL TYPE 1° 

1088 IFLAG(2)= ‘WEIBULL ‘ 

1098 TFLAG(3)='LOG EXTREMAL‘ 

1108 DEF='F(hs)=Pr (Hs<hs)= ' 

1118 FORM(1)="EXP(-EXP(-(hs-EPSI) /PHI))’ 
1128 FORN(2)='1-EXP((-(hs/BETA) ##ALPHA) ) ’ 
1138 FORN(3)="EXP(-((BETA/hs) ##ALPHA) ) ' 
1148 TITLE(1)=‘LEAST SQUARES RESULTS‘ 
1158 DATA RET /5.0,18.8,25.0,58.0, 182. 0/ 
1168 EULER=.5772156649 

1178 C2=.7796968 

1188 16 PRINT, PRINT RESIDUAL TABLES(Y/N) ?° 
1198 READ(5,17) LOGIC 

1208 17 FORMAT (AL) 

1218 IF(LOGIC .NE. ‘N' .AND. LOGIC .NE. ‘Y') 60 TO 16 
1228 DO 25 [=1,N 

1238 DO 30 K=1,3 

1240 YACT (I,K) =¥ (1) 


1258 38 CONTINUE 

1268 25 CONTINUE 

1278 

12800 INITIALIZE VARIABLES FOR LEAST SQUARES FIT OF THE DISTRIBUTIONS 
1298 5X=8 


1308 5Y=8 

1318 SXx=@ 

1328 SLX=@ 

1338 SLLY=0 

1348 SLXX=@ 

1358 SLLOY=@ 

1368 SXLLY=8 

1378 SLXLLY=8 

1388 TOOBIG=B 

1398 

14@8C CALCULATE SUMS FOR THE LEAST SQUARES METHOD 
1418 DO 40 J=1,N 

1428 SX=5X+HS (J) 

1438 SY=SY+VACT (J, 1) 

1448 SXX=SXX+HS (J) #42 

1458 SLX=SLX+AL06 (HS (J) ) 

1460 SLXX=SLXX+ (ALOGB (HS (J) )) #42 

1470 SLLY=SLLY-ALOG(-ALOG(YACT (J,1))) 

1480 SLLOY=SLLOY+ALO6(-ALOG(1.0-YACT(J,1))) 
1498 SXLLY=SXLLY-HS (J) #ALOG(-ALOG(YACT(J,1))) 
1508 SLXLLY=SLYLLY-ALOG(HS (J) ) #ALOG(-ALOG(YACT(J,1))) 


1518 40 TOOBIG=TOOBIG+ALO6 (HS (J)) #(ALOG(-ALOG(1.@-YACT(J,1)))) 
1528 
1538C CALCULATE SLOPE AND INTERCEPT OF EACH "PLOTTED LINE" 


1548 AUL)=(N#SXLLY-SX#SLLY) / (N#SXX-SX##2) 
1558 A(2)=(N#TOOBIG-SLX#SLLOY) / (N#SLXX-SLX##2) 
1568 A(3)=(N#SLXLLY-SLU#SLLY) / (N#SLXX-SLX##2) 
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1578 B(1L)=(SXX#SLLY-SXLLY#5X) / (NESXX-SX##2) 

1580 B(2)=(SLXX#SLLQY-TOOBIG#SLX) / (N#SLXX-SLY##2) 

1598 B(S)=(SLXX#SLLY-SLXLLY#SLY) / (N#SLXX-SLX##2) 

16@@C CALCULATE PARAMETERS OF EACH DISTRIBUTION FROM SLOPE AND INTERCEPT DATA 
1618 PHI=1.0/A(1) 


1620 EPST=-B(1)/A(1) 

1638 C=A(2) 

1640 SIGHA=EXP(-B(2) /A(2)) 
1658 U=A(3) 

1668 SIGMAZ=EXP(~B(3) /A(3)) 
1678 


168@C ASSIGN ARRAYS ALPHA AND BETA THE PARAMETERS OF EACH DISTRIBUTION 
1698C FOR EASY PRINTOUT OF DATA 


1708 ALPHA(1)=EPSI 

1718 BETA(1) =PHI 

1720 ALPHA(2)=C 

1730 BETA(2)=SI6MA 

1748 ALPHA (3) =U 

1758 BETA(3)=SIGHA2 

1760C CALCULATE APPROXIMATE PROBABILITY AS ESTIMATED BY DISTRIBUTION 
1770 DO 108 J=1,N 

1780 YEST(d,1)=FA(HS(J)) 

1798 YEST(J,2)=F2(HS(J)) 

1800 YEST(d,3)=F3(HS(J)) 

1818 188 CONTINUE 

1828 

1830C CALCULATE AVERAGE PROBABILITY AND CORRELATION COEFFICIENTS 
1848 DO 118 K=1,3 

1850 YAVG(K) =SY/FLOAT(N) 

1868 ST (K)=8 

1878 §B(K)=@ 

1888 118 CONTINUE 

1898 

1988 DO 128 K=1,3 

1918 DO 138 J=1,N 

1928 ST(K)=ST(K) + (YACT (I,K) -YEST (I,K) ) ##2 

1938 SB(K) =SB(K) +(YACT (1 ,K)-YAVG(K) ) ##2 

1948 138 = CONTINUE 

1958 CORR (K)=1.8-ST(K) /SB(K) 

1968 128 = CONTINUE 

1978 

198@C CALCULATE DATA FOR RETURN PERIOD TABLES 

1998 DO 57 J=1,9 

2008 PROB=1.8-1.0/ (LAMBDA#RET (J) ) 

2018 IF(PROB .LE. @) PROB=. 008008! 

2028 CHS(J,1)=-ALOG(-ALOG (PROB) ) #PHI+EPS! 

2030 CHS(J,2)=(-ALOG(1.@-PROB))## (1.8/0) #S1GHA 
2048 CHS (J,3)=SIGMA2/ ((-ALOG(PROB) )##(1.8/U)) 
2058 57 CONTINUE 

2068 

2070C CALCULATE MEAN AND VARIANCE FOR EACH DISTRIBUTION 
2080 MEAN (1) =EPSI+EULER#PHI 

2090 VAR (1)=1,6449341 #PH] ##2 
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21088 PARA=1.0+1.8/C 


2118 CALL GAMMA (PARA, WME) 

2120 HEAN (2) =S1GMA#WNE 

2130 FACI=SIGNAS#2¢WME#E2 

2140 PARA=1.0+2.0/C 

2150 CALL GANNA(PARA,WV2) 

2168 FAC2=SIGNA®#2#HV2 

2178 VAR(2)=FAC2-FACI 

2180 PARA=1.0-1.0/U 

2198 CALL GAMMA(PARA,HPC) 

2280 MEAN (3) =SIGMA2#HPC 

2218 PARA=1.0-2.0/U 

2228 CALL GANMA(PARA,HPD) 

2230 VAR (3) =S1GNA2##2#HPD-NEAN (3) ##2 

2240 

2250C WRITE OUT THE DATA FOR EACH DISTRIBUTION 

2260 WRITE(6,135) TITLE(AI 

2278 135 FORMAT(///, 28% ,A26,///) 

2288 DO 150 K=1,3 

2298 KTEMP=K 

2300 WRITE (6,168) IFLAG(K) ,DEF ,FORM(K) 

2310 168  FORMAT(15X,A3B,//,1%,A17,2%,A34) 

2328 IF( K «EQ. 1) WRITE(6,159) EPSI,PHI 

2330 159 FORMAT(1X,"EPSI="6X,F10.3,/,1X,"PHI=", 74 FLO. 3) 
2348 IF( K .6T. 1) WRITE(6, 161) ALPHA(K) BETA(K) 
2350 161  FORMAT(1X,"ALPHA=",5X,F10.4,/, 1X, "BETA=", 6X,F1@. 3) 
2360 WRITE (6,162) MEAN(K) , VAR(K) 

2378 162  FORMAT(1X,*MEAN=",OX,F10.3,/,1X, “VARTANCE=", 2X ,F1Q, 3) 
2380 IF( LOBIC .EQ. 'N’) GO TO 207 

2390 DO 178 1=1,N 

2400 DUNS (1) =YACT (I,K) 

2410 DUN2 (1) =YEST (I,K) 

2420 HOLD1=CORR(K) 

2438 L2=N 

2440 178 CONTINUE 

2450 CALL RESIDUAL (HS, DUM! ,DUN2, HOLD! ,L2) 


2468 287 WRITE(6, 208) 

2470 208  FORMAT(7X,"RETURN PERIOD TABLE",/,6X,"YEAR",13X,"Hs") 
2480 DO 211 J=1,5 

2498 WRITE(6,212) RET(d) ,CHS(J,K) 

2588 212  FORMAT(1X,F9.2,8X,F9.2) 

2518 241 CONTINUE 

2528 WRITE (6, 165) 

2538 165 = FORMAT(////) 

2540 158 CONTINUE 

2558 383 = RETURN 


2568 END 

2578 

2588C SUBROUTINE TO HELP PRINT OUT DATA 

2598 SUBROUTINE RESIDUAL(X,YACT, YEST,CORR,N) 
2680 DIMENSION X(N), YACT(N) , YEST(N) , DIFF (208) 
2618 SSR=8 

2628 DO 18 [=1,N 
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2630 DIFF (1)=(YACT(I)-YEST (1) ) ##2 

2640 10 SSR=SSR+DIFF (1) 

2600 WRITE (6,15) 

2668 15 FORMAT(//,1X," = XVALUE YVALUE YEST 
2670 DO 25 I=1,N 

2688 WRITE (6,20) X(1),VACT(I) ,YEST(I) SQRT(DIFF(1)) 
2690 20 FORMAT(IX FIL 4,FiL.4,FUL.4,FiL.4,/,) 

2788 25 CONTINUE 

2718 WRITE (6,48) CORR, SSR 

2728 40 FORMAT(/,1X,°NON-LINEAR CORRELATION 15",5X,F10.8,/, 
2730 & 1X,°SUM SAR RESIDUALS 15°, 10X,F10.5,//) 

2748 RETURN 

2730 END 

2768 

2770C SUBROUTINE TO EVALUATE THE GANMA FUNCTION 

2780C PROGRAM ADJUSTS ALPHA TO BE BETWEEN 1.8 AND 2.8 
2798C AND THEN MULTIPLIES BY GF 10 COMPENSATE 


2800 SUBROUTINE GAMMA (ALPHA, AREA) 
2810 DOUBLE PRECISION C(25) , SUN 
2820 6F=1.08 

2838 IF(ALPHA) 1,2,3 

2848 

2858 2 PRINT, TROUBLE IN GAMMA’ 
2848 AREA=1. 

2870 60 TO 200 

2888 


2B90C FOR GANHA OF A NEGATIVE NUMBER 
2988 3 M=INT (ALPHA) 


2910 EPSI=ALPHA-FLOAT (M) 

2920 IF( M .EQ. @) GF=6F/ALPHA 
2930 IF( M .EQ. @) ALPHASALPHA+I.@ 
2940 IF( H .EG. @) 60 TO 100 
2950 IF( MEQ. 1) GF=1.0 

2960 IF( M .E0. 1) 60 TO 108 
2978 DO 10 1=2,M 

2980 10  GF=GF#(FLOAT(I~1) +EPSI) 
2998 ALPHA=1.O+EPSI 

3000 80 TO 100 

3010 


3020C FOR GAMMA OF A POSITIVE NUMBER 
3038 1 M=INT (ALPHA) 


3040 EPS1=ALPHA-FLOAT(N) 

3050 IF( M .EQ. @) GF=1.0/(EPSI#(EPSI+1.8)) 
3040 IF( H .EB. @) ALPHA=EPS1+2.8 

3078 IF( MH .E8. @) 60 70 108 

3088 DO 20 1=1,2-" 

3098 J=t+ (1-1) 

3108 20 BF=6F/(EPSI+FLOAT (J) ) 

3118 ALPHA=EPS1+2.0 

3128 


31300 COEFFICIENTS FOR SERIES EXPANSION OF THE GAMMA INTEGRAL 


DIFF ",/,) 


31400 SEE HANDBOOK OF MATHEMATICAL FUNCTIONS BY ABRAMOWITZ AND SEGUN 


3150 108 «© C(1)=1. BBBRBORRAROBRABD 
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3168 
3178 
3180 
3198 
3200 
3218 
3220 
3238 
3240 
3258 
3268 
3278 
3288 
3290 
3308 
33108 
3328 
3330 
3340 
3358 
3368 
3378 
3380 
3398 
3488 


C(2)=.5772156649815329 
C(3)=-, 65587887 15282538 
C(4)=-, 0420026358340952 
C(5)=, 1665386113822915 
C(6)=-, 0421977345555443 
C(7)=-. 089621971527887 
C(8)=,087218943246663 
C(9)=-.0011651675918591 
C(10)=-. 0882152416741 149 
C(11)=, 0801280582823882 
C(12)=-, 0008201348547807 
C(13)=-, 0080012504934821 
C(14)=, 008801 1330272320 
C(15)=-. 6880802056338417 
C(16)=6. 116895E-09 
C(17)=5. B02B875E-89 
C(18)=-1, 1812746E-09 
C(19)=1.043427E-18 
C(28)=7, 7B23E-12 
C(21)=-3, 6968BE-12 
C(22)=5, 1E-13 

C(23)=-2. Q6E-14 
C(24)=-5. 4E-15 
C(25)=1.4E-15 


3410 SUH SERIES 


3420 
3430 
3440 


3458 58 


3460 208 


3478 
3480 
3498 
3588 
3518 
3520 
3538 
3548 
35508 
3568 
3578 
3580 
3598 
3688 
3618 
3620 
3630 
3648 
3658 
3668 
3678 
3680 


SUN=0.0 

DO 58 K=1,25 
SUN=SUM#C (K) #(ALPHA®#K) 
AREA=GF/SUN 

RETURN 

END 


SUBROUTINE EXPCT(W,V2,CUT,ALPHA, BETA, LAMBDA) 
DIMENSION ALPHA(3) ,BETA(3) 

DOUBLE PRECISION BU 

REAL LAMBDA 

FDA (X)=-(ALOG(-ALOB(X) ) #PHI) +EPSI 
FD2(X)=((-ALOG(1-X) )## (1/01) )#B1 
FD3(X)=B2/ ( (-ALOG(X) ) ##(1/A2)) 

PDF 4 (X) =EXP(-EXP (-(X-EPSI) /PHI) ) #EXP (-(X-EPSI) /PHI) /PHI 
PDF2(X)=Ala(X##(Al~1)) #EXP(~(X/B1) ##A1) / (BL FEAL) 
PDF3(X) =A2#(B2##A2) #EXP(~(B2/X) ##A2) / (X##(AD+L)) 
CDF (X) =EXP(-EXP( (EPSI-X) /PHI) } 
COF2(X)=1.0-EXP(-(X/B1) ##AL) 
CDF3(X) =EXP (-(B2/X) ##A2) 
B(X)=Wa(1.0-EXP(V28(X-CUT) )) 

PRINT, ‘SELECT A DISTRIBUTION... ' 

PRINT, ‘EXTREMAL TYPE I...1' 

PRINT, ‘WEIBULL... esseceee2’ 

PRINT, 'LOG-EXTREMAL......3' 

PRINT, SELECT 1, 2, OR 3° 

READ, 1D 


Ci2 


IF(ID .LT. 1 .OR. ID .GT. 3) 60 TO 72 


EPSI=ALPHA(1) 

PHI=BETA(1) 

AL=ALPHA (2) 

B1=BETA(2) 

A2=ALPHA(3) 

B2=BETA(3) 

BL=CUT 

IF(ID .EQ. 1) BUSFDI(1-. 0000001) 
IF(ID .EQ. 2) BUSFD2(1-. 0000001) 
IF(ID .EQ. 3) BUSFD3(1-. 0800001) 
SUN=0 

D=BU-BL 

K=-{ 

DO 10 I=1,101 

k=-K 

IF( K .LT. 8) C4 

IF( K .6T. @ )C=2 

IF( 1 EQ. 1 .OR, I EQ. 101)C=1 
ADD=FLOAT (1-1) #D/100.8 
XV=BL+ADD 


IF(ID .€0. 1 .AND. EXP(-(XV-EPSI)/PHI) .67. 82.0) 60 TO 10 
IF(ID .EQ. 2 .AND. ((XV/BL)##A1) .67. 82.8) GO TO 18 
IF(ID .E@. 3 .AND. ((B2/XV)##A2) .GT. 82.8) GO TO 10 


IF( 1D .EQ. 1) FACL=PDFI (XY) 
IF( ID .E8. 2) FACL=PDF2(XV) 
IF( 1D .€8. 3) FACL=PDFS(XV) 
FAC2=6 (XV) 

IF(XV LT. CUT) FAC2=0 
SUR=SUM+C#F ACI #FAC2 
CONTINUE 

SUR=SUM/ 308. 8D 

WRITE (6,14) SUN#LAMBDA 


FORMAT(/,1X,°EXPECTED ANNUAL LOSS IN MILLIONS OF DOLLARS 15°,F14.7) 


RETURN 
END 


C13 


woth oy 
atten Negi a 


ove 
1 | 
Be ee ean ve a 
Bites dhe, scene 
® Asia lom ah cs V7 shia 
; Ane ru “ithe Aneeh : ms an - 
Bisa 6¢) wa ed - y : 
eit ign W atin e 6 | i ? 5 
o ee a = Aa - “js, ot Leung Ak ay aa 
peer act 7 a te pre ab acs te gay 
as n aidan TP Te, Ayah 
* ieee gates 
it 8g, ' 
i ee be 2 
7 ae ; ipl vet Ai oe 
= 4 ee 
af i RT ts i 7 - a ae 4 
af Vii tet he , * “9 ik ian g ¥ Wey 7 
a 2 i : } ab Pt ' - 7 Ts Ff 
or, 4 ¢ wh 1 : = 4 ae cy alfa Netz 
a” : . - oe { : pm 4 4G wi: i. 4 ng, 4 Ris 4 
eas: a rt aa ha de UT an tei | 
Ri ae ! ak - —_ ine deals 
aa | m5 j Oe  s ; 7 — : 
: - ; Marat ie hs ace, “| Tete ease tat ane aR 
= av 1 i ee ee Me ee We? 
0 ; ; FAA RY Tht fee x ls COO PORN CT hy @ 
7 7 : : e Ay if th a eis ? r 
ty 7 i, ay Fa ik ae i 
; 1 4 = : m Wel hoe! my si as f : 
1G fv he 7 
- Facade ewan vies. is 
(aa: he Pie deena 
~ ~" ; ee ) 
; : ; hg Sk iba ; é 
: TAL Sah he Fi aT GY Eh We 
i ty 4 ub DOES Acie BR eee GASES | Peat oa Bi AASbe oh 
a 7 ? * oe ets a 
- : _ - i 
: : : 7 
7 : 7 y 
* _ = 
a - -s - : eS 
7 - L . 7 - . 3 
. - aS 2 : a 
rr : - : a 
t 
= : - | — 
a e : : ese ' 
; ; rm a a : 
eo - ay 
- : = 1" 
a - - _ - ‘“ ser = 1 
7 a = — ' A 
rT’ ‘ - , : a 


APPENDIX D: COMPUTER PROGRAM BWDAMAGE 


Estimation of Rubble-Mound Breakwater 
Armor Layer Expected Damages 


Program purposes 
1. The computer program BWDAMAGE estimates the expected annual damage 


to a rubble-mound breakwater, both in cost and percentage of the armor layer 
displaced. It also estimates the interval that repairs of a specified damage 
level could be scheduled. This information is vital to breakwater optimiza- 
tion analyses which require an estimate of expected repair costs to fully de- 
fine the life cycle costs of the structure. 
Program capabilities 

2. The program assumes that the damage to a particular breakwater con- 
figuration is primarily a function of the type of armor unit that has been 
placed. The breakwater of interest has been designed to suffer a minimum dis- 
placement (erosion) of the primary armor during a design event represented by 
a Significant wave height Hy, . The level of damage acceptable at this level 
of incident wave energy is in current practice defined by the resolution of 
scale model tests and is on the order of 1-5 percent. BWDAMAGE assumes that 
the damage caused by more severe incident events represented by a significant 


wave height H, is. predictable by a function of the following form: 
H [s_.ca/H,-1)] 
BD\ g— } = D(H, )e (D1) 
Hy d 


where S, is the "reserve stability factor." Sh ‘is a relative measure of 
the rate at which breakwaters built with a particular type of armor unit will 
suffer damage with increasing wave heights. A higher Sp value implies a 
higher rate of damage. The program applies the “D(H, = Hq) and S, values 
presented in Table 3 of the main text. The values in Table 3 were determined 
with monochromatic scale model tests which made no account of rocking and the 
probable associated armor unit breakage. The predictions using the above for- 
mula should not be used in final design decisions without careful verification 
in laboratory tests with irregular waves and attention to armor unit rocking. 
The armor units in the program's library include quarrystone (rough, uniform), 
quadripods, tribars, and dolosse. The program also allows input of an armor 


unit title, @D(Hy) , and S, for any other type of unit. 
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3. The expected annual, or long-term average annual, damages are esti- 
mated by the following formula, given the parameters which define the incident 
wave climate by an annual cumulative probability distribution of significant 


wave heights F(H.) as follows: 


5 
H_\[arc) 
BD \ _ eS Nao Se 
ae : vie wl? | Hi dH, (D2) 


where 
E{%D/yr} = the expected annual damage 


A = the Poisson parameter or average number of extreme 
events per year 


(dF(H,)/dH,] = the probability density function corresponding to 
F(H.) 
s 


4, The relation above implicitly assumes that Hy is a significant 
wave height representing some design sea state since its probability of ex- 
ceedance in any year and that of all higher values of H, in bD(H./H 4) are 
determined by a distribution of significant wave heights. It further assumes 
that the number of extreme events per year is a random variable which can be 
represented by a mean value and is independent of the individual significant 
wave heights representing the intensity of these storms. BWDAMAGE requires 
the user to input A and the ce and 9 parameters of the Extremal Type I 


cumulative probability distribution of incident significant wave heights where 


[(e-H,)/o] 
F(H,) =e (D3) 


The associated probability density function is thus 


ae ; P(H.) laa fad 


dH 


(D4) 


Input of the design significant wave height Hy and dA the average number of 
extreme events per year A are required so the expected annual damages can be 


estimated as 


H 
25 H dF(H_) 
BD\ . eS 5 | Seas 
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where Hoe = the significant wave height with an annual probability of 
exceedance of 0.0000001 computed by the program using F(H,). The program 
also requires that a representative repair cost per unit volume of armor layer 
and the total volume of the armor layer per unit trunk length by input. This 


allows the expected annual cost of repairs E{$D/yr} to be estimated as 


ef 8 = vol aa efit (D6) 


5. The interval at which a specified level of damage will occur is 
estimated in two ways. The first simply divides the specified level of per- 
cent damage by the expected annual amount to give an average repair interval 
in years. This is considered to be much more appropriate since it includes an 
account of all the events addressed in the computation of expectation and in- 
directly measures the real world accumulation of damage from successive 
storms. The other method involves solving the @D(H,/Hq) function for the 
wave height H, which would cause the specified level of damage. The speci- 
fied Extremal Type I PCH) is then applied to determine the return period 


RT of this H, as 


1 


RT(H,) oi x[T = FD] (D7) 


This second method is much less conservative and will predict an interval on 


the order of A times as long as that predicted by the first method. 


Program input 
6. Written in FORTRAN IV as implemented on the US Army Engineer Water- 


ways Experiment Station (WES) Honeywell DPS-8 mainframe system, BWDAMAGE is 
completely interactive in its present form. A BASIC version written for the 


IBM PC is also available. 


Sample Interactive Session 


7. The following sample interactive session demonstrates the required 


and optional user input: 
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PROGRAM BWDAMAGE 12/85 VERSION 


POUUU TCT TOC TITIES TTT TTT TT TS SS ST Te of 
"BWDAMAGE” ESTIMATES THE EXPECTED ANNUAL DAMAGE 
TO THE ARMOR LAYER OF A RUBBLEMOUND BREAKWATER, 
BOTH IN COST AND PERCENTAGE DISPLAYED. IT ALSO 
ESTIMATES THE INTERVAL THAT REPAIRS OF A SPECI- 
FIED DAMAGE LEVEL COULD BE SCHEDULED. 

POUPTE TTC CT ICT IIIT P PSST PPT ST CT SS TEE TS So 


* 


Row Oe 
ee 


ENGLISH OR METRIC UNITS(E/M)? 
=E 


«2. QUARRYSTONE (NON-BREAKING WAVES) 
. 2» GUARRYSTONE (BREAKING WAVES) 
...eQUADRIPODS (NON-BREAKING WAVES) 
4....TRIBARS (NON-BREAKING WAVES) 
5...-DOLOSSE (NON-BREAKING WAVES) 
6....DOLOSSE (BREAKING WAVES) 
7...OTHER 

SELECT NUMBER OF ARMOR UNIT 

at 


1 
2 
3 


INPUT HEIGHT OF DESIGN WAVE IN FEET 
=15.0 


SELECT A DISTRIBUTION... 
EXTREMAL TYPE I...1 
WEIBULL... + eeeeeee 2 
LOG-EXTREMAL......3 
SELECT 1, 2, OR 3 

=1 


INPUT EXTREMAL TYPE I EPSILON, AND PHI 
=-2.27,3,216 


INPUT AVERAGE NUMBER OF EXTREMAL EVENTS PER YEAR, 
THE POISSON ‘LAMBDA’ PARAMETER 
=4 


INPUT THE VOLUME OF ARMOR LAYER IN CUBIC FEET PER LINEAR FOOT 
=321.3 


INPUT THE COST OF THE ARMOR IN DOLLARS PER CUBIC FOOT 
=3.98 


1....QUARRYSTONE (NON-BREAKING WAVES) 


EXPECTED DAMAGE PER LINEAR FOOT PER YEAR IS @.42% 
EXPECTED REPAIR COST PER LINEAR FOOT PER YEAR IS $ Sino L 
DO YOU WANT TO PREDICT A REPAIR INTERVAL(Y OR N)? 

inPur % DAMAGE TO ARMOR LAYER AT TIME OF REPAIRS 


THE AVERAGE REPAIR INTERVAL FOR 5.008% DAMAGE 
BASED ON @.42% PER YEAR EXPECTED DAMAGE IS 11.8 YEARS 


THE RETURN PERIOD OF THE STORM CAUSING 5.08 % DAMAGE IS 75.81 YEARS 


ANOTHER RUN (Y/N)? 
=N 
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Sample Program Listing 


8. A sample program listing for BWDAMAGE (FORTRAN version) is given below: 


10C PROGRAM "BWDAMAGE” 12/85 VERSION 

20C DESIGN BRANCH-COASTAL ENGINEERING RESEARCH CENTER 
3@C U.S. ARMY ENGINEERS WATERWAYS EXPERIMENT STATION 
40C VICKSBURG, MS 39180-0631 

5@C FOR FURTHER INFORMATION CONCERNING THE APPLICATION 
60C OF "BWDAMAGE", CALL.» ceeeoees 


70C ORSON P. SMITH 601-634-2013 FTS:542-2013 OR 
8@C ROBERT B. LUND 601-634-2068 FTS:542-2068 OR 
90C DOYLE L. JONES 681-634-2069 FTS:542-2069 

100 

110C FORTRAN 4 HONEYWELL DPS-8 


120C REF: "COMPUTER PROGRAM WAVDIST" CETN-I- 

130C REF: “EXTREMAL STATISTICS IN WAVE CLIMATOLOGY” BY BORGMAN AND RESIO 

14@C REF: "SHORE PROTECTION MANUAL” 1984 4TH ED., 2 VOLS. 

150C REF: "COST EFFECTIVE OPTIMIZATION OF RUBBLEMOUND BREAKWATER" BY O. SMITH 
160 


170 

188C All = REGRESSION COEFFICIENT FOR DAMAGE EQUATION 
198C A22 = REGRESSION COEFFICIENT FOR DAMAGE EQUATION 
208C EPSI = EXTREMAL TYPE 1 LOCATION PARAMETER 

2108C PHI = EXTREMAL TYPE 1 SCALE PARAMETER 

226C Al = WEIBULL LOCATION PARAMETER 

230C Bl = WEIBULL SCALE PARAMETER 

240C A2 = LOG-EXTREMAL LOCATION PARAMETER 

250C B2 = LOG-EXTREMAL SCALE PARAMETER 

260C ID = IDENTIFIES DISTRIBUTIONS 


270C LAMBDA= POISSON ‘LAMBDA’ PARAMETER 


280C BL = HEIGHT OF DESIGN WAVE 

290C BU = _ UPPER LIMIT OF INTEGRATION 

300C L = £ OR M, ENGLISH OR METRIC UNITS 

310C K = NUMBER OF ARMOR UNIT USED 

320C V1 = VOLUME OF ARMOR LAYER PER LINEAR FOOT (OR METER) 
330C V2 = COST OF ARMOR IN DOLLARS PER CUBIC FEET (OR METERS) 
340C DP = PERCENTAGE OF DAMAGE TO ARMOR LAYER AT TIME OF REPAIRS 
350C 

360C INITIALIZE VARIABLES,STRINGS,AND FUNCTIONS 

370 DIMENSION A11(7),A22(7) 

388 COMMON EPSI,PHI,A1,B1,A2,B2,1D,LAMBDA 

392 REAL LAMBDA 

400 CHARACTER#1 L,ANS 

410 CHARACTER#6@ SL(9) 

420 CHARACTER*37 DUM,ST(7) 

430 CHARACTER#*6 UNIT (2) 

440 CHARACTER#12 VOL(2) 

450 CHARACTER#5 LEN(2) 

460 DOUBLE PRECISION BL,BU 

470 F(X) =TL#EXP(T2#(X-1)) 

480 DARAVALI/Sce28 Se ySiates enter 

490 DATA A22/6.95,3.65,6.,4.87,1.68,3.55,0./ 

500 UNIT(1)= ‘METERS’ 

510 UNIT (2)= ‘FEET 
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520 VOL(1)2‘CUBIC METERS’ 


538 VOL(2)=‘CUBIC FEET’ 

548 LEN(1)='METER’ 

558 LEN(2)=‘FOOT’ 

568 SL (1) = 1 ETHIE HHH HEHEHE TE TAHA EEE HEE SEHEHEHEE E E E E E 
578 SL(2)='* "BWDAMAGE" ESTIMATES THE EXPECTED ANNUAL DAMAGE * 
588 SL(3)='* TO THE ARMOR LAYER OF A RUBBLEMOUND BREAKWATER, * 
598 SL(4)='* BOTH IN COST AND PERCENTAGE DISPLAYED. IT ALSO +* 
600 SL(5)='* ESTIMATES THE INTERVAL THAT REPAIRS OF A SPECI- * 
618 SL(6)='* FIED DAMAGE LEVEL COULD BE SCHEDULED. * 
628 SL(7)=SL(1) 

638 SL(8)=' 

640 SL(9) = 

650 DUM='INPUT THE VOLUME OF ARMOR LAYER IN 

660 ST(1)='1....-QUARRYSTONE (NON-BREAKING WAVES) ’ 

670 ST(2)='2....-QUARRYSTONE (BREAKING WAVES) ' 

688 §T(3)='3....QUADRIPODS (NON-BREAKING WAVES) ‘ 

690 ST(4)='4,...TRIBARS (NON-BREAKING WAVES) ‘ 

708 ST(5)='5....DOLOSSE (NON-BREAKING WAVES) ‘ 

710 ST(6)='6...-DOLOSSE (BREAKING WAVES) ' 

720 ST(7)='7...0THER’ 

730 WRITE (6,43) 

748 43 FORMAT(//,1X,"PROGRAM BWDAMAGE 12/85 VERSION",///) 
750 1 IuU=1 

760 DO 5@ 1=1,9 

770 WRITE(6,52) SL(I) 

788 52 FORMAT (1X, A6Q) 

798 5 CONTINUE 

800 9 PRINT, ‘ENGLISH OR METRIC UNITS(E/M)?° 

818 READ(S,7) L 

8208 7 FORMAT (A1) 

838 IF(L .NE. ‘E’ .AND. L .NE. ‘M’) GO TO 9 

848 IF(L .EQ@. 'E’ ) TUs2 

850 PRINT,’ 

868 

870C SELECT ARMOR UNIT 

880 DO 53 I=1,7 

898 WRITE(6,4) ST(1) 

900 4 FORMAT (1X ,A37) 

918 53 CONTINUE 

928 3 PRINT, SELECT NUMBER OF ARMOR UNIT’ 

930 READ,K 

948 EEG K ali: 1 JOR K wGits, 7) GOTO S 

958 IF (K .NE. 7) GO TO 41 

960 PRINT, INPUT NAME OF ARMOR UNIT’ 

978 READ(S,71) ST(7) 

988 71 FORMAT (A2@) 

998 PRINT, INPUT COEFFICIENTS FOR DAMAGE EQUATION-' 

10008 PRINT, %D(Hd) AND Sr FOR ZD(H/Hd) = %D(Hd)explSr(H/Hd - 1)] 
1018 READ, A11(7),A22(7) 

1020 

1038 41 PRINT,’ 

1048 PRINT, ‘INPUT HEIGHT OF DESIGN WAVE IN ‘,UNIT(IU) 
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1058 READ, BL 


1060 IF(BL .LE. ® ) GO TO 41 

1070 PRINT,’ ’ 

1080C INITIALIZE VARIABLES FOR INTEGRATION 

1090 18 TLEALL(K) 

1100 T2=A22(K) 

1110 CALL EXPCT(T1,12,BL,BU,SUM,D) 

1120 PRINT, ' 

1130 

1140C CONVERT FRACTIONAL DAMAGE TO DOLLARS DAMAGE 

1150 IF(IU .EQ. 1)PRINT, ‘INPUT THE VOLUME OF ARMOR LAYER IN CUBIC METERS 
1160 & PER LINEAR METER‘ 

1170 IF(IU .EQ. 2)PRINT, ‘INPUT THE VOLUME OF ARMOR LAYER IN CUBIC FEET 
1180 & PER LINEAR FOOT’ 

1190 READ, V1 

1208 PRINT, ’ 

1210 PRINT, ‘INPUT THE COST OF THE ARMOR IN DOLLARS PER CUBIC ‘,LEN(IU) 
1220 READ, V2 

1230 SUM=SUM#V1#V2 

1240 WRITE(6,30) ST(K) ,LEN(IU) ,D 


1250 38 FORMAT(//,1X,A37,/,1X,"EXPECTED DAMAGE PER LINEAR ",AS," PER 
1260 & YEAR IS",F6.2,"%") 


1270 WRITE (6,32) LEN(IU),SUM 
1288 32 FORMAT(/,1X,"EXPECTED REPAIR COST PER LINEAR ",AS," PER YEAR IS $",F8.2,///) 
1290 IF(SUM .LE. ®) GO TO 31 

1300 PRINT, ‘DO YOU WANT TO PREDICT A REPAIR INTERVAL(Y OR N)?‘ 

1310 READ(5,33) ANS 

1328 33 FORMAT (A1) 

1330 IF(ANS.EQ.‘'N’) GO TO 31 

1340 93 PRINT, INPUT % DAMAGE TO ARMOR LAYER AT TIME OF REPAIRS’ 

1350 READ, DP 

13468 IF(DP .LT. 0.0 .OR. DP .GT. 100.8) GO TO 83 

1370 HDP=BL#(1+(ALOG(DP/T1)/T2)) 

1380 RT=0.0 

1390 IF(ID .EQ. 1 .AND. (HDP-EPSI)/PHI .GT. 82.8 ) GO TO 247 

1400 IF(ID .EQ. 1 .AND. EXP(-(HDP-EPSI)/PHI) .GT. 82.0 ) GO TO 247 
1410 IF(ID .EQ. 2 .AND. ((HDP/B1)##A1).GT. 82.@ ) GO TO 247 

1420 IF(ID .EQ. 3 .AND. (-(B2/HDP)##A2) .GT. 82.8 ) GO TO 247 

1430 IF(ID .EQ. 1) PHDP=EXP(-EXP(-(HDP-EPSI) /PHI)) 

1440 IF(ID .EQ. 2) PHDP=1.0-EXP(-(HDP/B1) ##A1) 

1450 IF(ID .EQ. 3) PHDP=EXP(-(B2/HDP) ##A2) 

1460 IF( PHDP .GT. .9999) GO TO 247 

1470 RT=1./(LAMBDA#(1.-PHDP)) 

1480 247 REPAIR=DP/D 

1490 WRITE (6,36) DP,D,REPAIR 

1500 36 FORMAT(//,1X, ‘THE AVERAGE REPAIR INTERVAL FOR ,F.2,'% DAMAGE’, 
1510 & /,‘BASED ON ',F6.2,‘% PER YEAR EXPECTED DAMAGE IS ’,F5.1,' YEARS’) 
1520 IF( PHDP .GT. .9999) GO TO 31 

1530 IFC RT .LE. af) GO TO 31 

1540 WRITE (6,37) DP,RT 

1558 37 FORMAT(/,1X,‘THE RETURN PERIOD OF THE STORM CAUSING ',F6.2,' % 
1560 & DAMAGE IS',F7.2,' YEARS’,/) 

1570 31 WRITE (6,321) 
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1588 321 FORMAT(///,1X,"ANOTHER RUN (Y/N) ?") 


1590 READ(5,35) L 

1600 35 FORMAT(AL) 

1610 TFL EQ. *¥°) G0, 104 

1620 IF(L .NE. 'N’) GO TO 31 

1630 STOP 

1640 END 

1658 

1660 SUBROUTINE EXPCT(T1,12,BL,BU,SUM,D2) 

1670 COMMON EPSI,PHI,A1,B1,A2,B2,1D,LAMBDA 

1688 DOUBLE PRECISION BL,BU 

1698 REAL LAMBDA 

1700 FD1(X)=-(ALOG(-ALOG(X)) *PHI) +EPSI 

1710 FD2(X) =((-ALOG(1-X)) ## (1/01) ) #B1 

1728 FD3(X) =B2/( (-ALOG(X)) ##(1/A2)) 

1730 PDF1(X) =EXP (-EXP (-(X-EPSI) /PHI) ) EXP (-(X-EPSI) /PHI) /PHI 
1748 PDF2(X) =AL#(X##(AL-1)) #EXP(-(X/B1) ##A1) /(BL##AL) 
1750 PDF3(X) =A2* (B2##AZ) EXP (-(B2/X) #¥AZ) / (X## (A241) 
1768 CDF1(X) SEXP (-EXP( (EPSI-X) /PHI)) 

1778 CDF2(X)=1.Q-EXP(-(X/B1) ##AL) 

1780 CDF3(X) =EXP (-(B2/X) #*A2) 

1798 G(X) =TL#EXP(T2*(X-1)) 

1800 70 PRINT, SELECT A DISTRIBUTION... ’ 

1810 PRINT, ‘EXTREMAL TYPE I...1' 

1828 PRINT, ‘WEIBULL. ...easeeee a 

1830 PRINT, ‘LOG-EXTREMAL......3° 

1840 PRINT, ‘SELECT 1, 2, OR 3° 

1850 READ, ID 

1860 IF(ID .LT. 1 .OR. ID .GT. 3) GO TO 70 

1870 IF(ID .EQ. 1) WRITE(6,104) 

1880 104 FORMAT(/,1X,"INPUT EXTREMAL TYPE I EPSILON, AND PHI ") 
1898 IF(ID .EQ. 1) READ,EPSI,PHI 

1990 IF( ID .EQ. 2) WRITE(6,114) 

1910 114 FORMAT(/,1X,"INPUT WEIBULL ALPHA AND BETA" 

1928 IF( ID .EQ. 2) READ,A1,Bi 

1930 IF( ID .EQ. 3) WRITE(6,124) 

1948 124 FORMAT(/,1X,"INPUT LOG-EXTREMAL ALPHA AND BETA" 
1950 IF( ID .EQ. 3) READ,A2,B2 

1960 WRITE(4,5) 

1978 5 FORMAT(/,1X,"INPUT AVERAGE NUMBER OF EXTREMAL EVENTS PER YEAR,", 
1980 & /,1X,"THE POISSON ‘LAMBDA’ PARAMETER") 

1990 READ, LAMBDA 

2000 IF(ID .EQ. 1) BU=FD1(1-. 0000001) 

2010 IF(ID .EQ. 2) BU=FD2(1-.0002001) 

2028 IF(ID .EQ. 3) BU=FD3(1-.0000001) 

2030 SUM=0 

2040 D=BU-BL 

2050 IF( DLE. 8) 60 TO 11 

2068 K=-1 

2080 DO 10 I=1,101 

2090 K=-K 

2108 IF( K .LT. @) C4 

2110 IF( K .6T. @ )C=2 

2120 IF( 1 .EQ. 1 .0OR. I .EQ. 101)C=1 

2130 ADD=FLOAT(I-1)#D/100.9 

2140 XV=BL+ADD 

2150 IF(ID .EQ. 1 .AND. EXP(-(XV-EPSI)/PHI) .GT. 82.8) GO TO 10 
2162 IF(ID .EQ. 2 .AND. ((XV/B1)##A1) .GT. 82.8) GO TO 10 
2178 IF(ID .EQ. 3 .AND. (-(B2/XV)##A2) .GT. 82.0) GO TO 18 
2180 IF( ID .EQ. 1 ) FACL=PDF1(XV) 

2190 IF( ID .EQ. 2.) FACL=PDF2(XV) 

2200 IF( ID .€Q. 3.) FACL=PDF3(XKV) 

2210 FAC2=G(XV/BL) /100.@ 

2220 IF( FAC2 .LT. @ ) FAC2=0 

2230 IF( FAC2 .GT. 1 ) FAC2=1 

2240 SUM=SUM+C#FACL¥FAC2 

2260 18 CONTINUE 

2270 11 SUM=SUM/300.Q*D#LAMBDA 

2280 D2=SUM*100.0 

2290 RETURN 

2300 END 
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APPENDIX E: NOTATION 


Empirical runup coefficient 
Empirical runup coefficient 


Empirical coefficient for prediction of economic losses 
due to wave attack 


Average eroded cross sectional breakwater area 

Empirical runup coefficient 

Breakwater crest width 

Empirical runup coefficient 

A parameter of the Weibull probability distribution 
Empirical coefficient for wave transmission by overtopping 
Dimensionless breakwater damage 

Breakwater damage rate (armor displacement per unit time) 


Sediment size from a graded sample for which 15 percent is finer 
by weight 


Sediment size from a graded sample for which 85 percent is finer 
by weight 


"Nominal" or equivalent cube dimension of the mean stone weight 
of a graded sample based on size gradation analysis 


Mathematical expectation 

Statistical probability density function of x 
Freeboard 

Cumulative probability density function of x 
Acceleration of gravity 

Eroded breakwater crest height 

Original breakwater crest height 

Wave height 

Critical wave height regarding wave transmission 
Design wave for armor unit sizing 

Incident wave height 


Maximum wave height below which economic losses can be neglected 
due to wave attack 


Significant wave height 
Transmitted wave height 


Empirical armor stability coefficient, as applied in the Hudson 
formula 


E1 


Armor stability coefficient, as applied in the DHI-Iribarren 
formula 


Coefficient of wave transmission by overtopping 

Layer coefficient 

Wave length corresponding to the period of peak energy density 
Number of armor units comprising the layer thickness 


Armor stability number, as applied in the original Iribarren 
formula 


Number of waves in a specific model test 
Layer thickness 

Statistical correlation coefficient 
Runup height above still-water level 
Reserve stability coefficient 


Dimensionless damage, equivalent to the number of Dn50 cubes 
eroded from a cross sectional width of Dn50 


Time duration of wave exposure 
Wave period of peak energy density 


Average wave period, derived from a count of "zero crossings" 
with respect to an arbitrary mean 


Wave period 

Flow velocity 

Total armor layer volume per unit breakwater trunk length 
Weight of a primary armor unit 

The mean stone weight of a graded sample 

Armor material density relative to sea water = (0, - pi )/0., 


A parameter of the Extremal Type I and Weibull probability 
distributions 


Angle from horizontal of the seaward slope of a breakwater 


Poisson parameter, denoting the mean number of extreme storms per 
year 


Coefficient of static friction 
Kinematic fluid viscosity 

Surf similarity parameter 

Mass density of the armor material 
Mass density of sea water 


Wave directional spreading parameter 
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$Bret 
$Benefits 


$D 
¢1St 
$L 
$L' 


$Linax 
$Total 


%D 
D* 
$/Vol 


Natural angle of repose of armor material 
Wave angle 


A parameter of the Extremal Type I and Weibull probability 
distribution 


Net incremental benefits 

Total incremental benefits 

Cost of breakwater repair 

Construction cost 

Economic losses due to wave attack 

Economic losses due to attack by transmitted waves 
Maximum conceivable economic losses due to wave attack 


Total incremental project feature (breakwater) cost, including 
associated economic losses 


Percent of breakwater armor displaced 
A critical value of breakwater damage by displacement 


Cost per unit volume per unit breakwater trunk length 
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